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Preface 


This quarterly publication provides archival reports on developments in programs 
managed by JPL’s Office of Telecommunications and Data Acquisition (TDA). In space 
communications, radio navigation, radio science, and ground-based radio astronomy, it 
reports on activities of the Deep Space Network (DSN) and its associated Ground Com- 
munications Facility (GCF) in planning, in supporting research and technology, in imple- 
mentation, and in operations. Also included is TDA-funded activity at JPL on data and 
information systems and reimbursable DSN work performed for other space agencies 
through NASA. The preceding work is all performed for NASA’s Office of Space Track- 
ing and Data Systems (OSTDS). 

In geodynamics, the publication reports on the application of radio interferometry 
at microwave frequencies for geodynamic measurements. In the search for extraterrestrial 
intelligence (SETI), it reports on implementation and operations for searching the micro- 
wave spectrum. The latter two programs are performed for NASA’s Office of Space 
Science and Applications (OSSA). 

Finally, tasks funded under the JPL Director’s Discretionary Fund and the Caltech 
President’s Fund which involve the TDA Office are included. 

This and each succeeding issue of the TDA Progress Report will present material in 
some, but not necessarily all, of the following categories: 

OSTDS Tasks: 

DSN Advanced Systems 

Tracking and Ground-Based Navigation 
Communications, Spacecraft-Ground 
Station Control and System Technology 
Network Data Processing and Productivity 
DSN Systems Implementation 
Capabilities for New Projects 
Networks Consolidation Program 
New Initiatives 
Network Sustaining 
DSN Operations 

Network Operations and Operations Support 
Mission Interface and Support 
TDA Program Management and Analysis 
GCF Implementation and Operations 
Data and Information Systems 

OSSA Tasks: 

Search for Extraterrestrial Intelligence 
Geodynamics 

Geodetic Instrument Development 
Geodynamic Science 


Discretionary Funded Tasks 
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VLBI measurements of time delay and fringe frequency at 2.29 GHz on baselines of 
10 4 km between Deep Space Network stations have been used to determine the positions 
of the milliarcsecond nuclei in 207 extragalactic radio sources. Estimated accuracies gen- 
erally range from ~0".l to ~1".0, in both right ascension and declination, with all sources 
having uncertainties < 4" in both coordinates. The observed sources are part of an all-sky 
VLBI catalog of milliarcsecond radio sources. Arcsecond positions have now been deter- 
mined for 752 of these sources. Arcsecond positions serve as a useful starting point in the 
construction of high-precision VLBI reference frames and are also important for unambig- 
uous determination of optical counterparts to compact radio sources. 


I. Introduction 

A survey is underway to develop an all-sky catalog of radio 
sources with milliarcsecond components at 2.29 GHz. This is 
being accomplished by searching for compact components in 
known extragalactic sources with intercontinental VLBI base- 
lines (Refs. 1 and 2). We have used the VLBI measurements of 
time delay and fringe frequency at 2.29 GHz to determine the 
positions of the milliarcsecond nuclei. Previously we deter- 
mined the positions of 546 nuclei (Ref. 3). In this article we 
determine the positions of an additional 206 nuclei and in- 
clude an improved position for one source (M 104) from 
Ref. 3. Estimated accuracies in this article range from 0"1 to 


3'.'8 in both right ascension and declination. The addition of 
the Australia to South Africa baseline has allowed the survey 
to continue south of -45° declination. In total we have deter- 
mined arcsecond positions for 84% of the detected sources, 

Arcsecond positions serve as a useful starting point in the 
construction of a high precision VLBI reference frame as well 
as allow unambiguous optical identifications. The positions 
from this article and Ref. 3 are presently being used as a first 
step in the formation of a precision reference frame of about 
200 sources in which relative radio positions should be deter- 
mined to milliarcsecond accuracy (Ref. 4). 
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II. The Observations 

The observations were performed with pairs of antennas on 
either California-Australia, California-Spain, or Australia-South 
Africa baselines (see Table 1) in 16 separate observing sessions 
between 1978 and 1982. A list of experiments appears in 
Table 2. The observations were performed at 2.29 GHz with 
MK II VLBI recording systems as described in Ref. 3. The 
length of each observation was a few minutes with 38% of the 
sources being observed more than once. 

III. Method of Position Determination 

The details of position determination for this set of sources 
closely follow the analysis given in Ref. 3, except for three 
differences: 

(1) The data for the 207 sources reported here were re- 
duced using the Caltech/JPL VLBI Mark II correlator 
and post-correlation software, whereas the data in 
Ref. 3 were reduced using the NRAO Mark II correla- 
tor and post-correlation software. 

(2) Increased accuracy in modeling of the troposphere and 
earth rotation parameters yielded a smaller scatter in 
the time delay and fringe frequency residual. Typical 
residual values for the rms time delay scatter were 
~10 ns and for the rms fringe frequency scatter 
~0.1 mHz. The observed scatter was assumed to be an 
estimate of the random measurement error. Earth rota- 
tion corrections were obtained from final smoothed 
values of UT1-UTC and polar motion from Bureau 
International de l’Heure (BIH). 

(3) Due to improved modeling of the troposphere and 
earth rotation parameters, 0"3 rather than O'.'S was 
added in a root-sum-square manner to formal estimates 
of right ascension and declination uncertainty in each 
experiment to account for unmodeled effects. 

As in Ref. 3, observations of sources with well-known posi- 
tions allowed instrumental delays and frequency offsets to be 
determined. At least 3 such calibration sources per experiment 
were spread in time among the sources whose positions we 
wished to determine (see Table 2). A list of the 44 calibration 
sources utilized appears in Table 3. Twenty-eight of the cali- 
bration source positions have been determined with VLBI 
(Ref. 5) and can be referred to the FK4 reference frame 
(Ref. 6) with an accuracy of 0”l. Other calibration source 
positions came from 1) Ref. 7, 8 sources, O'.'l accuracy; 
2) Ref. 3, 5 sources, 0'.'5 accuracy, 3) Ref. 8, 2 sources, l" 
accuracy, and 4) D. L. Jauncey et al. (private communication), 
1 source, O'! 1 accuracy. 


IV. Results 

The calculated positions of 207 sources and the corre- 
sponding uncertainties are shown in Fig. 1 . The positions are 
referred to the equinox of 1950.0, and elliptical aberration 
terms are included so as to agree with past astronomical con- 
vention. The source positions and position uncertainties for 
sources which were multiply observed were estimated from a 
weighted average. 

Seventy-eight sources were observed two or more times, 
and the scatters in their position estimates are consistent with 
the estimated uncertainties. Figure 2 displays a histogram of 
the ratio of the weighted rms scatter of the individual position 
coordinate estimates from the mean position of each multiply 
observed source, to the weighted rms estimate of position 
coordinate uncertainty. For the 78 multiply observed sources, 
the rms value of this ratio is 0.4 for right ascension and 0.5 for 
declination. Thus there is excellent agreement between the 
multiple observations of each source. 

In addition to the 207 sources listed in Fig. 1, we also 
reobserved 22 sources whose positions were reported in Ref. 3. 
For these 22 sources, Fig. 3 shows a histogram of the absolute 
value of the difference in each position coordinate between 
the two separate observations divided by the RSS error of the 
position coordinate uncertainties. There is excellent agreement 
between the common sources. The calculated ratio has an 
RMS value of 1 .3 for right ascension and 1 .0 for declination. 
For one source, Ml 04, an erroneous position was given in 
Ref. 3. The position for Ml 04 cited in Fig. 1 should be used 
instead. 

For 25 sources, we could compare our position estimates 
with other position estimates of better or similar accuracy 
(Refs. 7, 9, and 10). Figure 4 displays the ratio of the value of 
the difference between our source position estimate and the 
other catalog value to the RSS of the position uncertainties of 
both catalogs. For these 25 sources, the rms value of this ratio 
is 0.5 for right ascension and 0.9 for declination. For both 
right ascension and declination, the value of the ratio never 
exceeds 3.2. The bias offsets between our position estimates 
relative to those of the other catalogs for these 25 sources are 
0?005 ±0?008 for right ascension and -0"l 1 ±0'.’06 for decli- 
nation. Hence, these position estimate comparisons along with 
the multiple observation comparisons indicate our position 
uncertainty estimates are realistic. 

Figure 5 shows histograms of the number of sources versus 
estimated position uncertainty for declination and right ascen- 
sion, respectively. Estimated accuracies range from ~0'.'l to 
3'.'8 in both right ascension and declination. 


2 



V. Conclusion 

Positions for the milliarcsecond nuclei of 207 extragalactic 
sources have been determined to an accuracy of ~0'.'l to 3"8 
in both right ascension and declination. The reliability of the 


determined positions has been demonstrated by testing the 
repeatability of multiple observations on the same source and 
by comparing the results with other radio catalogs. Arcsecond 
positions have now been determined for 752 milliarcsecond 
nuclei. 
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Table 1. Participating observatories 


Location 

Designation 

Diameter, n 


Baseline Length 


Kilometers 

Wavelengths 

Tidbinbilla, Australia 

DSS 43 

64 








— 10.6 X 10 3 

8.1 X 10 7 

Goldstone, California 

DSS 14 

64 





DSS 13 

26 

_ 

84k in 3 

6.4 X 10 7 

Madrid, Spain 

DSS 63 

64 



Parkes, Australia 

PK 

64 


— 9.7 X 10 3 

7.5 X 10 7 

Hartebeesthoek, So. Africa 

HT 

26 

- 




Table 2. Experiment list 



Experiment 

Number of Position 
Stations Calibrator Source 

Observations 

Number of Source 
Positions Determined 

Yr. 

Mo. 

Day 

78 

05 

30/31 

13/43 

4 

27 

78 

06 

29/30 

13/43 

8 

10 

80 

02 

01 

14/43 

4 

7 

80 

02 

27/28 

14/43 

3 

10 

80 

03 

03 

13/63 

8 

25 

80 

03 

12/13 

13/63 

13 

26 

80 

03 

28 

13/63 

4 

25 

80 

04 

24 

PK/HT 

9 

17 

80 

04 

25 

PK/HT 

11 

14 

80 

04 

25/26 

PK/HT 

10 

14 

80 

06 

19 

14/43 

6 

9 

81 

05 

08 

14/43 

3 

1 

82 

02 

14 

43/HT 

8 

25 

82 

02 

17 

43/HT 

3 

33 

82 

02 

19 

43/HT 

8 

38 

82 

04 

20 

43/HT 

6 

43 






Table 3. List of calibration sources (all positions from Ref. 5 unless 
otherwise noted) 


P0003-00 2 

0016+73 1 

PO 104-408 

P0 106+01 

PO 113-118 

P033 2-403 

P0402-362 

P0438-43 

P0537-441 

DA193 

P0605-08 

0615+82 1 

0716+71 1 

P0723-008 3 

P0727-11 


P0743-006 3 

4C7107 

OJ287 

3C212 3 

4C39.25 

PI 104-445 

P1127-14 

PI 144-379 

1 150+81 1 

3C273 

3C274 3 

P1333-33 3 

3C309.1 

DW1555+00 

P1610-77 4 


3C345 

NRA0530 

P1741-038 

1749+70 1 

1803+78 1 

OV-236 

1928+73 1 

P1933-400 

2007+77 1 

P21 34+004 

P2215+02 2 

OY-172.6 

P2245-328 

P2345-16 


‘ Ref. 7 
2 Ref. 8 
2 Ref. 3 

D. L. Jauncey, et al., private communication 



SOURCE NAME 


RIGHT ASCENSION 


DECLINATION 



HR 

MIN 

SEC 

ERROR 

DEG 

MIN 

SEC 

ERROR 

P 0002-478 

0 

2 

2. 978 

0. 033 

-47 

53 

1. 81 

0. 30 

0014+81 

0 

14 

4. 473 

0. 151 

+81 

18 

28. 71 

0. 22 

0018+72 

0 

18 

34. 482 

0. 107 

+72 

56 

3. 93 

0. 24 

P 0019-00 

0 

19 

51. 665 

0. 020 

- 0 

1 

41. 88 

0. 32 

P 0022-423 

0 

22 

15. 413 

0. 020 

-42 

18 

41. 12 

0. 21 

0027+70 

0 

27 

17. 026 

0. 083 

+70 

21 

6. 26 

0. 31 

P 0034-01 

0 

34 

30. 555 

0. 020 

- 1 

25 

37. 78 

0. 36 

P 0035-02 

0 

35 

47. 176 

0. 020 

- 2 

24 

9. 45 

0. 30 

P 0036+03 

0 

36 

44. 176 

0. 020 

+ 3 

3 

24. 54 

0. 43 

P 0047-579 

0 

47 

48. 238 

0, 064 

-57 

54 

46. 83 

0. 30 

P 0055-01 

0 

55 

1. 572 

0. 020 

- 1 

39 

39. 40 

0. 34 

P 0056-572 

0 

56 

38. 640 

0. 041 

-57 

15 

22. 07 

0. 30 

P 0131-522 

1 

31 

5. 633 

0. 022 

-52 

15 

26. 40 

0. 17 

P 0133-203 

1 

33 

13. 577 

0. 025 

-20 

24 

4. 01 

0. 44 

0149+71 

1 

49 

20. 814 

0. 057 

+71 

0 

20. 76 

0. 24 

B2 0149+33 

1 

49 

40. 027 

0. 024 

+33 

35 

46. 85 

0. 31 

0153+74 

1 

53 

4. 339 

0. 054 

+74 

28 

5. 69 

0. 21 

0159+72 

1 

59 

13. 059 

0. 057 

+72 

18 

29. 24 

0. 23 

P 0202-76 

2 

2 

1. 005 

0. 422 

-76 

34 

26. 34 

0. 30 

0205+72 

2 

5 

26. 907 

0. 050 

+72 

15 

16. 40 

0. 22 

GC 0206+35 

2 

6 

39. 339 

0. 032 

+35 

33 

41. 36 

0. 40 

P 0208-512 

2 

8 

56. 970 

0. 019 

-51 

15 

7. 53 

0. 17 

0212+73 

2 

12 

49. 935 

0. 050 

+73 

35 

40. 15 

0. 21 

P 0214-522 

2 

14 

17. 229 

0. 029 

-52 

14 

5. 49 

0. 21 

P 0219-637 

2 

19 

37. 871 

0. 062 

-63 

43 

59. 84 

0. 21 

P 0220-349 

2 

20 

49. 570 

0. 028 

-34 

55 

4. 67 

0. 35 

GC 0223+34 

2 

23 

9. 735 

0. 027 

+34 

8 

1. 56 

0. 33 

P 0226-038 

2 

26 

21. 984 

0. 021 

- 3 

50 

57. 25 

0. 42 

P 0229-398 

2 

29 

51. 990 

0. 017 

-39 

49 

0. 15 

0. 18 

P 0230-790 

2 

30 

28. 852 

0. 139 

-79 

1 

0. 77 

0. 21 

P 0232-04 

2 

32 

36. 512 

0. 022 

- 4 

15 

8. 90 

0. 65 

P 0235-618 

2 

35 

37. 782 

0. 034 

-61 

49 

13. 47 

0. 21 

P 0244-452 

2 

44 

4. 579 

0. 023 

-45 

12 

13. 34 

0. 21 

GC 0248+43 

2 

48 

18. 504 

0. 027 

+43 

2 

56. 98 

0. 34 

P 0252-549 

2 

52 

0. 273 

0. 023 

-54 

54 

2. 27 

0. 17 

P 0302-623 

3 

2 

48. 181 

0. 034 

-62 

23 

3. 71 

0. 21 

P 0308-611 

3 

8 

51. 298 

0. 027 

-61 

9 

58. 24 

0. 17 

P 0312-77 

3 

12 

56. 328 

0. 205 

-77 

3 

0. 15 

0. 30 

P 0316-444 

3 

16 

13. 353 

0. 056 

-44 

25 

11. 24 

0. 30 

DW 0326+27 

3 

26 

56. 037 

0. 026 

+27 

46 

0. 31 

0. 33 

P 0332-403 

3 

32 

25. 244 

0. 026 

-40 

18 

24. 04 

0. 30 

P 0334-546 

3 

34 

36. 057 

0. 029 

-54 

40 

16. 95 

0. 21 

OE 367 

3 

40 

14. 804 

0. 029 

+36 

12 

44. 47 

0. 42 

P 0355-66 

3 

55 

27. 892 

0. 095 

-66 

54 

11. 43 

0. 22 

P 0355-483 

3 

55 

52. 505 

0. 027 

-48 

20 

49. 14 

0. 21 

GC 0402+37 

4 

2 

29. 854 

0. 029. 

+37 

55 

26. 85 

0. 42 

0403+76 

4 

3 

59. 196 

0. 067 

+76 

48 

52. 75 

0. 25 

P 0414-341 

4 

14 

16. 908 

0. 028 

-34 

10 

26. 06 

0. 37 

3C 119 

4 

29 

7. 937 

0. 027 

+41 

32 

8. 24 

0. 31 

P 0431-512 

4 

31 

4. 407 

0. 043 

-51 

15 

42. 11 

0. 23 

P 0437-454 

4 

37 

30. 667 

0. 025 

-45 

28 

12. 26 

0. 22 

P 0446-519 

4 

46 

35. 030 

0. 049 

-51 

56 

13. 75 

0. 23 

P 0448-392 

4 

48 

0. 452 

0. 017 

-39 

16 

15. 74 

0. 18 

P 0450-469 

4 

50 

27. 871 

0. 024 

-46 

58 

16. 45 

0. 21 

P 0454+039 

4 

54 

8. 935 

0. 020 

+ 3 

56 

14. 71 

0. 41 

P 0454-81 

4 

54 

18. 188 

0. 176 

-81 

5 

54. 13 

0. 21 
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Fig. 3. Comparison of position estimates with 22 sources common 
with Ref. 3. For each source, A is the ratio of the absolute value of the 
difference between the source coordinate estimate in this set of data 
and the source coordinate estimate cited in Ref. 3 to the weighted 
rms estimate of the uncertainty in that position coordinate. 


Fig. 2. Comparison of position estimates for 78 multiply observed 
sources. For each source, S is the ratio of the weighted rms scatter 
of individual position coordinate estimates about the mean value to 
the weighted rms estimate of the uncertainty in that position 
coordinate. 
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A 

Fig. 4. Comparison of position estimates with other radio positional 
catalogs. For each source, A is the ratio of the difference between 
our source coordinate estimate and the other catalog value to the 
RSS of the uncertainties of both catalogs. 



POSITION UNCERTAINTY, arcseconds 


Fig. 5. Histogram of number of sources versus estimated position 
uncertainty: (a) right ascension; (b) declination 
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Submilliarcsecond VLBI Observations of the Close Pair 
GC 1342+662 and GC 1342+663 

D, D. Morabito 

T racking Systems and Applications Section 


Differential VLBI has been simultaneously performed on the source pair GC 1342+662 
and GC 1342+663 (4.4-arcminute separation) using S-band on the Goldstone/ 
Madrid baseline. These measurements were acquired on two separate observing sessions: 
30 December 1982 and 14 May 1983. The change in separation of GC 1342+662 relative 
to GC 1342+663 between the two epochs was 0.03 ±0.08 milliarcsecond. The differences 
of the relative position measurements between epochs of GC 1342+662 relative to 
GC 1342+663 were -0.29 ± 0.05 milliarcsecond in right ascension and. 0.14 ±0.09 milli- 
arcsecond in declination. These measurements demonstrate submilliarcsecond accuracy 
and repeatability. The discrepancies outside of the formal uncertainties could be attrib- 
uted to the intrinsic properties of the sources such as structure and to a lesser probability , 
proper motion. These discrepancies could also be attributed to excursions in UT1-UTC of 
about four times the quoted BIH uncertainty. 


I. Introduction 

VLBI reference frames of several milliarcsecond precision 
are being used by the Deep Space Network for spacecraft navi- 
gation (Refs. 1, 2, and 3). Quasars from these reference 
frames have served as reference beacons in differential VLBI 
experiments with the Voyager spacecraft, the Pioneer Venus 
Orbiter (Ref. 3), and the Viking Orbiter (XX Newhall, private 
communication). 

The potential limiting accuracy attainable for differential 
VLBI positional measurements, using these reference frames, 
can be reached by simultaneously observing sources with small 
angular separations. These close pairs of sources must have 
angular separations smaller than the beamwidth of each 
antenna, and therefore can be simultaneously tracked by both 
antennas of a VLBI baseline, greatly enhancing the cancella- 
tion of measurement errors. By observing several close pairs 


of radio sources over a several-year period, limits on reference 
frame stability can be determined. 

Observations and results demonstrating the submilliarc- 
second accuracy of this technique at 2.3 GHz using the DSN 
baseline, Goldstone/Madrid, on the close pair GC 1342+662 
and GC 1342+663 are presented. These two sources have 
been identified as quasars (Ref. 4), and have originally been 
identified as VLBI sources in the DSN full-sky survey program 
(Ref. 5). Total-flux densities measured during April 1983 
indicate that both of these sources have flat spectra (R. Perley, 
private communication). 

II. Observations 

The close pair GC 1342+662 and GC 1342+663 was ob- 
served at 2.29 GHz with the Goldstone/Madrid baseline. It 
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was observed by DSS 14 (64 m) and DSS 63 (64 m) on 30 
December 1982 for 6 hours and 40 minutes, and by DSS 13 
(26 m) and DSS 63 (64 m) on 14 May 1983 for 6 hours. The 
Goldstone/Madrid baseline has a length of 8400 km and a 
fringe spacing of 3 milliarcseconds. The observations were 
performed with right circular polarization. The receiver chain 
consisted of an S-band traveling wave maser followed by a 
special phase-stable S-band VLBI receiver, which converted the 
RF signal to an IF of 50 MHz. The NRAO Mark II VLBI 
recording system then recorded a 1 .8-MHz data bandwidth by 
digitally sampling at a 4-Mbits/s rate (Ref. 6). Digital sampling 
and phase stability of the receiver chain were controlled by 
hydrogen maser atomic clocks. The centroid position of the 
source pair was used to point the antenna. X-band observa- 
tions were not performed since the 4.4-arcminute angular 
separation of this source pair exceeded the half-power beam- 
width of the limiting antenna (2 arcminutes), thus making 
simultaneous observations useless. 

The data were correlated on the Caltech/JPL Mark II 
correlator. The postprocessing software yielded 60-second 
points of fringe phase for each source. The phases were then 
input to the JPL Orbit Determination Program (Refs. 7, 8, 
and 9), where an accurate phase model was applied. The 
phases were then corrected for cycle ambiguities and differ- 
enced. This difference phase was input to a least-squares algo- 
rithm where an integral number of cycles was determined and 
then added to each of the difference phases. The next iteration 
of the least-squares algorithm was a two-parameter fit that 
yielded the relative position offsets in right ascension and 
declination. 

III. Results 

Figure 1 shows the postfit difference-phase residuals for the 
30 December 1982 experiment, which have an rms scatter of 
0.028 cycles (90 juarcseconds). Figure 2 shows the postfit 
difference-phase residuals for the 14 May 1983 experiment, 
which have an rms scatter of 0.026 cycles (80 /rarcseconds). 
The scatters reported here are in agreement with predicted 
values and are dominated by the solar plasma. The scatter on 
the individual 60-second integrated phase points is about 
0.003 cycles and is consistent with expected noise levels due 
to system noise and ionosphere over these time scales. 

The correlated amplitudes were converted to source 
strengths using measured values of system temperatures and 
antenna efficiency parameters. The resulting correlated flux 
densities at 2.3 GHz are shown in Fig. 3 for GC 1342+663 
and in Fig. 4 for GC 1342+662. Peak correlated flux densities 
on 30 December 1982 were 0.35 Jy for GC 1342+662 and 
0.75 Jy for GC 1342+663, and on 14 May 1983 they were 
0.35 Jy for GC 1342+662 and 0.60 Jy for GC 1342+663. 


The correlated flux density behavior for GC 1342+663 is 
essentially flat (see Fig. 3) indicating a lack of detectable 
structure over the observed hour-angle range. There is, how- 
ever, a significant variation of correlated flux density over the 
interferometric hour angle for GC 1342+662 (see Fig. 4) 
indicating complex structure over the milliarcsecond resolu- 
tion beam of the interferometer. This variation is consistent 
with a model of two point sources with equal flux densities 
(0.175 Jy each) separated by 1.1 milliarcseconds and aligned 
along a ~29-degree position angle (defined north through 
east). The position angle is coincidentally aligned with the 
position angle between the two sources GC 1342+662 and 
GC 1342+663, suggesting the possibility of an intrinsic rela- 
tionship between the sources. 

The measured position in arclength of GC 1342+662 
relative to GC 1342+663 referred to epoch 1950.0 is given in 
Table 1 for both experiments. The arclength differences of 
these relative position measurements between experiments 
were -0.29 ±0.05 milliarcseconds in right ascension and 
0.14 ±0.09 milliarcseconds in declination. The results agree 
to within 5.8 and 1.5 times the rss of the individual standard 
deviations for right ascension and declination respectively. 
Marcaide and Shapiro (Ref. 10) had similar few-sigma discre- 
pancies between epochs of their relative position measure- 
ments of 1038+528A and 1038+528B. The change in separa- 
tion of the GC 1342+662 relative to GC 1342+663 between 
the two epochs was 0.03 ±0.08 milliarcsecond, implying that 
the above discrepancies are effectively a rotation in right 
ascension and declination space. Because most of the discrep- 
ancy is in right ascension, timing errors are immediately sus- 
pect. However, such errors would require deviations from BIH 
values of UT1-UTC of about four times the quoted uncer- 
tainty (Ref. 1 1). 

Source structure is the most likely cause of the measured 
discrepancy in the differential position measurements between 
epochs for GC 1342+662 and GC 1342+663. Proper motion 
is a possible cause for the discrepancy, but not very probable 
assuming current cosmological assumptions. To properly 
understand and distinguish between these effects, more obser- 
vations need to be performed over a several-year time period. 

Marcaide and Shapiro (Ref. 10), Shapiro et al. (Ref. 11), 
and Gorenstein et al. (Ref. 12) have performed similar experi- 
ments on other close pairs. Table 2 displays a summary of 
these experiments along with the one discussed in this paper. 

IV. Conclusion 

The improved accuracy and repeatability of simultaneous 
differential VLBI using the close pair GC 1342+662 and 
GC 1342+663 with the Goldstone/Madrid baseline has been 
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demonstrated on two separate observing sessions. The change 
in separation of GC 1342+662 relative to GC 1342+663 
between the two epochs was 0.03 ±0.08 milliarcsecond. 
The differences of the relative position measurements of 
GC 1342+662 relative to GC 1342+663 between the two 
epochs were -0.29 ±0.05 milliarcseconds in right ascension 


and 0.14 ±0.09 milliarcseconds in declination. Discrepancies 
between the relative measurements outside the formal uncer- 
tainties could be attributed to either source-structure effects 
or proper motion. Further monitoring of this close pair along 
with other close pairs will be useful for testing and placing 
limits on reference frame stability. 
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Table 1. Relative position measurements 


Source 

Epoch 

Relative 

right ascension, 3 
arcseconds 

Relative 

declination, 

arcseconds 

Total 

separation, 

arcseconds 

GC 1342+662 

30 Dec 1982 

-139.58033 

-225.17070 

264.83216 



±0.00004 

±0.00006 

±0.00006 

GC 1342+662 

14 May 1983 

-139.58062 

-225.17056 

264.83219 



±0.00003 

±0.00007 

±0.00006 


a In arclength referred to declination of GC 1342+662. 


Table 2. Comparison with other investigations 


Investigators 

Source pair 

Angular 

separation, 

arcminutes 

Measurement 

accuracy, 

Marc 

Ref. 11 

3C 345 
NRA0512 

30 

500 

Ref. 12 

0957+561 A 
0957+56 IB 

0.1 

50 

Ref. 10 

1038+528A 

1038+528B 

0.6 

10 

This paper 

GC 1342+662 
GC 1342+663 

4.4 

60 



0.10 



Fig. 1. Postfit residual difference phase, 30 December 1982, baseline DSS 14/DSS 63 



Fig. 2. Postfit residual difference phase, 14 May 1983, baseline DSS 13/DSS 63 




CORRELATED FLUX DENSITY, Jy CORRELATED FLUX DENSITY, 


(a) Dsc 30< 1982 



INTERFEROMETER HOUR ANGLE 


Fig. 3. Correlated flux density vs interferometer hour angle for GC 1342+663: 
(a) first experiment; (b) second experiment 



INTERFEROMETER HOUR ANGLE 


Fig. 4. Correlated flux density vs interferometer hour angle for GC 1342+662: 
(a) first experiment; (b) second experiment 
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Application of a sensitive hydrogen gas detector has led to the observation of a signifi- 
cant increase in the efficiency of atom production of a magnetically enhanced dissociator. 
An increase of about 7000% in the production of atoms in the hyperfine state desired for 
maser operation is observed at certain values of the applied magnetic field. A preliminary 
test of the dissociator in the maser has led to an observed increase of 6. 18 dB in the maser 
power. 


I. Introduction 

The development of the magnetically enhanced hydrogen 
gas dissociator was initiated in response to a need for improv- 
ing the operational life of the hydrogen maser. An improve- 
ment in the efficiency of the dissociator would result in re- 
duced gas consumption, extending the operational life of 
vacuum pumps. This would reduce the relatively large number 
of maser failures which result from the breakdown of Vacion 
pumps, as often encountered in the field. 

Previous investigations involving the application of an axial 
dc magnetic field on the rf excited dissociator had led to en- 
couraging results (Refs. 1 and 2). The results of these investi- 
gations were based on the determination of the flux of the 
atomic hydrogen effusing from the dissociator with and with- 
out the applied magnetic field. These measurements, however, 
could only yield relative values for the atomic flux; the abso- 
lute value of the hydrogen atom flux could not be obtained 
because of the nature of the measuring technique (Ref. 2). 

In the present work we report on our observation of the 
atomic hydrogen flux as measured directly with a hydrogen 
atom detector. The use of this detector permitted a further 


investigation of the dependency of the atomic flux on the 
strength of the dc magnetic field. The latter results, together 
with the details of the hydrogen atom detector are subjects of 
this article. 

Section II of this article presents the details of the hydro- 
gen atom detector and Section III provides the result of inves- 
tigations using the detector. The study will be summarized in 
Section IV followed by recommendations for future inves- 
tigations. 

II. The Hydrogen Atom Detector 

An efficient detection of atomic hydrogen may be accom- 
plished through measurement of the heat released when two 
hydrogen atoms combine to form a molecule. The amount of 
heat produced in the reaction H + H -*■ H 2 is large enough 
(4.38 eV/reaction) for detection via the measurement of 
change of resistance of a suitable wire. It can be shown that 
for a length 2 a of wire immersed in hydrogen the change in 
the resistance A R is given by: 

(RanyEaD) 

^ (yrkd) 
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Here a is the temperature coefficient of the wire, k is the 
thermal conductivity, n is the number of atoms per unit area 
per unit time, y is the probability of recombination for H on 
the particular material of the wire, 2 E is the heat of recombi- 
nation, d is the diameter of the wire, and D is the atomic beam 
diameter. 

In our studies the expected flux of atoms out of the disso- 
ciator is of the order of 10 15 atoms per second. The atom 
detector therefore had to be designed with a sensitivity high 
enough to measure such a flux accurately and efficiently. The 
apparatus, designed with the above constraint, consisted of 
two short segments of a thin (2.5 /am in diameter) platinum 
wire arranged in an ordinary dc bridge circuit (see Fig. 1). The 
platinum wire was prepared by etching the surface of a 
Wollaston wire with nitric acid to expose the thin platinum 
core. 

Only one of the segments was exposed to the hydrogen 
beam, so the difference in the voltage across this segment and 
that across the second segment could be amplified with a dif- 
ferential amplifier. In this way the detector was sensitive to 
a hydrogen flux as small as 10 11 . The detector and its elec- 
tronic elements were carefully shielded from the rf power used 
to excite the dissociator. This eliminated a relatively large 
noise problem which was initially encountered with the 
detector. 

Since the probability of the recombination of H on plati- 
num changes with surface contamination, provisions were 
made in the detector to periodically clean the surface of the 
wire through passage of a relatively large (2 mA) current. Thus 
the sensitivity of the detector was kept at a constant value 
throughout each run. 

The platinum wire detector was placed in a stainless steel 
vacuum chamber, evacuated with a Vacion pump. The de- 
tector was placed on a plate attached to a rotatable feed- 
through which could be used to move the plate and the 
detector in the horizontal direction. The rf excited source was 
placed external to the vacuum on a flange so that the holes in 
its exit collimator defined a vertical beam of atoms aligned to 
intersect the platinum wire. The source, rf exciter, and sole- 
noid providing the dc field are identical to those described 
previously (Ref. 2). 


III. Results and Measurements 

As a first test of the hydrogen atom detector we attempted 
to replicate the previously reported results obtained with the 
magnetically enhanced dissociator. Under the normal condi- 


tions of source pressure and rf power the application of the dc 
magnetic field resulted in an increase in the source efficiency 
which agreed closely to the values reported in Refs. 1 and 2. 

The increased efficiency of the atom production with the 
magnetically enhanced dissociator was large enough (up to 
70% for certain values of pressure and source power) that an 
optimization of the beam optics for application in the maser 
seemed appropriate. We therefore modified the dissociator test 
apparatus to support a state selecting hexapole magnet at the 
exit end of the source. This magnet selects atoms of hydrogen 
that are in the hyperfine ground state having F = 1 and m F = 1 
and 0 and focuses them on the entrance hole of the storage 
bulb in the maser. The remainder of the atoms in the F = 0 
state are defocused and therefore do not enter the storage 
bulb. This selection is necessary because all atoms that are not 
in the F = 1 , rrif = 1 state do not participate in the maser oscil- 
lation, and yet they contribute to line broadening through spin 
exchange collisions. 

The observed atom production efficiency of the magne- 
tically enhanced dissociator with the state selecting hexapole 
magnet in line was quite unexpected. There appeared to be a 
very significant increase in the production of atoms focused by 
the state selector. While the flux of atoms which are focused 
by the hexapole magnet is typically of the order of 10 15 
atoms per second, the application of the dc magnetic field on 
the source increased the flux of focused atoms by as much as 
7000% for a particular value of the applied field. Beyond this, 
a plot of atomic flux vs. magnetic field revealed pronounced 
structures which in certain cases implied a decrease to zero 
for the flux of focused atoms. This result is presented in 
Fig. 1. 

As indicated in Fig. 2, there are a number of features in the 
plot of flux vs. field. There are two broader peaks at about 
150 and 500 Gauss, and a rather sharp peak at about 230 
Gauss. A number of valleys are also present including one at 
about 30 Gauss, which indicates a decrease in the flux of the 
atoms to zero value. These results were checked for reproduci- 
bility many times and were verified each time. They did 
however indicate changes in the amplitude of the peaks as 
parameters of rf power and hydrogen pressure were changed. 

In view of the remarkable observation presented in Fig. 2, 
we tested the measuring apparatus to ensure the reliability 
of results and absence of artifacts. We also examined possible 
changes in the focusing properties of the hexapole magnet as 
a result of the presence of the applied dc field. In all cases we 
reached the conclusion that the results obtained were not due 
to any artifacts associated with the technique of measurement 
or the measuring apparatus. 
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The observed results then clearly indicate that for certain 
values of the dc magnetic field the total efficiency of atom 
production of the source is increased. But more importantly, 
the results indicate a much larger increase, at certain field 
values, in the production of those atoms in the desired hyper- 
fine state. 

This last contention was tested by retrofitting a maser with 
the magnetically enhanced dissociator and measuring the 
maser output power. Because of the physical requirement that 
the solenoid producing the dc field be outside the vacuum 
chamber, the dissociator had to be placed external to the vac- 
uum chamber of the maser. This produced an increase in the 
path length of hydrogen atoms by a factor of nearly two, as 
compared to the path length of atoms when the dissociator is 
in its usual configuration inside the vacuum chamber. 

The result of this test indicated that the maser oscillation 
ceased when the dc field was at a value which indicated zero 
atom production in the test bed (about 30 Gauss). As the field 
increased the maser power increased and reached a level 
6.18 dB above its value without the magnetic field for fields 
larger than about 120 Gauss. The 6.18 dB increase in the 
power level, however, remained constant for all field values 
above 120 Gauss. 

While this observation verified that the efficiency of atom 
production in the desired hyperfine state is significantly in- 
creased with the application of a dc field, it failed to establish 
a clear correlation between the observed peaks of Fig. 1 and 
the maser power. We believe that a clear understanding of this 
observation will be possible only when modifications are made 
to ensure that the path length of the atoms in shortened to 
be equivalent to the length in a normal maser configuration. 

At this point it may be appropriate to discuss the possible 
causes of the behavior of H-atom flux with applied dc field 
when the hexapole magnet is placed in the beam path. A pos- 
sible explanation of the observed result is the preferential 
pumping of the levels of the ground and excited electronic 
states of hydrogen by the rf power and the cyclotron radia- 
tion of electrons circulating in the applied magnetic field. For 
example, the peak at 230 Gauss could be attributed to pump- 
ing of the 25, F = 1 , m f = - 1 state, which is forbidden to make 
a transition to the desired ground state, to the 2P, F = 1 , 
rrif = 1 state, which is allowed to make the transition. While we 
have similar explanations for other features in the spectrum 
depicted in Fig. 2, we believe that a totally conclusive explana- 


tion should await further studies with the source in the modi- 
fied test chamber. 


IV. Summary 

In this article we have described some details of a hydrogen 
atom detector which has enabled a comprehensive investiga- 
tion of the atom production of the magnetically enhanced 
dissociator. Our observations have revealed that while a sig- 
nificant increase in the production of atoms is obtained, the 
magnetically enhanced dissociator can increase at a greater 
proportion the number of atoms in the desired hyperfine state. 
This observation was further supported when a large (415%) 
increase was produced in the power output of the maser, 
which was operated in a test configuration with the magneti- 
cally enhanced dissociator. Finally, a preliminary explanation 
was put forth to describe the physical processes which can lead 
to the features observed in the spectrum of atomic flux vs. 
the magnetic field. 

While the above results are to be considered quite preli- 
minary, they do imply rather significant contribution to 
H-maser technology. As discussed in previous works (Refs. 1 
and 2), the use of the axially applied magnetic field with the 
dissociator can enable maser operation at a significantly re- 
duced hydrogen throughput. This will increase the operational 
life of the vacuum pumps without degrading the maser power 
output. The observation of the significant increase of the 
atoms in the desired hyperfine state further implies that for 
a fixed output power of maser (and therefore a fixed number 
of desired atoms in the storage bulb)' the magnetically en- 
hanced dissociator can lead to significantly reduced line broad- 
ening due to spin exchange collisions. This may be brought 
about since the applied field can evidently increase the ratio 
of “desired atoms” to “unwanted atoms” in the storage bulb 
of the maser cavity. Therefore significantly higher line Q’s 
should be expected with the use of the magnetically enhanced 
dissociator, when all other parameters are kept fixed. 

The encouraging results of this investigation have led us to 
believe that the above investigations should continue to fully 
determine the ultimate influence of the enhanced source on 
the operation of the H-maser. Currently further studies are 
underway to shorten the optical path of the atoms to the 
normal length as well as to pursue investigations with the test 
chamber. 
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This article reports the status of the X-band uplink development program from a 
systems ' viewpoint. 

The hardware, designed and built under the program, is now in place at DSS-13. Sys- 
tem stability testing is underway. 

An X-band receive capability is being implemented on both the Galileo and the Venus 
Radar Mapper spacecraft. Experiments will be conducted that will serve both to demon- 
strate the capabilities of the X-band uplink system and to enhance science objectives on 
those missions. These objectives are briefly described. 


I. Introduction 

The DSN undertook a phased development program in 
1979, under RTOP 64, to add an X-band uplink capability to 
the deep space network. The first phase was to design and 
build highly stable ground station subsystems for installation 
and checkout at DSS-13. A second phase was to arrange a 
flight demonstration of a complete X-band uplink system to 
verify predicted performance in the space environment. 

In order to plan and coordinate the activities, a systems 
work unit was established. Under it, various tasks were accom- 
plished, utilizing a design team representing the various hard- 
ware subsystems, space environment disciplines, and the user 
community. The main efforts were directed towards the fol- 
lowing tasks: 

(1) Quantifying the requirements imposed on the overall 
system by potential users of the X-band uplink, par- 
ticularly in the navigation and radio science areas. 


(2) Establishing error budgets for the telecom functions of 
telemetry, ranging and Doppler, including error esti- 
mates, with tolerances, for ground and spacecraft 
hardware, and the natural media effects due to the 
troposphere and interplanetary and ionospheric 
charged particles. 

(3) Maintaining extensive liaison with NASA flight proj- 
ects to develop mutually beneficial X-band uplink 
experiments. 

(4) Preparing an Experiment Requirements Document, 
defining a flight demonstration of the X-band uplink, 
including preparatory ground testing and inflight cali- 
brations as well as experiment sequences. 

Much of the early systems work was summarized in TDA 
Progress Report 42-62 (Ref. 1). 


24 



Subsequently, the NASA spacecraft of the ISPM mission, 
which was initially to have carried the X-band uplink demon- 
stration experiment, was cancelled. This, report tells of the 
developments which have taken place since then, including 
the status of DSS-13 X-band uplink activities and plans of the 
DSN now to support X-band uplink operations on the Galileo 
and Venus Radar Mapper missions. 

II. Status of X-Band Uplink Activities at 
DSS-13 

A. Hardware Status 

The ground station hardware (Ref. 2) has been built and 
tested at the subsystem level, and is now installed at DSS-13. 
Included are the following subsystems: 

(1) The X-band uplink exciter, X/X and X/S Doppler 
extractors and mixers and test signal translators which 
coherently translate the X-upffnk frequency to X- and 
S-downlink frequencies for in-station tests. 

(2) The 20-kW klystron power amplifier, with buffer 
amplifier driver, complete with power supplies, con- 
troller and protective circuitry. 

(3) A new microwave antenna feed design to accommo- 
date S- and X-band uplinks and downlinks. 

(4) An X-band diplexer, specially designed to minimize 
noise temperature contributions. 

(5) Block III receiver modifications. In order to realize the 
potential advantages of the X-band uplink in reducing 
phase jitter due to charged particle variations in the 
solar wind and the ionosphere, the station will use a 
Hydrogen-maser timing reference and will incorporate 
a temperature-controlled Multi-mission Receiver 
(MMR) front end. Concepts consistent with unat- 
tended operations have also been employed. 

The components that have been redesigned in the X-band 
uplink R & D effort, and are now in place at DSS-13, are indi- 
cated in the block diagrams. Fig. 1 shows the station as it 
would be configured to support an X-band demonstration with 
the Galileo mission, while Figs. 2 and 3 show the hardware in 
more functional detail for the X-band uplink/X-band downlink 
and X-band uplink/S-band downlink paths respectively. 

B. DSS-13 System Test Status 

Individual subsystems were checked out prior to instal- 
lation at DSS-13 to verify that they met their design 
specifications. 


In addition, a stepwise test plan was formulated to evaluate 
the total X-band uplink system. The approach emphasizes 
stability tests as the most informative. The testing is currently 
underway and will continue through 1984. The test sequence 
is as follows. 


1. Baseline exciter subsystem stability. Exciter subsystem 
stability tests will serve as the baseline. Special tests were con- 
ducted prior to installation at DSS-13, in the Interim Fre- 
quency Standards Test Facility, using laboratory Hydrogen- 
masers as frequency references. Phase stability was measured 
at the X-band uplink frequency and at the downlink X- and 
S-band frequencies coherently derived from the uplink using 
the test signal translators that are part of the exciter subsys- 
tem. The phase data was processed to yield a square root Allan 
variance record versus sampling time. 

The results of these measurements are encouraging. Fig. 4 
shows a typical plot of the X-band uplink frequency stability 
over sampling times from 1 to over 20,000 seconds. The mea- 
surements have established that these exciter subsystem ele- 
ments operate within the stability bounds specified for their 
design. Typical (A///) values below 3 parts in 10“ 15 were 
obtained for integration times over 1000 seconds. 


2. X-band uplink stability at DSS-13. The exciter stability 
measurements are being repeated at the station at this writing. 
Because of the hardware’s location in the antenna feedcone, a 
portable, highly stable microprocessor data system had to be 
designed for these in situ measurements. This has been built 
and tested, verifying an acceptably close match in stability to 
the laboratory instrumentation described previously. 

The uplink frequency stability measurements will next be 
extended to various points moving progressively from the 
exciter to the buffer amplifier and the 20-kW klystron power 
amplifier. These measurements will identify additional noise 
increments added along the path to the antenna. 

The measurements will be performed first by open loop, 
and then by closed feedback path back to the exciter from, for 
example, the output of the klystron power amplifier. The 
feedback loop has been included to compensate for relatively 
slow phase changes in subsystems that are especially subject to 
variations in ambient temperature. The loop should essentially 
clamp the phase’s excursions to those of the exciter, and its 
success depends on the stability in the feedback path and the 
gain in the loop. Note that the loop time constant should be 
long enough to permit passage of command and ranging 
modulations. 
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3. Receiver testing. The next tests will bring in portions of 
the receiver. For these tests, downlink frequencies will be gen- 
erated using the X-up/X-down and X-up/S-down frequency 
translators operating on a sample of the X-band uplink ex- 
tracted via a directional coupler. The receiver frequencies 
would be reinserted in the downlink path, also via couplers, 
after being adjusted to the proper level. Stability measure- 
ments would be made at these RF frequencies down to the 
point where downconversion occurs. At this point, viz the 
outputs of the X- and S-band masers, the DSS-13 station con- 
figuration is representative of the configuration which will be 
implemented in the High Efficiency (HEF) operational sta- 
tions, DSS-15, DSS-45 and DSS-65. Therefore, stability mea- 
surements made at DSS-13 at the maser output should also be 
representative of those expected at the operational stations. 
From this point, however, the configurations differ. The HEF 
stations, shown in Fig. 5, will route closed-loop data for radi- 
ometric measurements, such as Doppler, through conventional 
34-m station hardware consisting of an X- to S-downconverter 
followed by a Standard Block III S-band receiver. 

In the HEF stations a separate signal path will utilize the 
temperature-controlled and phase-calibrated Multimissioa 
Receiver (MMR) RF frontend. This path leads to an open-loop 
receiver within an Occultation Data Assembly (ODA) for use 
by radio scientists to collect data for further processing in 
experiments such as the extraction of Doppler signatures for 
gravitational wave searches to be conducted on the Galileo 
mission. 

In contrast, DSS-13 will utilize only a closed-loop receiver 
at this time, but will enhance the Block III receiver perfor- 
mance by adding modifications. The receiver modifications 
consist of a MMR frontend and an upconverter/downconverter 
section prior to interfacing with a Standard Block III receiver. 
Figs. 1 and 5 compare the configurations at DSS-13 and the 
HEF stations respectively. 

4. Doppler stability tests at DSS-13. The next step is to 
bring the entire closed-loop DSS-13 system into play in an 
end-to-end stability test. The most meaningful point in the sys- 
tem to make this test is felt to be the output of the Doppler 
mixer, at either the 5 MHz or 1 MHz bias frequency point. 
This measurement would involve all of the transmit and 
receive hardware at the station, except the antenna, and would 
be indicative of the system noise “floor” for the station 
components. 

5. Inflight demonstration. The final test of system stability 
will be accomplished using an actual space mission. This will 
first occur during the Galileo flight and will bring into play the 
elements of the total system not available before, including the 


station antenna, the total spacecraft receive/transmit system, 
and the intervening media. 

Emphasis in the initial inflight demonstration test would 
again be on two-way Doppler stability using DSS-13 and 
Galileo. In order to obtain the best achievable conditions, the 
experiment would be performed at Galileo’s first opposition. 
This view period has two advantages from a stability stand- 
point. It would provide a “quiet” environment (nighttime 
viewing) and would occur while the spacecraft is still close to 
earth; hence, it would feature strong signal conditions. The 
test should, in effect, establish a new system noise “floor,” but 
this time for the total end-to-end system, complete with space- 
craft and media. 

Subsequent X-band uplink technology demonstration 
sequences would test other telecommunications functions, 
such as commanding and two-way coherent telemetry. 

C. Supporting Analysis: Doppler Error 
Source Modeling 

As noted, end-to-end tests are planned during 1984 to 
measure the phase stability of the DSS-13 X-band uplink sta- 
tion configuration. Preceding this, an analysis task was under- 
taken to guide the interpretation of the test results. In it, 
equations were derived for the phase transfer functions relat- 
ing the contributions of error sources in the system to various 
output points including the X-band and S-band biased Doppler 
mixer outputs. A follow-up task is recommended for the 
future. This will involve obtaining and assigning phase power 
spectra for the major phase error sources such as VCOs, fre- 
quency synthesizers, cables and the H-Maser reference so that 
accurate quantitative predictions can be made. 

The first part of this analysis, the derivation of the transfer- 
function equations, has been completed and was reported in 
TDA Progress Report 42-76 (Ref. 3). It considered two modes. 
One, for station testing, uses a test signal translator (a zero- 
delay device) as an uplink/downlink turnaround device. The 
second configuration models a system with an actual space- 
craft at a significant round trip light time from the DSN sta- 
tion. The results indicate that system tests at the station can 
be useful in localizing many of the phase errors contributed 
by the station hardware, but cannot accurately reflect the net 
Doppler error that will occur during actual spacecraft tracking. 
Aside from space environment factors, lack of fidelity in simu- 
lating real mission conditions is due to two main differences 
between test and flight conditions: 

(1) The round-trip times between an uplink transmission 
and subsequent downlink reception in a mission cannot 
be simulated in station tests. Therefore, in the station 
tests, phase error sources that are common to both the 
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uplink signal and the Doppler extractor will be highly 
correlated because they are generated from the same 
reference. These errors will be largely cancelled out in a 
process such as the conventional two-way Doppler 
measurement, wherein a Doppler reference and an esti- 
mate of a received signal are subtracted. However, 
some errors that enter only one part of the system, 
such as the phase error contributions of the X-band 
and S-band receiver VCOs and frequency synthesizers 
(DCOs) will be accurately identified in station tests. 

(2) Translator turnaround devices do not simulate the 
spacecraft transponder function accurately. A tran- 
sponder produces a downlink frequency by multiplying 
the uplink frequency by the required turnaround ratio 
(880/749 for X/X and 240/749 for X/S). On the other 
hand, the test signal translators, used in the DSS-13 
implementation, pick off samples of the uplink fre- 
quency from selectable points in the exciter power 
amplifier chain, and, using mixers, simply add or 
subtract coherent components required to give the 
downlink X- or S-frequencies. These components of 
the downlink frequencies are derived directly from the 
H-maser and DCO references and are not subject to 
errors generated in the power amplifier section of the 
uplink. 

Thus, an X- or S-downlink generated by a test signal 
translator would have essentially the phase noise of the 
uplink. This is not the same as found in a downlink 
generated by a transponder where the phase noise 
would be multiplicatively increased by the X-down/ 
X-up turnaround ratio (880/749) or decreased by the 
S-down/X-up turnaround ratio (240/749). This error, 
of course, could be eliminated by using a transponder 
for station tests, when available. 

The above conclusions are somewhat qualitative, but do 
point out what to expect. Within the restrictions noted, the 
phase transfer functions derived in the analysis can be used to 
calculate quantitative results, given models for the phase 
power spectrum of each of the phase error sources. 


III. Flight Projects 

Both the Galileo and Venus Radar Mapper (VRM) projects 
have made decisions to fly an X-band uplink receive capability. 
Support had been provided to them prior to, and since, that 
time by the X-band uplink system work unit. This consisted 
primarily of technical inputs, such as estimates of radiometric 
accuracies using X-band vs S-band uplinks, and preliminary 
estimates of spacecraft hardware impacts and costs. 


The X-band uplink sequences during these two missions will 
serve two purposes. On the one hand they will compare the 
performance of telecom functions using an X-band uplink 
system to the predictions of the systems design team. And 
they are expected to enhance selected flight project science 
objectives. Because science enhancements, rather than oper- 
ational objectives, are primary for the X-band uplink on these 
two missions, present plans call for single string (nonredun- 
dant) implementations of X-band receive hardware on both 
spacecraft. 

Fortuitously, the science experiments on Galileo and VRM 
present the opportunity to test an X-band uplink system in 
two different natural environments, using different station 
configurations. The gravitational wave search by Galileo will 
be conducted during solar oppositions (nighttime viewing) at 
long range. Doppler signatures with slow variations are 
expected. Open-loop receivers and extensive processing will be 
used. 

The VRM gravity experiment, on the other hand, will 
necessarily be conducted during daytime viewing, with the 
sun-earth-probe angle often quite small. The Doppler signature 
in this experiment will be acquired using the closed-loop 
receiver and the sampling time will be short, i.e., typically 
5 seconds. 

Some of the objectives achievable through the use of an 
X-band uplink are summarized: 

A. Galileo: X-Band Uplink Demonstration and 
Gravitational Wave Search 

The objectives to be accomplished with the use of X-band 
uplink on the Galileo Orbiter can be included in three general 
areas. They are: 

(1) A demonstration of improved performance in the 
various telecommunication links relative to S-band 
uplink. 

(2) A demonstration of the DSN capability to maintain 
flight operations with operational, high-power X-band 
transmitters. 

(3) A significant enhancement in the capability to perform 
the Gravitational Wave Experiment. 

The plan for the demonstration of the advantages of 
X-band over S-band in telecommunication link performance is 
essentially the same as that which was previously intended for 
the ISPM U.S. spacecraft. Since these objectives were defined 
in detail in a previous report, (Timpe, C. T., private communi- 
cation) they will be summarized here. The telecommunications 
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links on which improved performance will be demonstrated 
with X-band uplink are: 

(1) Two-way coherent carrier tracking (Doppler stability). 

(2) Command threshold. 

(3) Two-way ranging. 

(4) Two-way coherent telemetry data-rate. 

The expected improvements in performance will result 
from the following two main factors: 

(1) The net effect of a higher frequency on a one-way link 
(having the same size antennas, the same transmit 
power, and the same receive sensitivity) is a perfor- 
mance gain by a factor of approximately the square of 
the frequency ratio. 

(2) The susceptibility of radiated energy to phase fluctua- 
tions induced by charged particles in the ionosphere, 
and in the interplanetary plasma decreases with higher 
frequencies. 

The enhancement in the capability to perform the Gravita- 
tional Wave Experiment will result from the improvement in 
Doppler fractional frequency stability with X-band uplink. 
The system goal for overall fractional frequency stability on 
the ISPM ESA spacecraft with S-band uplink is 3 X 10“ 14 (the 
Allan variance on A/// for a 1000 second time interval). The 
corresponding goal for Galileo with X-band uplink is 5 X 
10~ 15 , or nearly an order of magnitude improvement. The 
two spacecraft will perform a joint gravitational wave (GW) 
search during their coincident solar opposition with the Earth 
in October of 1987. It is anticipated that GW activity, which 


may not be discernible in ISPM Doppler data by itself, may be 
very significant in verifying the existence of GW activity 
present in the Galileo Doppler data which can be time- 
correlated. 

B. Venus Radar Mapper: X-Band Doppler Data for the 
VRM Gravity Experiment 

In the past, S-band Doppler data have been the primary 
source for gravity field information. The experimenter relied 
strictly on what the navigation transponder produced. How- 
ever, the major error in the S-band data is the effect from 
interplanetary media (plasma). The plasma produces ampli- 
tudes and frequencies that can be interpreted as gravity 
anomalies, when features near the resolution limits are being 
investigated. 

With an X-band Doppler system these errors and high fre- 
quency noise can be reduced by a factor of 5 to 8 and much 
better resolution and statistical results on all features can be 
obtained. Presently, a gravity map of Venus has been extracted 
from the Pioneer Venus Orbiter data. A more detailed map 
will be obtained from Venus Radar Mapper data due to the 
X-band Doppler system and the smaller eccentricity of its 
orbit. 

Gravity data are generally correlated with topographic 
results and then inferences are made on internal structure. For 
example, thickness of a crust or depth of isostatic adjustment 
are calculated at many areas. These in turn provide insight as 
to a general global structure. Also many topographic features 
can be classified by their gravity signatures, and, by using 
other geological evidence, inferences can be made on construc- 
tion and age. 
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Fig. 1. X-band uplink functional configuration for DSS-13 
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Fig. 2. Detailed functional block diagram for the X-band uplink/X-band downlink hardware at DSS-13 

(Example frequencies are for channel 14) 
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Effects of NRZ-M Modulation on Convolutional 

Codes Performance 

L. Deutsch, F. Pollara, and L. Swanson 
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Non-Retum-to-Zero Mark (NRZ-M) modulation is often used to resolve data sense in 
suppressed carrier telemetry systems. This is because such systems are subject to half 
cycle slips that result in complementing the encoded data stream. The performance of 
coded telemetry systems with NRZ-M is sensitive to the order in which the various opera- 
tions are done. This means that a system that demodulates the NRZ-M waveform and 
then decodes will perform differently from a system that does the decoding first. In this 
report, the performance of the NASA standard (7, 1/2) convolutional codes is determined 
for several systems using NRZ-M. In addition, several different demodulation schemes for 
NRZ-M are considered. It is shown that, even for the best soft-decision method examined, 
there is a 2.7 dB loss at a decoded bit error rate of 5 X 10~ 3 if the NRZ-M demodulation 
occurs before rather than after Viterbi decoding. 


I. Introduction 

Emerging coding standards for interagency projects, such as 
those adopted by the Consultative Committee for Space Data 
Systems, will include Non-Return-to-Zero Mark (NRZ-M) 
modulation (Ref. 1) on suppressed carrier links. This type of 
modulation, used to avoid serious errors in case of half-cycle 
carrier synchronization slips, is defined as follows. 

Suppose that { x n } is a binary stream of bits from (0,1). 
This conventional binary representation of data is sometimes 
called “Non-Return-To-Zero Level” or NRZ-L modulation. 
The NRZ-M encoder output O',,} is defined by 

y - x ©.y . 

where © denotes addition modulo two. 

Decoding (i.e., translation back to the original NRZ-L 
binary stream) is done by noting that x n - y n ®y„- 1 - 


The present DSN uses residual carrier transmission exclu- 
sively, which is subject only to full-cycle slips. In the DSN case 
NRZ-L modulation is therefore appropriate. 

If a suppressed carrier system were used with NRZ-L modu- 
lation, a half-cycle slip, which corresponds to an inversion of 
the data stream, could easily destroy large blocks of data. 
With NRZ-M the inversion would cause only isolated errors. 
If there is an outer block code, such as the NASA standard 
Reed-Solonron code, then these isolated errors are likely to be 
corrected by the block decoder. 

This article studies the effects of NRZ-M modulation on a 
(7, 1/2) convolutionally coded system. (The coding and decod- 
ing for such systems is not discussed here but a good presenta- 
tion may be found in Ref. 2.) In particular, two different 
schemes for implementing NRZ-M will be examined. The first, 
called Inner NRZ-M, introduces NRZ-M encoder and decoder 
just before and after the channel, so that the NRZ-M modula- 
tion operates on convolutionally encoded bits. The second 
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scheme introduces encoder and decoder before the convolu- 
tional encoder and after the Viterbi decoder, respectively, and 
will be called Outer NRZ-M. (Outer NRZ-M is an option only 
in the case of transparent codes [Ref. 2] such as the NASA 
standard convolutional codes.) 

A detailed description of these two schemes is in Sec- 
tions III and IV. Their performance with various decoding 
methods was determined by analysis and software simulation. 
These are described in Sections V and VI. The conclusions of 
the study follow immediately in Section II. 

II. Conclusions 

The addition of NRZ-M outside the convolutional channel 
(Fig. 1) causes only a slight drop in performance on a Gaussian 
channel whether symbols are hard- or soft-quantized, and so is 
a reasonable solution to the problem of half cycle slips in the 
case of suppressed carrier tracking. 

The other possible implementation of NRZ-M, inside the 
convolutional channel (Fig. 2) causes far more performance 
degradation, and should be avoided. 

Numerical results are shown and explained in Section VI. 

III. NRZ-M Modulation 

This type of modulation, widely used to avoid phase ambi- 
guities of half cycles or multiples, is also and more properly 
referred to as “differential encoding” (Ref. 1). 

Suppose that { x n } is the binary data stream to be encoded. 
Then the encoder output {y n } is defined as in Section I by: 

y = x (±)y , 

y n n J n-1 

If we consider the sequence {x n } to be bipolar data (i.e., each 
x n is taken from the set {- 1 ,+ l } rather than (0,1 }), then we 
may rewrite this expression as 

y - -x y , 

This notation is used for the remainder of this article. 

Let the received sequence be { y n }, where y n is the result 
of passing y n through a binary symmetric channel. In order to 
restore the original data, we must appropriately decode the 
received sequence {y n }. Let the NRZ-M decoder output be 
{£„}, where 

x = -y y , 

n 'nr n— 1 


Then, if there are no errors in the received sequence, i.e., 
y n = y n for each n, we have 

A. A A 

x = - y y . 

n ' n' n—l 

= x y v t 
n' n-r n - 1 

= X 

n 

as required. Note that if both_y fJ _ 1 and y n are complemented, 
we still have x n = x n . This invariance under signal inversion is 
the property that makes NRZ-M an attractive algorithm for 
suppressed carrier systems. 

The encoder/decoder algorithm just described can be 
directly used to implement the Outer NRZ-M scheme. 

The same is true for the Inner NRZ-M, if the received data 
is hard quantized. Otherwise, for a soft quantized receiver, we 
need to modify the NRZ-M decoder in some way. 

With NRZ-L data, alternate symbols from the convolutional 
encoder are usually inverted. This increases received symbol 
transition probability, resulting in improved symbol synchroni- 
zation (Ref, 3). In the case of Outer NRZ-M encoding, such a 
system would increase symbol transitions in the same way, 
with no change in the results of this report. In the case of 
Inner NRZ-M encoding, alternate symbol inversion schemes do 
not guard against long strings without transitions. 

IV. Soft-Decision NRZ-M Decoding 

The implementation of Inner NRZ-M encoding/decoding 
requires either a new Viterbi decoder that operates on NRZ-M 
encoded symbols, or an algorithm for the NRZ-M decoder to 
send soft quantized symbols to the Viterbi decoder. Conven- 
tional NRZ-M decoding described in Section III produces hard- 
quantized symbols as an output, but use of these in the Viterbi 
decoder would cause a substantial loss in Viterbi decoder per- 
formance. An optimal system based on channel probabilities 
is computationally cumbersome. 

One system which has been suggested is 

x n = sgniy^yyj 

= -(^ViVlVjl W 

This algorithm is the one used in the Symbol Synchronizer 
Assembly (SSA) in the DSN, according to the Symbol Syn- 
chronizer Assembly Technical Manual, Operations and Mainte- 
nance, and it agrees in sign with hard-quantized NRZ-M 
decoding. 
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Another system which agrees in sign with hard-quantized 
NRZ-M decoding is 


-yjn-l 


( 2 ) 


This second system, which appears more natural to us, per- 
formed better in the Viterbi decoder simulations (see Sec- 
tion VI). 


V. A Theoretical Comparison of the Inner 
and Outer Schemes 

In the case that received information over the channel is 
assumed to be hard quantized, a mathematical analysis of 
inner and outer NRZ-M systems at high SNRs is possible. This 
analysis consists of considering channel error sequences that 
produce Viterbi decoder bit errors. The probability of decoded 
bit error is then approximated by a sum over these sequences 
of their probability of occurrence multiplied by the number of 
bit errors caused by each one of them. 

Consider the Inner NRZ-M scheme. Also, assume that the 
sequence that is sent is the all zero sequence. Since both 
Viterbi decoding and NRZ-M decoding are linear operations, 
there is no loss of generality in this assumption as long as the 
system is synchronized. Since we are only considering hard 
quantization in this analysis, the channel may be modeled as a 
binary symmetric channel (BSC) with transition probability p. 

In the NRZ-L system, the Viterbi decoder (either JPL 
[Ref. 4] or Goddard conventions) is capable of correcting any 
error pattern of Hamming weight less than five (Ref. 4). A 
computer search was performed to determine sequences that 
cause bit errors when they are passed first through an NRZ-M 
decoder and then the Viterbi decoder. Since a weight two 
channel error pattern can only become at most a weight four 
pattern after NRZ-M decoding, all such patterns are corrected 
by the inner NRZ-M system. All weight three channel error 
patterns of length up to 16 symbols were checked by using a 
software Viterbi decoder. It was found that there are weight 
three sequences that produce bit errors in both the JPL and 
Goddard decoders if they were first NRZ-M decoded. In fact, 
there were sequences that produce eight and five bit errors, 
respectively, for these two codes. A complete listing of these 
sequences is shown in Table 1 . 


e w . ' E i>(",)p 3 (i-p )‘ 3 

i= 1 

The channel transition probability p is given by 

P = Q{\/2EJN 0 ) 

= Q(VeJn o ) 

where E s is the energy in a channel symbol and Q is the 
Gaussian error function 


r oo 

Q(x) = — — r / exp {-v 2 jl)dv 

J x 

For large x we may use the approximation (Ref. 5) 

Q(x) = — exp(-x 2 /2) 

\[2ti x 

Let a = E b /N a . Then, for large a we obtain 


'b., 


E b(v i } 

j= i 


(l/27ra) 3/2 exp(-3a/2) 


As a becomes large, this probability behaves like its expo- 
nential term. Hence, asymptotically, 


P bit = exp(-3a/2) (3) 

for both the JPL and Goddard codes. In comparison, these 
codes, without NRZ-M on the channel, behave asymptotically 
like 


^bit = ex P(~ 5a / 2 ) ( 4 ) 

Suppose that the bit SNR required to achieve a certain error 
rate for the conventional system is |3 and that the SNR re- 
quired to get the same performance in the Inner NRZ-M sys- 
tem is a. Then by Eqs. (3) and (4), 

exp(-3a/2) = exp(- 5/3/2) 


or 


a = ( 5 / 3 )/3 


At very high SNRs, the probability of a weight three 
channel error sequence occurring at any time is approximately 
p 3 . Suppose the weight three error sequences are {v ( } 
0=1,2, ...,«) and that the number of decoded bit errors that 
each produces is h(v ; ). Then for very high SNRs, the probabil- 
ity of decoded bit error is given by 


This means that the inner NRZ-M scheme should behave about 
2.2 dB worse than the NRZ-L system in the limit as SNR be- 
comes infinite. For “reasonable” signal-to-noise ratios, the dif- 
ference is even greater because of the terms other than the 
exponential terms. In particular, to obtain bit error probability 
0.005, a = 2.04)3, about a 3 dB loss. 
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If NRZ-M decoding is added after the Viterbi decoder, then 
the bit errors produced by the decoder are at most doubled. 
This constant factor is not significant in asymptotic behavior. 
Therefore the Inner NRZ-M scheme should also behave 2.2 dB 
worse than Outer NRZ-M for very high SNRs. 


VI. Simulation Results 

In order to determine the performance of various telemetry 
system configurations involving NRZ-M and convolutional 
coding, software simulations were carried out. The software 
modules included a convolutional encoder, Viterbi decoder, 
NRZ-M encoder, NRZ-M decoders for each of the NRZ-M 
decoding algorithms described above, and a Gaussian channel 
simulator. These were written in C-Language and run on a 
VAX 1 1/750 computer under the UNIX operating system. 

The configuration used for the simulation is shown in 
Fig. 3. The “generator” output is a stream of 0’s and l’s. The 
“coder” is a convolutional encoder implementing one of the 
standard NASA (7,1/2) codes. The “display” routine compares 
the delayed information bits with the decoded bits and gener- 


ates the statistics. In addition there are NRZ-M coders and de- 
coders which may be inserted in various places. 

Figure 4 shows all configurations used in the simulation. 
Results are shown in Figs. 5 and 6, with decoded bit error rate 
as a function of the bit signal-to-noise ratio, E b /N 0 , for each 
scheme. The simulations assume perfect carrier and subcarrier 
tracking and perfect Viterbi decoder node synchronization. 
The Viterbi decoder path memory length was 32 bits. 

The probability of error P E is measured on a random se- 
quence of N = 1,000,000 bits or until the 95% confidence 
interval C, defined as 

C = 1.96\/P e ( 1-P e )IN 

is less than 5% of P E . Therefore the true P E is inside the range 
P E ±C with 95% confidence. Table 2 shows a comparison of 
the various schemes for bit error probabilities of 1 0 — 2 , 
5 X 10 -3 , and 10~ 3 . 

As expected, the Outer NRZ-M scheme outperforms the 
Inner NRZ-M scheme, for both the JPL and Goddard codes, 
the latter being slightly better. 
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Table 1. A list of Viterbi-decoded bit error sequences caused by the 
Inner NRZ-M system as a result of weight three channel symbol error 
sequences. Hard quantization is assumed. 


JPL (7, 1/2) code ^ 

Goddard (7, 1/2) code 

Viterbi-Decoded 

Number 

V iterbi-Decoded 

Number 

Bit Error 

of 

Bit Error 

of 

Sequence 

Errors 

Sequence 

Errors 

1001000000000000 

2 

1001000000000000 

2 

1000000000000000 

1 

1000000000000000 

1 

1011000000000000 

3 

1100100000000000 

3 

1100000000000000 

2 

1100000000000000 

2 

1100100000000000 

3 

1011000000000000 

3 

1010100000000000 

3 

1110010000000000 

4 

1010000000000000 

2 

1111001000000000 

S 

1111001000000000 

5 

1110000000000000 

3 

1111110000000000 

6 

1111000000000000 

4 

1011011000000000 

5 

1101100000000000 

4 

1111000000000000 

4 

1010010000000000 

3 

1100111111000000 

8 

1010000000000000 

2 



1010100100000000 

4 


T able 2. E b /N 0 required to achieve a given P e for different schemes 


P E 


Goddard Code 



JPL Code 


NRZ-L 

Soft 

Outer 

NRZ-M 

Soft 

Inner 

NRZ-M 

Soft a 

Inner 

NRZ-M 

Soft 

NRZ-L 

Soft 

Outer 

NRZ-M 

Soft 

Inner 

NRZ-M 

Soft 3 

Inner 

NRZ-M 

Soft 

10 _1 

0.5 

0.65 

3.1 

4.0 

0.45 

0.55 

3.1 

4.0 

10 -2 

1.8 

1.9 

4.55 

5.65 

1.85 

1.95 

4.55 

5.6 

5 X 10 -3 

2.1 

2.2 

4.85 

5.95 

2.15 

2.25 

4.85 

5.95 

a Method of Eq. (2) 
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Fig. 1. Outer NRZ-M coding system 
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Fig. 2. Inner NRZ-M coding system 



Fig. 3. Block diagram for simulation of NRZ-M with convolutional codes 
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NRZ-L, HARD NO)SE 



OUTER NRZ-M, SOFT NOISE 



OUTER NRZ-M, HARD NOISE 



INNER NRZ-M, HARD NOISE 



INNER NRZ-M, SOFT 



Fig. 4. Configurations used in simulation 


























BIT ERROR PROBABILITY 



Fig. 5. Performance of Goddard (7, 1/2) convolutional code under 
various NRZ scenarios 


40 


BIT ERROR PROBABILITY 



Fig. 6. Performance of JPL (7, 1/2) convolutional code under 
various NRZ scenarios 
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Required E^/N 0 for Given Bit 
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Instead of using the criterion of maximum free distance, d f , or the maximum d f with 
minimizing a few first distance profiles, we searched short constraint length rate 1/N 
convolutional codes using a new criterion of minimizing required bit energy-to-noise 
density ratio, E h jN 0 , for a given value of desired bit error rate (BER ), for the goodness of 
a code. The considered channel was binary antipodal signaling over additive white 
Gaussian noise and no quantization at the channel output. For the BER calculations, the 
transfer function bounding technique was used. Partial searches were performed using 
some known facts and a very useful new idea that “good codes generate good codes. ” 
That is, for a given constraint length K, good rate 1/(N + 1) codes can be found by 
extending the code generator matrices of good rate 1 /N codes. The code search results are 
tabulated for 3 < K < 7 and 2 <N < 8. For many pairs of K and N, the new codes are 
shown to save 0.1 to 0.4 dB in the required EJN 0 compared to previously reported 
codes. Additionally, the benefits of coding bandwidth expansion are confirmed with our 
new codes. 


I. Introduction 

The use of a short constraint length convolutional code 
together with Viterbi decoding has been very popular for 
several applications. The class of time-invariant, nonsys- 
tematic, rate l IN convolutional codes has been studied much 
more extensively than any other class, partially due to its ease 
of analysis. The bit error rate (BER) at the Viterbi decoder 
output is well bounded by the well known transfer function 
bound (Refs. 1 and 2): 


BER < c • T(D, Z) 
O dZ 


D = D , Z = 1 
o 


( 1 ) 


where c Q and transfer function T(D, Z) depend on the code 
and the type of channel used. D 0 , which has a value between 0 
and 1, is the union-Bhattacharyya distance of the coding 
channel (everything inside the encoder-decoder pair). The right 
hand side of Eq. (1) is often represented as a series expansion 
(Refs. 1 and 2) as 
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BER < c • V a. • D‘ 

O l O 




( 2 ) 


G(n) = (G(n,l),...,G(n,k) G(n, K)), 

n=l,2,...,N (3) 


where d f is the free distance of the code. 

A large number of good rate l IN convolutional codeshave 
been found and reported (Refs. 3, 4, 5, etc.). Up to now, all 
the researchers have used the maximum criterion or the 
criterion of maximum d y together with minimizing the first 
few a- s, for determining goodness of a code in their code 
search procedures. These criteria are very valuable, if the 
values of operating D 0 ’s are much smaller than 1 , since in such 
cases Eq. (2) can be well approximated by the first few terms. 
Hence, we may consider the previously reported codes to be 
good when the required BER is extremely small. However 
when the required BER is in the moderate range of 10~ 2 to 
10 -6 , such criteria may not be good since much more than 
the first few times are required for a good approximation to 
Eq. (2). 

For uses with systems which require BER in the moderate 
range, we searched for good rate 1 IN convolutional codes by 
using a new criterion of minimizing the bit energy to noise 
density ratio, E b /N , required for a desired value of BER. For 
the evaluation of BER, we directly used the transfer function 
bound, Eq. (1). In all cases considered we assumed the use of 
binary antipodal signaling over the additive white Gaussian 
noise (AWGN) channel with no quantization at the channel 
output. For the desired values of BER, we picked 10~ 3 and 
10 "«. 


In the next section, the code structures and the transfer 
function bounding techniques are briefly reviewed for famili- 
zation with our notation. Then the partial code searching 
techniques used in this study are presented, where the very 
important idea that “good codes generate good codes” is 
introduced and explained. The code search results are given for 
3 < K < 7 and 2 KN K 8 and are compared with previously 
reported codes. As an example, for the K = 7 and N = 4 case, 
our new code can save 0.4 dB in required EJN 0 compared to 
any previously reported code when the desired value of BER is 
10 -6 . We also confirm the benefits of coding bandwidth 
expansion with our new codes. 


II. Preliminaries and Notations 


Notice that in many reports (e.g., Refs. 4, 5), the code gen- 
erator G is represented by ( G(i ), . . . , G(n), . . . , G(N)) and 
G(n)’s are in octal forms. We call this “a regular representa- 
tion” of the code generator. For later use, we also define 

g(k) = (G(i, G(n, k),... , G (N, k)), 

k =1,2 K (4) 

The «th bit in tth output vector y* (see Fig. 1) for n = 1, 2, 
. . . , N and t = 1 , 2, . . . is given by 

K 

y' n = £ G («> k) • x ( ~ k+l (5) 

k = l 


where g is a mod-2 summation, x { e{0,l}, t= 1, 2, . . . and 
x f = 0 for t < 1 by convention. The “present state” at time t, 
S { , is defined as 


S t = (x t ~ K+ 1 jc f -i) (6) 

Notice that there are 2 K ~ 1 (= M) distinct states for any value 
of N. 

To find code transfer functions, one often uses state 
diagrams where nodes represent states and directed branches 
represent state transitions. The metric value on the directed 
branch from a state to another state (or to itself), assuming the 
existence of such a transition, is the product of D to the power 
of the Hamming weight of the output vector and Z raised to 
the Hamming weight of the input bit, when a binary input 
channel is used. 

As an illustration, a K = 3, N = 2 convolutional encoder, 
with one of the best code connections, is depicted in Fig. 2. 
For this code, 


G = 


1 1 1 
1 0 1 


G(l) = (1, 1, 1), G( 2) = (1,0, 1), 
g(\f = (1, 1), g( 2) = (1,0), 


A typical nonsystematic, constraint length K, rate 1/7V SO) = (1, I) 

convolutional encoder is shown in Fig. 1 . The code connection 

box is often represented by an jVX K binary matrix G, which and the regular representation of this code generator is (7,5). 

is called code generator matrix. For a given pair of K and N, The number of states, M, is 4, and its state diagram is shown 

this code generator matrix G determines the performance. Let in Fig. 3. 


42 


~T{D,Z) = 


The transfer function T(D, Z ) can be represented (Refs. 2 and 
and 6) as 


T(D, Z) = B ■ (I - A)~ x • C (7) 


D 5 

(1-2 DZ) 2 


where / is the (M- 1) X (M- 1) unit matrix. The (M - 1) X (. M 
- 1) matrix A, the (M - l)-dimensional row vector B, and the 
( M - l)-dimensional column vector C can be obtained from 
the state diagram as 

A (/, /) = metric value on the branch from state i 
to state /, if there is such a transition 

= 0, otherwise 

B(j) = metric value on the branch from state 0 
to state /, if there is such a transition 

= 0, otherwise 

C(i) = metric value on the branch from state i 
to state 0, if there is such a transition 



The purpose of this report is to find good codes which 
minimize the required EjN 0 for desired BER = 10~ 3 or 10~ 6 , 
assuming the use of binary antipodal signaling over an AWGN 
channel and no channel output quantization. For this special 
coding channel, D a and c Q in Eq. (1) are given by (Refs. 1, 2) 

D o = exp (- EJN q ) (10) 


and 


c o = Q (2- d f - EjN o ) exp (ci f • EJNJ (11) 


where N Q is the one-sided noise power spectral density and E s 
is the received signal energy per channel bit, 

E s = r-E b (12) 


= 0, otherwise 


(= EjN , in our cases) and 


where neither state i nor state / is the all zero state (state 0). 
For Eq. (1), we need (Refs. 2 and 6) 


Q(w) = 



exp (- t 2 /2) • dt I \/2n 


(13) 




In the next section the techniques used for the partial code 
searches are described. 


+ B-(I-A)~ l ■ |§ -(I- AT 1 -C (9) 


For the previous example, 



0 D DZ 


0 

A = 

ZOO 

, B = [D 2 Z 0 0], C = 

D 2 


0 D DZ 


0 


(I-Af l 


1 

1 - 2JDZ 


1 - DZ 
Z( 1 - DZ) 
DZ 


D 

DZ 

1 - DZ 

DZ 2 

D 

1 -DZ 


T(D, Z) = 


P 5 Z 

(1 - 2 DZ) 


III. Code Searching Techniques 

For rate 1 /TV convolutional codes of constraint length K, 
the number of possible code generators is 2 K ' N . For example, 
for K = 6 and N= 4, the number is 2 24 which is over ten 
million. Furthermore, because of our new criterion, to test the 
goodness of code generators, we have to perform matrix 
inversions which require considerable amounts of computing 
time. These make exhaustive searches prohibitively difficult 
except for cases of very small K and/or very small N. For 
moderate sizes of K and N, partial searches are necessary. In 
order to find good codes with partial searches, some tech- 
niques are required for reducing the code generator space 
effectively. 

First, note that changes in the orders of G(n )’ s will not 
change the state diagram at all (Fact 1) since the Hamming 
weights of the output vectors do not depend on the orders of 
their elements. When we use this fact, the required search 
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space for code generators reduces by roughly a factor of A 
factorial. Also observe that reversing the order of the g(ky s 
gives the same performance (Fact 2) since this is equivalent to 
just redefining the state in reverse order, i.e., S* = (x f_1 , . . . , 
x f ~ K+1 ), instead of Eq. (6). This allows one to reduce the 
search space by roughly another factor of two. Notice that we 
can use these facts .for reducing the search space with no loss 
in the chance of finding the best code. There are some other 
known facts on equivalence relations between code generators. 
However, in our partial code searches, only Facts 1 and 2 were 
used since the others are much more difficult to employ. 

Recently, another useful observation was made in Ref. 7. 
That is, every known good convolutional code satisfies the 
following condition (in Ref. 7, it is called “flow conserva- 
tion”); for each state in the state diagram, summation of the 
exponents of D in all incoming branches must be equal to sum 
of those in all outgoing branches. In Fig. 3, we can easily 
check that the code satisfies this condition. We noticed that 
this condition is automatically satisfied in our cases of rate 
1/A if 

g(i) = g(K) = (1, 1, .... 1) (14) 

Notice that every reported good code in this class also satisfies 
Eq. (14). Hence we restricted our code searches to only those 
codes which satisfy Eq. (14). With this restriction, there may 
be some possibility of losing the optimum code. However, the 
use of Eq. (14) further reduces the required size of code search 
space by roughly a factor of 2 2N . 

For a given pair of K and N (at least one of them is small), 
we established the code space to be searched using Eq. (14) 
and Facts 1 and 2. Then catastrophic codes were deleted, as 
were codes for which d y was too small. We considered the 
value of df to be too small if 

d f < d fm - r (K • AO./IO 1 (15) 

where d^ m is the maximum free distance of (K, 1 IN) 
convolutional codes and r x" 1 is the smallest integer that is 
greater than or equal tox. For an example, consider the/f = 4 
and N = 3 case. The size of the reduced code space is 13 by 
using Eq. (14) and Facts 1 and 2, while the size of the original 
space is 4096 (=2 4 ‘ 3 = 2 12 ). After deleting 5 catastrophic 
codes, the following are the code generators to be considered 
further, in regular representations with a natural ordering of 
the largest element to the left: 

(17,17, 15), (17, 15, 15), (17, 15, 13), (17, 15, 11), 

(15,1 5, 13), (15, 15, 11), (15, 13, 11), and (15, 11, 11) 


In this case, d^ m is 10. Hence (15,1 1,1 1) is deleted since its <3^- 
is 7. 

For the pair of K and N, we estimated two values of EJN 
at which we believed the best code(s) could achieve the desired 
values of BER. Then for each of the remaining codes, we 
calculated the values of BER using Eq. (1) at those two values 
of E b /N 0 . Then we listed the codes with an ordering of the 
best one to the top. We consider a code to be better than 
another if the sum of the two common logarithmic values of 
calculated BER’s is the smaller. Such a listing is illustrated in 
Table 1 for the above example of the K = 4 and N = 3 case. 

By examining these listings for several values of A with a 
given (small) K , we made a very useful observation that good 
codes generate good codes. That is, for a given K, we can find 
good rate 1/(A + 1) codes by extending the code generator of 
good rate 1/A codes. For example, from the best (4,1/2) code, 
(17,15), we can obtain (1 7,1 7,1 5), (1 7,1 5,15), (1 7,1 5,13) and 
(17,15,11). Notice that all of these (4,1/3) codes are shown in 
Table 1 . For more insight on this idea, still for K = 4, some of 
the upper parts of the listings for A = 2 to 5 are shown in 
Fig. 4, with lines indicating that the left codes (with smaller A) 
generate the right codes. We limited the number of listings in 
Fig. 4 so as not to complicate the figure. Since good codes are 
generated by good codes, we do not have to use all the rate 
1/A codes for the generation of good rate 1/(A+ 1) codes. For 
the example of K = 4, the 5 best rate 1/3 codes in Table 1 are 
generated by the 2 best rate 1/2 codes. Also, the 8 best rate 
1/4 codes are obtained from the 4 best rate 1/3 codes, and so 
on. Hence when we use this idea, we can reduce the size of the 
code space very effectively. 

This idea that “good codes generate good codes” can be 
supported by the following fact: If a code# is generated by a 
noncatastrophic code A with free distance d f A , then code B is 
also noncatastrophic and its free distance, d JB , is greater than 
dj? A by at least 2. The reason for this fact is as follows: The 
value of the exponent of D on a branch in the state diagram 
for code B is always greater than or equal to the value of the 
exponent of D on the same branch in the state diagram for 
code A. Hence the code B cannot be catastrophic, since the 
code A is noncatastrophic. Additionally, the values of the 
exponents of D on the departing branch from state 0 and on 
the incoming branch to state 0 for code B are greater than 
those for code A by one [see Eq. (14) and recall that the rate 
of code B is 1 /(A + 1 ) while that of code A is 1 /A] , Therefore 
d fB - dj A is always greater than or equal to 2. 

Before presenting our code search results, we mention that 
we used the following approximation for the matrix inversion 
in Eq. (9); 
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oo p 

( i-A r 1 = Yj a ™ 88 S Am (16) 

m- 0 m-0 

and we picked the stopping number v such that 


T. y A m (i,j) < 10 5 for any m > v (17) 
i /' 

Notice that the left hand side of Eq. (17) approaches 0 mono- 
tonically as m increases. With the approximation of Eqs. (16) 
and (17), we have 5 to 7 digits of accuracy for most values of 
E b /N 0 of interest. 

IV. Results, Discussion, and Conclusions 

Table 2 summarizes our code search results which give the 
best codes in the sense of minimizing the required EjN 0 for 
desired BER = 10~ 6 and 10 -3 (among the searched codes). 
Also, previously reported codes (in Refs. 3, 4 and 5) are 
compared. Note that some of the previously reported codes 
are shown to be the best with our criterion as well. In 
particular, the codes reported in Ref. 3 are very good since 
only N =2 and 3 cases were considered and the criterion of 
maximum dj with minimizing the first few a- s in Eq. (2) was 
used. Note also that, as the code rate gets smaller (or N gets 
larger), the importance of maximum dj- diminishes, since the 
value of D a gets larger for a given value of EjN 0 [see Eqs. 
(10) and (12)]. An interesting point is that the (7,1/3) code 
reported in Ref. 3 happens to be the best for the desired BER 
= 10~ 6 , despite the fact that the larger d y codes were acci- 
dently overlooked (as noted in Ref. 4). In Fig. 5, BER versus 
E b lN 0 plots are given for the K = 7 cases. 

Using our new codes it is possible to save up to 0.4 dB in 
required EjN Q compared to previously reported codes [see 
the (7,1/4) case]. To see the gains visually, we plot the 
required E b /N () for BER = 10“ 3 and 10 -6 versus N curves in 
Fig. 6. Note that for our new codes, the required E b /N () 
becomes smaller as the code rate gets smaller (Fact A) for a 
given (large) value of K. This fact A is known as the gain due 
to coding bandwidth expansion. That is, we can save in 
required signal energy at the expanse of using the larger 
bandwidth associated with a lower code rate (with similarly 
complex codes 1 ). Note that for the previously found codes, 
this fact cannot be observed, since the codes are not 


'The complexity of ( K , 1 /AO convolutional code increases as K 
increases and also as N increases. But since K is a much more 
important factor for complexity than N , we will consider the codes 
with the same K to have similar complexities for any N, 


necessarily optimum in the sense of minimizing the required 
EjN 0 . Note also that for small values of K (3 and 4) this fact 
cannot be seen even with our codes. Therefore, we conjectured 
that the benefit of coding bandwidth expansion can be 
obtained only when enough complexity is allowed. 

Now consider the capacity and cutoff rate for further 
insight into Fact A. The capacity, C, in the dimension of code 
rate [information bits/channel bit] of this coding channel 
(binary antipodal signaling over AWGN, no channel output 
quantization) is given by (Ref. 2) 

C = [- ( 1 + In 2n)/2 - fp(y ) • In p(y) • dy] /In 2 

(18) 

where 

p{y) ~ [exp (- (y - b) 2 / 2) + exp (- (y + b) 2 12)] /V 8tt 

(19) 

and 

b 2 ! 2 = EjN o (20) 

The channel coding theorem (Ref. 2) says, with coding, 
arbitrary small error rate can be achieved provided that the 
data rate R b [information bits/sec] is less than the channel 
capacity [information bits/sec] of the coding channel. 
With dimensions of code rate, the above condition can be 
restated as 

r < C (21) 

That is, the choice of code rate must be smaller than C, in 

theory. Note that C is a function of EjN 0 , and C (EJN 0 ) is a 
monotonically increasing function of its argument. Hence its 
inverse function exists. Therefore Eq. (21) can be rewritten as 

EjN o >r l (r) (22) 

With Eq. (12), 

EjN o > C~' (r)lr (23) 

We call the right hand side of Eq. (23) “the required EjN 0 to 
achieve capacity.” These values were found and plotted in 
Fig. 7 as a function of coding bandwidth expansion factor (or 
1/code rate r). Also the cutoff rate of the coding channel R a 
[information bits/channel bit] is given by (Ref. 2) 

R (> = 1 - log 2 (1 + exp (- E JN o )) (24) 
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It is often said (e.g., Ref. 8) that a small enough error rate can 
be achieved, with practical coding , provided that R h is less 
than the cutoff rate Q in [information bits/sec] . Or, 
equivalently, 

' < R o <W (25) 

This can likewise be inverted to yield 

EJN q >R- o 1 (r)/r = - j • In [ 2 - 1] (26) 

We call the right hand side of Eq. (26) “the required EJN Q to 
achieve cutoff rate.” These values are also plotted in Fig. 7. 
These theoretical curves tell us that to be in the region below 
the capacity curve is impossible in theory; to be in the region 
between the capacity and the cutoff rate curves is possible in 
theory but difficult in practice; to be in the region above the 
cutoff rate curve is practically possible. Note that from both 
the capacity and the cutoff rate curves, we see that lowering 
the code rate gives the benefit of reducing the required E b /N Q , 
which is the theoretical view of Fact A. 


In the same figure, we have plotted the performances of our 
new codes for K = 5, 6 and 7 (from Fig. 6). From comparisons 
of the slopes of the curves for BER = 10~ 3 with those for BER 
= 10“ 6 cases, we conclude that the benefit of coding 
bandwidth expansion is greater when the desired BER is larger. 
Also, increasing the system' complexity (i.e. , increasing K) 
provides more gain when the desired BER is smaller, as seen by 
the larger spacing between the curves for the BER = 10“ 6 
cases compared to the BER = 10~ 3 cases. Note also the 
similarity between curves for our new codes with those derived 
from theory (especially the cutoff rate curve). 

In conclusion, we have found good ( K , l/N ) codes which 
minimize the required E b /N 0 for the desired BER = 10~ 3 and 
10 -6 for 3 < K < 7 and 2 <JV < 8. For the partial searches of 
codes, we used some known facts together with another very 
useful idea that “good codes generate good codes.” For many 
pairs of K and N, our new codes are shown to save 0.1 to 0.4 
dB in required E b /N Q for moderate values of required BER, 
compares to the previously reported codes. Also, we con- 
firmed the benefits of coding bandwidth expansion with our 
new codes, whereas the previously reported codes did not 
uniformly confirm this property. 
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Table 1. Examples of listing of code search results for the 



K 

= 4 and N = 

3 case 


G in octal 

d f 

Iog 10 (BER) 
at E b /N a = 

Sum of 
the Two 


6.0 dB 

3.5 dB 

Values 

(17, 15, 13) 

10 

-6.059 

-3.070 

-9.129 

(17, 15, 11) 

9 

-6.008 

-3.082 

-9.090 

(15, 13, 11) 

8 

-5.609 

-3.014 

-8.623 

(17, 15, 15) 

10 

-5.702 

-2.567 

-8.289 

(15, 15, 13) 

9 

-5.516 

-2.510 

-8.026 

(15, 15, 11) 

8 

-5.422 

-2.490 

-7.912 

(17, 17, 15) 

8 

-5.279 

-2.488 

-7.767 
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Table 2. Code search results 


Required E b /N a , dB 


K 

N 

Code Generator G in 
Regular Representation, octal 

d f 

for desired BER= 

Notes 

1 .1-1-6 

l.E-3 

3 

2 

7,5 

5 

6.706 

4.017 

1 , a, b 

3 

3 

7, 7,5 

8 

6.736 

4.140 

1 , a, b 

3 

4 

7, 7, 7,5 

10 

6.938 

4.354 

2 

3 

4 

7, 7,5,5 

10 

6.706 

4.017 

4, a, b 

3 

5 

7, 7, 7, 5, 5 

13 

6.669 

4.055 

3, a, b 

3 

6 

7, 7, 7, 7,5,5 

16 

6.736 

4.140 

3 

3 

6 

7, 7, 7, 5,5,5 

15 

6.706 

4.017 

4, a, b 

3 

7 

7, 7, 7,7, 5,5,5 

18 

6.664 

4.035 

3, a, b 

3 

8 

7, 7, 7, 7, 7, 5, 5, 5 

21 

6.685 

4.080 

3, a 

3 

8 

7, 7, 7, 7, 5, 5, 5, 5 

20 

6.706 

4.017 

4, b 

4 

2 

17, 15 

6 

6.180 

3.735 

1 , a, b 

4 

3 

17, 15, 13 

10 

5.958 

3.437 

1, a 

4 

3 

17, 15, 11 

9 

5.994 

3.427 

4, b 

4 

4 

17, 15, 15, 13 

13 

6.004 

3.511 

2 

4 

4 

17, 15, 13, 11 

12 

5.906 

3.286 

4, a, b 

4 

5 

17, 17, 15, 15, 13 

16 

5.991 

3.483 

3 

4 

5 

17, 15, 15, 13, 11 

15 

5.909 

3.332 

4, a, b 

4 

6 

17, 17, 15, 15, 13, 13 

20 

5.958 

3.437 

3 

4 

6 

17, 17, 15, 15, 13, 11 

19 

5.865 

3.325 

4, a, b 

4 

7 

17, 17, 15, 15, 15, 13, 13 

23 

5.974 

3.462 

3 

4 

7 

17, 17, 15, 15, 13, 13, 11 

22 

5.849 

3.306 

4, a, b 

4 

8 

17, 17, 17, 15, 15, 15, 13, 13 

26 

5.972 

3.456 

3 

4 

8 

17, 17, 15, 15, 15, 13, 13, 11 

25 

5.860 

3.332 

4, a 

4 

8 

17, 17, 15, 15, 13, 13, 11,11 

24 

5.906 

3.286 

4. b 

5 

2 

35,23 

7 

5.745 

3.495 

1 , a 

5 

2 

31, 23 

6 

5.845 

3.430 

4. b 

5 

3 

37, 33, 25 

12 

5.395 

3.115 

1 , a, b 

5 

4 

37, 35, 33, 25 

16 

5.303 

2.999 

2 

5 

4 

37, 35, 25, 23 

15 

5.298 

3.000 

4, a 

5 

4 

35, 31, 27, 23 

14 

5.317 

2.965 

4, b 

5 

5 

37, 35, 33, 27, 25 

20 

5.270 

2.923 

3 

5 

5 • 

37,35, 31, 27, 25 

19 

5.243 

2.924 

4, a 

5 

5 

37, 35, 33, 25, 21 

18 

5.297 

2.912 

4, b 

5 

6 

37, 35,35,33, 27, 25 

24 

5.291 

2.957 

3 

5 

6 

37, 35, 33, 27, 25, 23 

23 

5.21 1 

2.880 

4. a 

5 

6 

37, 35, 33, 27, 25, 21 

22 

5.294 

2.867 

4, b 

5 

7 

37,35,35,33,27,27,25 

28 

5.286 

2.955 

3 

5 

7 

37, 35, 33, 31, 27, 25, 23 

26 

5.21 1 

2.845 

4. a 

5 

7 

37,35,33,27,25,23,21 

25 

5.256 

2.839 

4, b 

5 

8 

37, 37, 35, 33, 33, 27, 25, 25 

32 

5.284 

2.949 

3 

5 

8 

37,35,35,33,31,27,25,23 

30 

5.211 

2.860 

4, a 

5 

8 

37, 35, 33, 31, 27, 25, 23, 21 

28 

5.280 

2.819 

4, b 
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Table 2 (contd) 


K 

N 

Code Generator G in 
Regular Representation, octal 

d f 

Required E b /N Q , dB 
for desired BER= 

Notes 



l.E-6 

l.E-3 


6 

2 

75, 53 

8 

5.310 

3.289 

1 

6 

2 

77,45 

8 

5.236 

3.242 

4, a 

6 

2 

75, 57 

8 

5.293 

3.211 

4, b 

6 

3 

75,53,47 

13 

4.918 

2.900 

1, a 

6 

3 

75,67,41 

12 

5.034 

2.854 

4, b 

6 

4 

75,71,67,53 

18 

4.836 

2.747 

2 

6 

4 

77, 73, 55,45 

18 

4.779 

2.729 

4, a 

6 

4 

77, 73,51,45 

17 

4.807 

2.719 

4, b 

6 

5 

75,73,71,65,57 

22 

4.826 

2.680 

3 

6 

5 

77,73,71,55,45 

22 

4.742 

2.660 

4, a 

6 

5 

75,71,65,57,53 

22 

4.753 

2.645 

4, b 

6 

6 

75,73,65,57,55,47 

27 

4.764 

2.616 

3 

6 

6 

77,73,67,55,51,45 

26 

4.694 

2.591 

4, a 

6 

6 

75, 71,65,57,53,47 

26 

4.704 

2.590 

4, b 

6 

7 

75,75,67,65,57,53,47 

32 

4.762 

2.630 

3 

6 

7 

77, 73,67, 63, 55, 51,45 

30 

4.696 

2.565 

4, a 

6 

7 

75,71,65,57,53,47,43 

29 

4.717 

2.564 

4, b 

6 

8 

75,73,67,65,57,57,51,47 

36 

4.728 

2.599 

3 

6 

8 

77,73,67,63,57,55,51,45 

35 

4.693 

2.552 

4, a 

6 

8 

77,73,67,63,55,51,45,41 

32 

4.772 

2.541 

4, b 

7 

2 

171, 133 

10 

4.802 

3.036 

1, a 

7 

2 

161, 133 

9 

4.818 

3.035 

4, b 

7 

3 

175, 145, 133 

15 

4.599 

2.706 

2 

7 

3 

171, 145, 133 

14 

4.489 

2.672 

1, a 

7 

3 

161, 135, 107 

13 

4.600 

2.666 

4, b 

7 

4 

163, 147, 135, 135 

20 

4.761 

2.814 

2 

7 

4 

175, 151, 133, 117 

20 

4.372 

2.520 

4, a 

7 

4 

173, 167, 135, 111 

20 

4.433 

2.511 

4, b 

7 

5 

175, 147, 135, 135, 131 

25 

4.562 

2.653 

3 

7 

5 

175, 151, 133, 127, 1 17 

25 

4.310 

2.444 

4, a 

7 

5 

175, 165, 151, 133, 117 

25 

4.350 

2.438 

4, b 

7 

6 

173. 163, 151, 137, 135, 135 

30 

4.394 

2.499 

3 

7 

6 

175, 171, 151, 133, 127, 117 

30 

4.286 

2.387 

4, a 

7 

6 

175, 165, 151, 137, 133, 117 

30 

4.307 

2.385 

4, b 

7 

7 

173, 165, 147, 145, 137, 135, 135 

36 

4.381 

2.456 

3 

7 

7 

175, 171, 155, 127, 123, 1 17, 1 13 

34 

4.266 

2.356 

4, a 

7 

7 

175, 171, 165, 151, 133, 127, 117 

35 

4.282 

2.353 

4, b 

7 

8 

173, 165, 153, 147, 137, 135, 135, 111 

40 

4.312 

2.387 

3 

7 

8 

175, 171, 165, 151, 133, 127, 1 17, 1 13 

39 

4.251 

2.327 

4, a, b 

Notes: 
1 : 
2: 
3: 
4: 
a: 
b: 

bound by Odenwalder (Ref. 3) 
l ound by Larsen (Ref. 4) 
bound by Daut, et al. (Ref 5) 
bound by the author 

Minimizes required E JN a for desired BER = l.E-6 
Minimizes required h'/jN () for desired BER = l.E-3 
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Fig. 3. The state diagram of the encoder in Fig. 2 


N = 5 

(17,15,15,13,11) 



Fig. 4. An example for the idea of “good codes generate good codes,” K = 4 



































Fig. 6. Required Eg /Ng versus N (= 


N (= 1/CODE RATE r) 


i rate r). Solid line indicates codes found previously; dashed line indicates codes found here 
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In implementing multifrequency PPM, a naturally arising question is: Let P be a fixed 
number; among all integer valued processes X v X 2 , X 3> . . . with E(\X n+1 ~ X n \ 2 ') 
which has the largest entropy? Earlier work by McEliece and Rodemich answered this 
question, but there is no obvious way to use this process to implement a code for 
multifrequency PPM. The present article describes an easy-to-implement process X v X 2 , 
with E(\X - X n \ 2 ) <,P, whose entropy is nearly as great as that of the 
McEliece-Rodemich process. 


I. Introduction 

McEliece and Rodemich, in a study of the feasibility of 
using multifrequency PPM for optical communication (Ref. 1), 
were led to consider the following problem. Among all 
stationary random processes {. . . , X_ , X Q , X 1 , . . .} taking 
values in the set {1, 2, . . . , JV} and subject to the constraint 

0 ) 

how large can the entropy be? McEliece and Rodemich 
showed that the maximum entropy is achieved by a Markov 
chain and gave an exact formula for the transition probabilities 
of the maximizing chain (Ref. 1). In order to apply these 
results to the original problem in optical communication, 
however, it would be necessary to encode a given binary data 
stream into a sequence of symbols from {1,2, . . . , N] which 
closely resembles a “typical sequence” from the optimal 
Markov chain. Given the complex form of the solution given in 
Ref. 1, such an encoding would very likely be difficult to 
implement. In this article we shall examine some approxima- 
tions to the solution given in Ref. 1 which are not far from 
optimal (about 0.25 bit) and which suggest practical ways to 
implement energy-efficient multicolor PPM. 


II. Approximations 

In this section we attempt to motivate the construction to 
be given in Section III, as follows. We consider a simplified 
version of the problem mentioned in the Introduction, viz, we 
enlarge the state space from {1,2, . . . , N} to Z, the set of all 
integers. For a given value of P this will increase the maximum 
possible entropy, but ithP « N 2 this increase will be small. In 
any case our object in this section is to compare the entropy 
vs. P relationship for three possible Markov chains on Z subject 
to Eq. (1): the optimal Markov chain and two others. In 
Section III we will modify one of the suboptinral chains 
discussed here to devise a practical modulation scheme for 
multicolor optical communication, whose performance is (in a 
sense to be given below) within 0.25 bits of the theoretical 
optimum. 

With the state setZ, the maximum entropy Markov chain 
subject to Eq. (1) is a Markov chain for which the increments 


are independent and identically distributed. The entropy of 
the resulting chain is just the entropy of the random variable 
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A n , and so to solve the original problem on Z we just need to 
maximize H( A) subject to the condition ^(A 2 ) This can 
be done using straightforward variational techniques (Ref. 2, 
Problem 1.8), and the maximum entropy is given parame- 
trically as follows: 

For each X > 0, define 

oo 

m (X) = 2 (3) 


The optimal distribution A n is unfortunately not well 
suited for adaptation to a practical modulation scheme. The 
exact distribution is in fact 

Pr {A = k} = e~ n x /m(X) (10) 

a nonuniform distribution on a countable set (0, ±1, ±2, . . .} 
of values. However, we can get a surprisingly large entropy by 
considering instead of Eq. (10) a much simpler random 
variable A( L ), which is uniformly distributed on {-L, -L + 1, 
..., 1 - 1 }: 


so that 


-m'(X) = n 2 e n *■ 


(4) 


Pr {A (i > = A} = 


_1_ 
2 L 


if -L <K <L - 1 


= 0 otherwise 


(ID 


Then among all random variables satisfying £"(A 2 ) the 
maximum entropy, H , is given by 


For the sequence {. . . , X_ v X Q , X x , . . . } whose increments 
X n+l - X n are i.i.d. with common distribution A^) , a simple 
calculation gives 


(5) 

where 


P = 


m '(V 

m(\) 


( 6 ) 


For small values of X the sums in Eqs. (3) and (4) are well 
approximated by the following integrals: 



Using these approximations in Eqs. (5) and (6), we obtain the 
approximation 


P = ^ + t -*) 2 

= E (A (l) 2 ) 


H = log (2 L) 

Thus for this particular Markov chain, the relationship 
between the entropy H and the “power” P is 

H = j logP+ j log 12+ l - log ( 1 - ~) (14) 


( 12 ) 

(13) 


Comparing Eqs. (9) and (14), we see that the difference in 
entropy between the optimal distribution of increments (Eq. 
10) and the suboptimal distribution (Eq. 1 1 ) is approximately 
1/2 log ne/6 = 0.255 bits. 


^max ^ 2 l0g/,+ 2 l0g27re ’ ' PLarge ( 9 ) 

which is extremely good even for small values of P, as 
exhibited in Table 1. In Table 1 we have tabulated, for a range 
of X’s, the corresponding value of P calculated from Eq. (6), 
the exact value of H calculated from Eq. (5), and the 
approximate value of H m&x from Eq. (9). We conclude that 
for P~> 1, there is no significant difference between the exact 
value of // max given by Eq. (5) and the approximation given 
by Eq. (9). 


As a comparison, we consider the Markov chain {. . . , X_ lr 
X Q , Xy . . .} in which the components X t are i.i.d., uniformly 
distributed on (1,2, . . . , /,}. In this case it is easy to calculate 
[cf. Eqs. (12) and (13)] , 

L 2 1 

/>=--! (15) 

H = log L (16) 
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from which follows 

H = |log/> + ^ log 6 + -^+0(P- 2 ) (17) 

Comparing Eqs. (17) and (14), we see that for a given value of 
P, a Markov chain with uniform and independent increments 
gives one-half bit more entropy than one whose components 
are independent and uniform. For a fixed value of H, the 
distribution (Eq. 11) requires about 1.53 dB more power than 
the optimal distribution, whereas a uniform distribution on 
the X’s requires 4.53 dB more power. 

Motivated by these results, in the next section we introduce 
an encoding process which maps a sequence of 0 ’s and l’s into 
a sequence of elements from {-L, -L + 1, . . . , L - 1} that 
closely resembles a typical sequence from Markov chain whose 
increments are described by Eq. (1 1). 


III. An Encoding Algorithm 

Motivated by the results of Section II, we propose a 
method of encoding a random stream of 0 ’s and l’s, say U v 
U v U 3 , . . . into a sequence X v X 2 , . . . of elements from the 
set {-N/2, -Nj2 + 1, . . . ,0,1, . . . ,N/2~ l},such that Eq. (1) 
is satisfied . 1 


The random variables Aj A 2 , ... are i.i.d., each being 
uniformly distributed between ~2 W_1 and 2 H ~ 1 - 1. Hence if 
we define the Markov chain {X' m } m>0 by 


*o =° 


X' = X' + A _ 

m+1 m m + 1 


( 21 ) 


it follows from the results of Section II that the resulting 
Markov chain has [cf. Eq. (12)] 


E(\X ' ~ X', 

v m m+1 


I 2 ) 


( 2 h ~ 1 ) 2 

2 


+ 


6 



( 22 ) 


by Eqs. (18) and (19). The entropy of {X^ } is H bits. Thus 
from Section II, if P is large, the difference between the 
entropy of this chain and that of the optimal Markov chain is 
only about 0.255 bits. 

However the chain {X ' m } is not satisfactory, since there is 
no guarantee that X' m will lie between -N/2 and N/2 - 1. The 
following definition remedies the situation. 


First choose L to be the largest integer such that 

T + 6 <P < 18 > 

and H to be the largest integer such that 

2 H ~ l < L (19) 

In our encoding algorithm the symbol X m + l will be deter- 
mined by X m and the ( m + l)-st block of // data bits, viz 
[U Hm +i> ■■ •> U H(m + i)l' These H bits are in fact used to 
determine an integer A m + 1 in the range [-2 H ~ l , 2 H ~ l - 1] 
via two’s complement arithmetic. For example, 

H 

Aj = £ u k 2 "~ k - u i 2 H (2°) 

k = t 


*For notational convenience, in this section we assume N is even, and 
let the state set be as described, instead of {l ,2, . . . , N} . 


Jr = min (X ,B) + A _ if X >0 

(23) 

= max (X ,-B- 1)-A if X <0 

where 

B =■ 7 - 2 W_1 (24) 

6 

The number B defined in Eq. (24) is the largest value for X' m 
that guarantees X^ +l < N/2 - 1. Similarly -B ~ 1 is the 
smallest value for X' m that guarantees X' m+l > -N/2. It 
follows that the chain {X m } defined by Eq. (23) will lie in the 
desired range -N/2 < X m <N/2. The entropy of still 

H bits, since H(X m+l X m ) = H for all m. The value of E{\X m - 
X m+ , 1 2 ) will be somewhat larger than the corresponding 
value for { X ^}, since when X m > B or X m < -B - 1 the 
difference X m + 1 ~X m will no longer be uniformly distributed 
on {-2 W_I , . . . , 2 H ~ l - 1}. However, since E (A m ) = - 1/2 
the rule [Eq. (23)] causes the chain to be attracted to 0, and 
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unlikely to lie near the boundaries. In the Appendix we make 
this precise and show that in fact 

where r is the unique solution in (0,1) to the equation 
2"- 1 

£ zk = 2" (26) 


Table 2 gives the value of these roots for H = 1,2, . . . , 7. It 
follows that for a fixed P, if N is sufficiently large, Eq. (1) will 
be satisified by the Markov chain {X m }. 

We conclude with a simple numerical example. 


Example: Let N = 64, P = 23. From Eqs. (18), (19), (24) 
we have L = 8, H = 4, B = 24. The data sequence u = (1000 
01 10 01 1 1 01 10 1 100 . . .) yields the increment sequence Aj 
= -8, A 2 = 6, A 3 = 7, A 4 = 6, A s = -4, and so by Eq. (23) we 
have X Q = 0 and (X t , X 2 , . . .) = (-8, - 14, -21, -27,-21, . . .). 
The entropy of this chain is// = 4 bits, and from Eq. (25) and 
Table 2 

^d^ + i'^J 2 ) O+(0.953817) 24 ) + j = 28.36 

In fact, an exact calculation of the steady state probabilities 
for this chain shows that P = 22.6. As comparison, we note 
that Eq. (9) shows that the largest possible entropy for a 
Markov chain with P = 22.6 is 4.3 bits. The performance of 
our algorithm in this example is very close to the 0.255 bit loss 
predicted in Section II. 


References 


1. McEliece, Robert J., and Rodemich, Eugene, “A Maximum Entropy Markov Chain,” 
Proc. 17th Annual Conference on Information Sciences and Systems, Johns Hopkins 
University, Baltimore, MD (1983), pp. 245-248. 

2. McEliece, Robert J., The Theory of Information and Coding, Encyclopedia of 
Mathematics and Its Applications, Addison-Wesley Publishing Company, 1977. 

3. Kleinrock, Leonard, Queuing Systems, Volume 2: Computer Applications, John Wiley 
and Sons, Inc., 1976. 


60 



Table 1. Comparison of H max to its approximation 


X 

P 

H (exact) 

max v ' 


1.00 

0.499 

1.0715 

1.0713 

0.95 

0.526 

1.0974 

1.0974 

0.90 

0.5551 

1.1247 

1.1247 

0.85 

0.5880 

1.1534 

1.1534 

0.80 

0.6249 

1.1838 

1.1838 


Table 2. Roots of Eq. (26) 


H 

r 

2 

0.414214 

3 

0.823408 

4 

0.953817 

5 

0.988325 

6 

0.997073 

7 

0.999268 



Appendix 


This appendix uses Kingman’s bound, below, to bound the 
value of E(X m+1 - X m ) 2 , where X v X 2 , . . . is the process of 
Section III. This is possible because Kingman’s bound gives a 
maximum possible value to P(X m >£) for |£|>.fi, while we 
have 


P(X m >V 


^ (N/2 - fi) « 

< —' r 


where 0<r< 1 and 


E " X J 1 1 ♦»« 



/=-a:+i 


Throughout this appendix, K will mean 2 H 1 , and N is an 
even integer larger than 2 K. 

Lemma 1: Let Y v Y 2 , ... be i.i.d. random variables, P(Y n 
= £) = 1/2 K, - K < ft <K. Let W Q = 0, and let W n+1 = max (0, 
W n + Y n+1 ). Let r, 0O< 1, be a root of 

t Z '~ 2K 

i=-k+i 


Proof: Given that X Q = 0, X v X^, ■ ■ ■ has exactly the same 
distribution as T v T 2 , . . . in lemma 2. Thus/ > (Tf m >5) 

Separately, for B < 2 <N/2, 

, tv/2-i 
j = 2 ~K + 1 

So 


Then P(X m > ft) < r s . Lemma 1 is a special case of Kingman’s 
bound (Ref. 3, p. 44). 


Lemma 2: Let Y v Y 2 , . . . be as in lemma 1. Let T 0 = 0 
and B = N/2 - K, and define T v T 2 , . . . by 

r +1 = miner ifT n >0 

= max(T„,-Z?-l)-r n+1 if 7 < 0 

Let r, again, satisfy 0 < r < 1 and 


P(X m =B)>P(X m =B+ >P(X m = N/2 - 1) 


Therefore 


— T r> = l 
2 K 


i=-K + i 


P(X >ft)< ■ r B 

K ' K 


Corollary: E(X m + , - Z m ) 2 <K 2 / 3 + 1/6 + (4 K 2 /3)r B . 
Proof: As shown in Section II, 

£■((* +1 - X ) 2 kfi- 1 <B)=22-+ 1/6 

vv ra + 1 m 7 I m 7 3 

Therefore 


Then for all n, and all £, we ha \eP(T n > £) <r e . 

Proof: The statement is trivial for £ <0. For any sequence 
Y v Y 2 , . . ., T n < IP^, where IP,, is the process of lemma 1, 
and so the statement is true for £ >0. 


E(X 


m + 1 



+ 


6 


+ 2 


Nf 2-1 

z 


e=s 


/>(* =£) 

K m ' 


Theorem: For the process X v X 2 , . . . described in 
Section III, the steady state probability P satisfies, for ft^B, 


X E( ( X n 


X y \X =£) 

m m ' 
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But 


At/ 2-1 


2 £ !*,=«) 


fi = B 


< 


/V/2-1 B / , */2-l 

2 E — '(k E «-# 

e=B \ <=8-* 


B At/2- 1 At/2-1 

7 ; £ £ 

A £=£ /=5-/C 


_ (2* + 1) (2K - 1) B 
3 


, 4 K 2 i 

< r 

3 
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A microprocessor based antenna controller system developed as part of the unattended 
station project for DSS 13 is described. Both the hardware and software top level designs 
are presented with a disucssion of the major problems encountered. Developments useful 
to related projects include a JPL standard 15 line interface using a single board computer, 
a general purpose parser, a fast floating point to ASCII conversion technique, and experi- 
ence gained in using off board floating point processors with the 8080 CPU. 


I. Introduction 

The design of a microprocessor antenna controller was per- 
formed as part of the DSN unattended station development at 
DSS-13. The objectives were to replace the existing antenna 
controller (a Modcomp II) with a low cost standard 8 bit mi- 
croprocessor controller similar to the microprocessor control- 
lers used in the other automated subsystems, and to replace 
the current assembly language software with a high level pro- 
gramming language utilizing a top-down structured design. 
Another objective was to replace the antenna data collection 
system with a similar design. 

II. System Description 

A. Hardware 

The system consists of two controllers as shown in Fig. 1. 
The controllers utilize only commercially available modules to 
simplify requirements for documentation and spares. Collect- 
ing analog and digital data from the various remote transducers 
and sending this data to the main controller located in the 
control room is the function of a data collection controller 
located in the antenna pedestal. The main controller collects 


data from local monitor points in the control room, provides 
operator interface, communication, display, and necessary 
logic and calculations for controlling the antenna. Communi- 
cation to the station controller is through a Star Switch 
utilizing a JPL standard 15 line interface. 

The remote data collection controller configuration is 
shown in Fig. 2. It consists of a single board computer (8080 
CPU) with 8 K-bytes of onboard PROM and 4 K-bytes of 
RAM. The CPU drives a short haul modem through a serial 
port and is connected to 24 digital channels via onboard paral- 
lel ports. An external digital interface uses opto-isolator 
modules to convert 24 volt antenna levels to T.T.L. levels. A 
32 channel differential input multiplexer and analog to digital 
converter provide interface to the analog channels. All input 
lines are protected from lightning discharges by filters as well 
as varistor type surge protectors. The following development 
modules, RAM module, floppy disk controller, communica- 
tion expansion module, and dual floppy disk drives allowed 
software development and testing in the remote data collec- 
tion controller, thus eliminating the need to use costly cross 
compilers and/or development systems. After checkout the 
program was burned into PROM and these development mod- 
ules were removed, leaving only the CPU and analog to digital 
converter cards. 
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The main controller configuration is shown in Fig. 3. Two 
CPUs are utilized on the same Multibus. CPU 1 provides inter- 
face to the angle encoders through its parallel ports and per- 
forms most of the tracking related calculations. CPU 2 pro- 
vides interface to the Star Switch through its parallel ports and 
the serial data from the data collection controller and modem 
through its serial port. CPU 2 performs the conversion of data 
to floating point and ASCII. It also handles the display and 
operator interface. In addition, CPU 2 performs the major part 
of the tasking logic and initiates appropriate action in the 
event of a fault detection. 

Analog and digital interfaces are provided in a similar man- 
ner as the remote data collection controller, except that an 
input-output expansion module is necessary since the parallel 
ports of both CPUs are utilized. Numerical computations are 
facilitated by three floating point processors dedicated to 
tracking computations, digital servo-loop computations, and 
data conversion tasks such as binary to floating point or 
ASCII. A video-keyboard module provides interface to the 
local maintenance operator keyboard and high resolution 
monitor. Near real time status and performance data are 
displayed in the upper portion of the maintenance screen 
while response messages and operator inputs are displayed in 
the lower portion of the screen. Upper screen data that is not 
within expected limits, or states, is displayed for the mainte- 
nance operator in reverse video and flagged within the program 
for possible automatic monitor action. A communications 
expansion module provides interface to a printer for mainte- 
nance data logging, or for making listings of the program 
during development. It also provides an interface for an 
additional terminal which allows separate control of each CPU 
for debugging purposes. A digital to analog converter module 
provides the output drive voltages to operate the antenna 
servo-hydraulic system. As with the remote data collection 
controller, software was developed in the controller using a 
floppy disk controller and two disk drives. 

B. Software 

The program is written in the PL/M (an Intel trademark) 
language. Currently the main controller program contains 
about 10,000 lines of source code with 200 lines of assembly 
language for high speed peripheral drivers. The compiled pro- 
gram requires 54 K-bytes including data storage. For compari- 
son, the minicomputer program contained 2,300 lines of 
FORTRAN and 4,600 lines of assembly language requiring 30 
K-words of memory (a Modcomp word is 2 bytes). 

The top level flow chart for the remote data collection con- 
troller is shown in Fig. 4. All of the flowcharts shown have 
been simplified for illustration purposes, and some modules 
may actually represent a combination of several. After an ini- 


tialization module, the program enters an endless loop contain- 
ing the routines REMDAT, which inputs data from the trans- 
ducers, and REMSEND, which formats, adds a checksum, and 
sends the data over the serial link to the main controller. A 
timer on board the CPU provides a 50 ms interrupt which is 
divided modulo 4 so that the main loop is executed every 
200 ms. The serial line operates at 9,600 baud. The data is 
formatted in two 73 character blocks of ASCII characters 
which can be displayed directly on a local maintenance termi- 
nal if desired for interface verification. 

The main controller top level software is shown in Fig. 5. 
Both CPUs begin execution of the same code from exactly 
the same memory location. Both have the ability to interro- 
gate themselves to find out which CPU they are and jump to 
their separate initialization routines. From there each enters a 
separate endless background loop except for the servicing of 
their interrupts. 

The level 2 flowchart for CPU 1 is shown in Fig. 6. The 
background modules perform the following functions: 


CHKSUM 

Verifies that the code has not changed 

SCHEDULER 

Looks at current time and initiates any 
tasks that may be scheduled 

CHKSTAK 

Generates a warning message if the stack 
is getting too large 

There are three interrupts used with CPU 1 . Interrupts 3 
and 4 are driven from the station clock at rates of 50 Hz and 
1 Hz, respectively, which control all the time dependent mo- 
tion of the antenna. Interrupt 2 is used to read a character of 
serial time from the station clock. The station clock is read 
initially to establish time. Thereafter the antenna clock is in- 
cremented in the 1 Hz interrupt routine. The 50 Hz routine 
counts each interrupt to insure exact synchronization with the 
1 Hz interrupt. If the count is not correct or if the antenna 
time does not match station time, an error message is gener- 
ated. The count is also used modulo 5 to call the following 
modules every 100 ms: 

READOUTS 

Inputs current position from the angle 
encoders 

UPDATE 

Updates the angle commands to the cur- 
rent time 

SERVO 

Closes the digital servo loop (type 1 or 
type 2 depending on start up conditions) 

D2AOUT 

Updates the current output drive voltage 
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After the interrupt 4 routine (1 Hz) steps the clock, the fol- 
lowing modules are called: 

ANAACTCNT Decrements the counters which delay 
action on analog channels until after 
specified delays 

TIMERS Decrements action timers 

CMPRATES Computes the angle rates by differenc- 

ing the position readouts 

There are three tracking modes and an idle mode. An azi- 
muth elevation mode allows pointing to a given azimuth eleva- 
tion angle pair. A sidereal rate mode is used for tracking ob- 
jects at sidereal rate. A three day fit mode is used for tracking 
objects such as spacecraft that are near sidereal rate and uses 
a second order fit of three values of right ascension and decli- 
nation. Depending on the mode, one of the routines AEPRED, 
SID, or FIT3 is called to compute the current pointing angle in 
azimuth elevation coordinates. Three more routines are called: 

CHKPRELIM Checks the physical movement limits of 
the antenna and stops motion if they are 
about to be exceeded 

REFC Computes a refraction correction using 

the corrected Berman model 

WINDAVG Calculates the average wind based on the 
last 2 minutes; if the average wind ex- 
ceeds its limit the antenna will automati- 
cally stow 

The CPU 2 level 2 software is shown in Fig. 7. The back- 
ground routines perform the following functions: 

SPEEDSELECT Selects either high or low speed mode de- 
pending on the distance to the desired 
point. An azimuth distance greater than 
10 degrees or elevation greater than 
3 degrees will select high speed. 

CONVERT Inputs the control room monitor chan- 
nels and converts all channels (including 
those from the data collection controller) 
to floating point and ASCII. This module 
also checks limits and writes data to the 
DATAOUT buffers. If monitor action is 
required, flags are set for the TASKREQ 
module. 

MONACT Performs the actual tasks selected by the 

logic in TASKREQ as the result of a mon- 
itor action. 


TSKACT Performs the actual tasks selected by the 

logic in TASKREQ as the result of an 
operator or scheduler input. 

DATAOUT Outputs data from the screen, Star 
Switch or printer buffers. 

Four interrupt routines handle the receipt of a character from 
the data collection controller, a Star Switch message, an opera- 
tor keyboard character, or a Star Switch message time out. 

Five help menus are available at the local display, the sta- 
tion controller display, or the remote display in the NOCC to 
provide the operator with assistance in remembering the com- 
mands available and their syntax. They are shown in Fig. 8. 
All operator and monitor events are logged on the station 
printer as well as reported through the Star Switch for even- 
tual recording on magnetic media. A typical printout of a Pio- 
neer 8 track is shown in Fig. 9. 


III. Related Applications 

This project produced several results that have application 
to other projects. A single board computer (available from 
Intel or National) was configured to function as a JPL stan- 
dard 15 line interface port. The module and software are appli- 
cable to any Multibus microprocessor and require no external 
circuitry other than a cable adapter from a ribbon to an MS 
type connector. This module is currently being installed in 
transmitter and maser controllers replacing the JPL designed 
quad 15 line interface requiring two Multibus slots. 

A general purpose parser was written to handle multiple 
line inputs of fields containing alphanumeric, integer, and 
fixed point values. The design accepts any nonalphabetic or 
nonnumeric character as a delimiter and considers consecutive 
delimiters as a single delimiter, thus providing free form in- 
puts. An ampersand character (&) is used to indicate the end 
of a physical line and the continuation of a logical line for long 
line inputs. 

Because of the requirement to update the local display with 
angles and analog data in near real time, a technique for fast 
conversion of floating point numbers to ASCII was developed. 
The technique is to separate the integer and fraction parts with 
the floating point processor into two 10 bit integers (the frac- 
tional part is scaled by 1 ,000). An assembly language routine 
then does two table look ups, using five bits at a time, to get 
4 bit binary coded decimal values. Use is made of the 8080 
BCD addition capability to combine the result. Conversion to 


66 



ASCII is then simply done by addition of a hexadecimal 30 to 
each BCD character. A floating point angle can be converted 
in about 300 microseconds. A faster technique, of course, is 
to simply use a table look up on all 10 bits. This was the origi- 
nal design before memory constraints required a different 
technique. 

Although not completed at this time, two features are being 
included to be used with the SETI project for performing spe- 
cial scanning requirements. One is a raster scan pattern speci- 
fied by overall size and scan rates. The second is a general 
three point fit capability to allow fitting a second order curve 
to any three points of right ascension-declination or azimuth- 
elevation. The current three day fit requires the points at zero 
hours on successive days. The more general method will allow 
the points to occur at unequal time intervals down to one 
second. 

A great deal of experience has been gained configuring a 
maximum loaded 8080 based system with multiple processors. 
Difficulty was encountered in debugging problems relating to 
critical timing relationships between the two CPUs. Present 
generation logic analyzers lack the capability of triggering on 
combined events in each CPU and back-tracing either CPU or 
the bus. 

Another inconvenience is the inability of the 8080 (and 
most other CPUs including 16 bit devices) to perform bit 
operations. The least addressable element is a byte. This means 
that to change a single bit in memory, an input-output port, or 
interrupt mask requires reading the byte, changing the desired 
bit, and storing the byte. If more than one CPU or an interrupt 
routine attempts to alter a bit in the same byte one of the 
changes may not be made. The Z80 CPU does allow bit opera- 
tions, which is desirable for control applications. 

All software was developed in the controllers or similar 
microprocessor computers. While this proved convenient from 
the standpoint of always having an editor or compiler avail- 
able, the time required for making changes in several modules, 
compiling, relinking, listing, and reloading typically was 
4 hours or more. Much of the debugging and patching was per- 
formed at the machine code level however, as one of the ad- 
vantages of PL/M is that the source statements are easily trace- 
able to the machine level code. 

Several undesirable features of the floating point processors 
were encountered and should be given consideration when 
choosing components for future systems. Stack oriented pro- 
cessors cannot be shared by multiple users (such as interrupt 
routines or other CPUs) if their state cannot be restored to the 
state at the time of the interrupt. The capability to save the 


processor stack and restore it existed, but the status register 
could not be restored. It would also be useful in many appli- 
cations to have access to the remainder following a fixed point 
division for modulo operations. Another annoying feature of 
the floating point processors was, if given operands that were 
out of range, an unpredictable result was generated rather than 
a closest approximation. In the refraction algorithm for 
example, it is necessary to evaluate small powers of e. If the 
argument is near zero, a very large result is obtained instead of 
unity. In time critical control applications, where it is not pos- 
sible to check the range of all operands before computation, 
a closest approximation would be desirable. 

The PL/M compiler available for the 8080 does not support 
floating point directly. Operands must be pushed onto the 
floating point processor stack, an operation code supplied, the 
result taken off the stack and stored. The execution time for 
getting numbers in and out of the stack is longer than many of 
the floating point operations. Some calculations can be 
chained to minimize this overhead but clever manipulation of 
a stack is one of the tasks that one hopes to avoid by using a 
high level language in the first place. The PL/M86 compiler 
for the 8086 CPU using an onboard floating point processor 
avoids this and allows floating point operations to be written 
algebraic form. 


IV. Current Status 

The system has been operational for several months for 
Pioneer 8 tracks and is now being used for all DSS 13 tracking 
operations. All features are not complete however, and work is 
continuing in the following areas: The general three point fit 
algorithm and raster scan for SETI are currently being imple- 
mented. The scheduler module to allow scheduled tasks to 
occur at a specified time and the capability of reading com- 
monly used source predicts from an on line file are not com- 
plete at this time. 


V. Conclusion 

While the objectives of this effort have been achieved, it is 
clear that this represents a data point on the maximum capa- 
bility of an eight bit microprocessor based system. Limits of 
memory size, execution speed, module space, and general com- 
plexity were encountered. All of these limits could be cor- 
rected by use of one of the more powerful 16 bit processors 
that have become available since the start of this project. The 
software is transferable to a 16 bit processor with minor 
modifications. 
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Fig. 1. DSS 13 antenna controller functional block diagram 



Fig. 2. Data collection controller configuration 
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Fig. 3. Main controller configuration 









Fig. 4. Data collection controller top level software 
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Fig. 6. CPU 1 level 2 software 
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Fig. 7. CPU 2 level 2 software 










HELP 

HELP MENUS - TYPE: 

HELP 0 - station level commands 

HELP 1 - cnf formats 

HELP 2 - str formats 

HELP 3 - maint mode commands 

HELP O 

ANTENNA CONTROLLER STATION LEVEL COMMANDS: 

IN IT : TURNS ON ELECTRONICS AND HYDRAULICS. 

CNF : INDICATES START OF SCHEDULE BLOCK DATA. 

OFF : STOPS ANTENNA, TURNS ELECTRONICS AND HYDRAULICS OFF. 
OPERATE: EXAMINE AND INITIATE ACTION SPECIFIED IN SCHEDULE BLOCK 
STOW : ANTENNA TO STOW POSITION? ELECTRONICS, HYDRAULICS OFF. 
STOP : SLOWS ANTENNA, APPLIES BRAKES. 

STR : REPORT STATUS OR DATA. 

CLEAR : RESET ANTENNA LOCKOUT AFTER A FAULT 
HELP : DISPLAY HELP MENU 

HELP 1 

ANTENNA CONTROLLER CNF MENU 
SID <name><ari9les> 

3DAY CnameXS sets of ansles> 

AZEL <anslesXnirap-up crod> 

DFOF <ans1es>, Edirect offsets] 

AEOF <ansles>, Caz/cos el 3 
CLDO, Cel ear- DFOF ] 

CLPO , Eel ears AEOF, RDOF3 

RDOF <ansles>. Era, dec offsets] 

RAEO Cansle rates>, Erate offsets! 

CLRO, Eel ear rate offsets] 

HELP 2 

ANTENNA CONTROLLER STR MENU 
STR CNF : CURRENT CNF BLOCK 
STR ANA : ANALOG DATA 
STR DIG : DIGITAL DATA 
STR ANG : ANGLE DATA 
STR ALL : ALL DATA 
STR TANG : TIME, ANGLES 
STR : SUBSYSTEM STATUS 

HELP 3 

MAINTENANCE COMMANDS 

10- SET BRAKE AZ 30-SET BRAKE EL 

11 - RELEASE BRAKE AZ 31 -RELEASE BRAKE EL 

12- HIGH SPEED AZ 32-HIGH SPEED EL 

13- LOW SPEED AZ 33-LOW SPEED EL 

14- ELECTRONICS ON 34-ELECTRG'NICS OFF 

15- PUMPS ON 35-PUMPS OFF 

7- I NCR RATE AZ 9-INCR RATE EL 

S-DECR RATE AZ O-DECR RATE EL 

19- IGNORE LIMITS AZ 39- IGNORE LIMITS EL 

20- SET LIMITS AZ 40-SET LIMITS EL 

21 - PR I NT CURRENT INCREMENTAL RATES 
23-PRINT RATE INCREMENT 


Fig. 8. Help menus 
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053 

15: 

40: 

58 

REM 


053 

15: 

42: 

28 

F5; INIT* 


053 

15: 

42: 

30 

INITIALIZING 


053 

15: 

42: 

42 

PUMP 75HP START RLY = ON 

053 

15: 

43: 13 

PUMP 125HP START 

RLY= ON 

053 

15: 

43: 

28 

F5;CNF; 


053 

15: 

43: 

53 

F553DAY PN8 246.659 -21.586 

053 

15: 

44: 

04 

INIT COMPLETE 


053 

15: 

44: 

11 

F5; 247. 584 -21.720 

053 

15: 

44: 

29 

F5? 248. 510 -21.849 

053 

15: 

44: 

36 

F5; OPER; 


053 

15: 

44: 

54 

MOVING TO POINT 


053 

15: 

46: 

19 

ON SOURCE 


053 

15: 

46: 

20 

OPER COMPLETE 


053 

15: 

49: 

03 

F5;STR? 


053 

15: 

49: 

03 

ANTENNA IS UNDER 

REMOTE CONTROL 

053 

15: 

49: 

03 

AZ =* 206.447 EL 

» 28.389 

053 

15: 

49: 

03 

AVG WIND * 9.29 

MPH 

053 

15: 

49: 

17 

F55STR; 


053 

15: 

49: 

17 

ANTENNA IS UNDER 

REMOTE CONTROL 

053 

15: 

49: 

17 

AZ * 206.507 EL 

= 28. 366 

053 

15: 

49: 

17 

AVG WIND = 9.29 

MPH 

053 

15: 

49: 

56 

F5; DFOF .087 -.01 

053 

16: 

09: 

02 

F5;str; 


053 

16: 

09: 

02 

ANTENNA IS UNDER 

REMOTE CONTROL 

053 

16:09: 

02 

AZ = 211.271 EL 

a 26.411 

053 

16: 

09:02 

AVG WIND * 10.7 

MPH 

053 

16: 

19:01 

F5; STR; 


053 

16: 

19: 

01 

ANTENNA IS UNDER 

REMOTE CONTROL 

053 

16: 19: 

01 

AZ = 213.548 EL 

a 25.322 

053 

16: 

19: 

01 

AVG WIND = 12.0 

MPH 

053 

16: 

24: 

22 

F5; STR; 


053 

16: 

24: 

22 

ANTENNA IS UNDER 

REMOTE CONTROL 


PS I 
PS I 
PS I 
PS I 
PS I 
PS I 
PS I 
PS I 
PS I 
MPH 
MPH 
DEG 
DEG 
GPM 


053 16:24:22 AZ = 214.738 EL « 24.712 


053 16:24:22 AVG WIND = 11.8 MPH 
053 16:24:52 F5;STR ANA; 

053 16:24:52 A32 EL LS MClTR DIF PRESS -17.6 
053 16:24:52 A33 EL HS MOTR DIF PRESS 13.2 
053 16:24:52 A34 EL SYSTEM PRESS 2880 
053 16:24:52 A35 EL HS PRESS 82.0 
053 16:24:52 A36 EL LS PRESS 2630 
053 16:24:52 A37 AZ HS MOTR DIF PRESS -11.7 
053 16:24:52 A38 AZ LS MOTR DIF PRESS 341. 
053 16:24:52 A39 AZ HS SYSTEM PRESS 2190 
053 16:24:52 A40 P>lt LS SYSTEM PRESS 1810 
053 16:24:52 A41 WIND SPEED SW TOWER 16.6 
053 16:24:52 A42 WIND SPEED SE TOWER 14.9 
053 16:24:52 A43 WIND AZ SW TOWER 349. 
053 16:24:52 A44 WIND AZ SE TOWER 302. 
053 16:24:52 A45 PUMP VOL 75HP LEFT 27.2 


Fig. 9. Typical operation record 
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NASA is participating in the development of international standards for space data 
systems. The purpose of this article is to provide a TDA assessment of the recommenda- 
tions for standards thus far developed. The Recommendations for Telemetry Coding and 
Packet Telemetry provide worthwhile benefit to the DSN; their cost impact to the DSN 
will be small. They will be of advantage to the NASA space exploration program. Their 
adoption should be supported by TDA, JPL, and OSTDS. 


I. Introduction 

NASA is participating in the development of international 
standards for space data systems. Work on some of the system 
elements is mature and is being readied for concurrence by the 
participating space agencies. 

The purpose of this article is to provide a TDA assessment 
of the recommendations for standards thus far developed. 
Recommendations for Telemetry Coding and Packet Teleme- 
try are considered in detail; others, which are less well devel- 
oped, are considered briefly. 

A. Summary of Findings 

The Recommendations for Telemetry Coding and Packet 
Telemetry provide worthwhile benefit to the DSN; their cost 
impact to the DSN will be small. They will be of advantage to 
the NASA space exploration program. Their adoption should 
be supported by TDA, JPL, and OSTDS. 

Six other Recommendations were also reviewed. They all 
are less mature than the Telemetry Recommendations. They 
deal with standardizable elements of the system, and it is 
expected that they will lead to adoptable standards after the 
appropriate additional work. 


B. Background 

The Consultative Committee for Space Data Systems 
(CCSDS) is in the process of developing a set of Recommenda- 
tions for standardization of space data system functions by 
space agencies. The purpose of the development is to enable 
economical cross-support between space agencies. From 
NASA’s standpoint, it is expected that the Standards as issued 
by NASA in accordance with the Recommendations will also 
facilitate support of NASA’s own missions by the Deep Space 
Network (and by TDRSS). 

The CCSDS has representatives of European, Japanese, and 
other national space agencies, including NASA. JPL has several 
active members on the CCSDS Panels. The process of develop- 
ing CCSDS standards Recommendations involves a sequence of 
increasingly mature documents: Concept Papers, White Books, 
Red Books, Blue Books. The details of the process are de- 
scribed in the CCSDS Recommendations Documents (e.g., 
CCSDS Recommendations for Space Data Systems Standards, 
Packet Telemetry, (Red Book) lssue-1, September 1983, 
page ii: “Document Status”). 

The Telemetry Coding and Packet Telemetry Recommenda- 
tions are the farthest along. They are being readied for formal 
concurrence in March, 1984. The Coding Recommendation is 
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based on mature technology. It is essentially the technique 
planned for use by Voyager 2 at Uranus. 

The Packet Recommendation will produce a “leading stan- 
dard.” That is, its past applications to space data systems are 
few and have only demonstrated the concept, but its future 
potential is considered good. For example, a relatively simple 
adaptive packet telemetry system was used for the low rate 
channels on Seasat A with good success (Ref. 1). The system 
was very effective in supporting the ground data handling and 
distribution functions. Also, the JPL ISPM spacecraft data sys- 
tem design was initiated, but not completed, with a packet 
telemetry structure closely related to that of the present 
Recommendations. The JPL Mars Geoscience/Climatology 
Observer Project has made a preliminary assessment of using 
packet telemetry and has concluded that it would be cost 
effective (private communication). 

In addition to the two Telemetry Recommendations, there 
are a number of others being prepared by the CCSDS Panels. 
They are: 

(1) Time Code Formats (Draft Red, Blue Books) 

(2) Packet Telecommand (White Book) 

(3) Radio Frequency and Modulation (Concept Paper) 

(4) Standard Format Data Units (White Book) 

(5) Electronic Communications (pre-Concept Paper) 

(6) Radio Metric and Orbit Determination (pre-Concept 
Paper) 

Recommendations (2), Packet Telecommand, and (3), 
Electronic Communications, are expected to have the greatest 
impact on DSN service. However, none of these six are further 
discussed here. 


II. Detailed Assessment of Impact of the 
Recommendations on the DSN 

A. Summary of Approach 

A block diagram of the DSN Telemetry System for the 
1988-90 period was developed (see Fig. 1). It is a direct evolu- 
tion of the mid-80’s MK IVA; it includes the currently planned 
SFOC/NOC concept. In this, the current split of functions be- 
tween the Mission Control and Computation Center (MC 3 ) 
and the DSN Network Operations Control Center (NOCC) is 
changed somewhat to reduce minor duplication, creating a 
Flight-Project implemented and operated Space Flight Opera- 
tions Center (SFOC) and a DSN Network Operations Center 


(NOC). Using the diagram, the elements of the system that 
were seen to be affected by the Telemetry Recommendations 
were identified by a study team from the TDA Office. 

Figure 1 depicts the DSN Telemetry System of the 1988-90 
SFOC/NOC era. It encompasses that portion of the end-to-end 
telemetry system for which the DSN has responsibility. The 
diagram identifies those changes to the DSN Telemetry System 
that will be required to support missions using the CCSDS 
Telemetry Coding and Packet Telemetry Recommendations. 
It shows telemetry frame sync and Reed-Solomon decoding 
being done at the Signal Processing Centers (vs. the NOC) in a 
new “DSN Standard Decoder.” The station location is not 
required by the Telemetry Recommendations, and may turn 
out not to be desired. In any case, the location chosen is inde- 
pendent of the requirements to support the standards. The 
Telemetry System and interface functions requiring changes 
related to the adoption of the Recommendations are identified 
in the boxes of Fig. 1. 

The marginal implementation and operating costs for the 
DSN to adopt the Telemetry Recommendations were esti- 
mated. Also, the operating characteristics of the DSN incorpo- 
rating the Telemetry Recommendations were assessed. 

Finally, the suitability of the Recommendations for the 
JPL end-to-end space data system (Ref. 2) was considered, 
from a TDA perspective, by study team members who have 
participated in the JPL effort to develop the CCSDS Recom- 
mendations. 


B. Assessment of Telemetry Coding 
Recommendation 

1. Description. The Recommendation is CCSDS Blue 
Book, Draft Issue-0, February 1984. The coding is for the link 
between the spacecraft and the ground signal processing cen- 
ter. Proper coding and decoding of the link improves its per- 
formance by a factor of 5 to 7 dB at a bit error rate (BER) of 
10~ s . The very low BER is needed to transport spacecraft 
science data that have been compressed and to enable effec- 
tive use of packet data streams (Ref. 3). 

The Recommendation proposes both an inner and an outer 
code. The inner code is a rate 1/2 constraint length 7 convolu- 
tional code with Viterbi decoding (Ref. 4). The outer code is 
a Reed-Solomon (255,223) block code with symbols from the 
256-element field interleaved to depth 5 (Ref. 5). The inner 
code can be concatenated with the outer code, or it can be 
used separately as in Ref. 4. The coding is compatible with, 
but does not require the use of, the proposed packet telemetry 
concept. Also, an uncoded link is allowed. 
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Two changes to the original Red-Book Recommendation 
proposed as a result of the TDA assessment have been 
accepted by the CCSDS Panel, and incorporated into the draft 
Blue Book, upon which the discussion of this article is based. 
They related to standardizing on soft decision Viterbi decod- 
ing so that deep space missions can be assured of cross- 
support, and putting the packet sync word in the Reed- 
Solomon code blocks so that coding and packets are indepen- 
dent to aid testing and to allow independent evolution. They 
also foreshadowed the possible need to go to constraint length 
7 rate 1/3 convolutional codes for some missions (Ref. 6). 

2. Impact on DSN. Minor software and hardware changes 
are required in the Viterbi decoder, the Telemetry Processing 
Assembly, and their monitor and control interfaces. Also, 
frame sync and Reed-Solomon decoding hardware and soft- 
ware are required in the DSN to support the Recommendation. 
However, those sync and decoding functions must be provided 
by the DSN to support the SFOC/NOC concept for committed 
mission support, independent of the adoption of the stan- 
dards. As we will see in a subsequent section of this article, 
the advantages of the Coding Recommendation outweigh the 
disadvantages. 

Table 1 gives a summary of estimated implementation and 
operations costs for adopting the Telemetry Coding Recom- 
mendation. The costs involved are small ($70K). Specifically, 
the delta cost of providing the standard decoding capability is 
limited to the cost of implementing the capability to reverse 
the order in which the data symbols are sampled in the maxi- 
mum likelihood-convolutional decoder. This single change pro- 
vides compatibility with the “blue book” standard for Viterbi 
decoding. 

Additional operational costs attributable to implementation 
of the Coding Recommendation are negligible, as testing and 
training of operators for a new decoder should be no more 
expensive than for maintaining the old capability. (Likewise, 
operational support of packet telemetry should be no more 
difficult to provide than with the present mission-unique inter- 
faces.) In fact, once a standard has been adopted, testing and 
training should become simpler, as support will be required for 
only one data interface for all projects. 

The Coding Recommendation is fundamentally sound. 
Some details will probably need modification in application. 
The basic technology embodied in it has had extensive devel- 
opment and flight/ground application. It is a well proven and 
effective concept for deep space missions. Concatenated 
convolutional/Reed-Solomon coding, in essence identical to 
that defined by the Recommendation, will be supported by 
the DSN and flight proven for deep space use by Voyager 2 at 
Uranus (Ref. 7). The recommended convolutional coding/ 


decoding has been used by Voyager 1 and 2 from launch and 
by TDRSS; it will be used by Galileo, ISPM (ESA), Giotto 
(ESA), and AMPTE (Ref. 8). Galileo and Giotto will use con- 
catenated Reed-Solomon coding very similar to that of the 
Recommendation (Ref. 9). 

Cross support of one agency’s spacecraft by another 
agency’s tracking network (Ref. 10) will be simplified by 
embracing the Recommendation. That will be a benefit to the 
DSN in the future. Ultimately, when all spacecraft not con- 
forming to the standards have expired, fewer code options will 
require sustaining and operational support by the DSN. That 
will reduce DSN support costs. 

Although technology developments during the 90’s may 
allow a 1 to 2 dB increase in link capability, the proposed 
standard is expected to have good longevity through the 90’s 
for deep space mission use. However, it may be that the stan- 
dard will be modified to permit constraint-length 7 rate 1/3 
convolutional codes. This would permit a gain of 0.5 dB 
(Ref. 6) at little additional complexity on spacecraft and 
ground. Thus, the Coding Recommendation was endorsed to 
OSTDS as desirable for the DSN and NASA. 

C. Assessment of Packet Telemetry 
Recommendation 

1. Description. The Recommendation assessed was CCSDS 
Blue Book, Draft, Issue-0, February 1984. The packets, con- 
sisting of data and an identifying header, are made up by the 
spacecraft individual data sources. They are asynchronously 
time multiplexed into the spacecraft data stream. They pass 
through the system intact to the user of the data. For transfer 
from the spacecraft through delivery to the Project’s mission 
control, the packets are carried in “telemetry transfer frames.” 
The frames carry information on the mission and instructions 
for extracting the packets. The frames are multiplexed syn- 
chronously into the data stream. On the ground, synchronizing 
to and reading of the frame header allow the data stream to be 
directed toward the proper spacecraft Project. Then the indi- 
vidual packets can be demultiplexed and read by the Project 
(Ref. 11). 

The concept enables multimission identification, extrac- 
tion, and distribution of the data in the stream. It also enables 
but does not require an adaptive data stream that reacts to 
events on board the spacecraft. That capability is expected to 
benefit spacecraft projects and simplify ground processing. 
The telemetry link efficiency with a data stream compatible 
with the Packet Recommendation should be comparable to 
that of current deep space data system designs. Indeed, calcu- 
lations show that the approximately 0.25 dB average link inef- 
ficiency which can be computed (details omitted) for pack- 
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etized telemetry is about the same as or lower than the ineffi- 
ciencies computed by the JPL Telecommunications Systems 
Section (private communication) for Voyager and other fixed- 
format time-division multiplexed telemetry. 

The Recommendation as we have said defines the format 
for the packets and the transfet frarrte. The detail in the 

1 ‘ k 

Recommendation document is sufficient to insure compatibil- 
ity between users. Several options are available within the 
Recommendation to handle packets that are longer than a 
single frame permits. The Recommendation is compatible 
with, but does not require, per se, the use of the channel cod- 
ing recommended. It does, however, require low error proba- 
bility in the channel, as can be provided by proper application 
of the coding recommendation. 

2. Impact on DSN. In concept, the packets, carried by 
their transfer frames, travel through the ground data system 
unobserved until they are demultiplexed from their frames at 
their destination. Therefore, the impact on the DSN is small. 

Currently (Ref. 12), the DSN NOCC finds the telemetry 
frame header and does some demultiplexing to provide data re- 
quired by DSN performance analysts. The simple algorithm 
used will not suffice to demultiplex an asynchronous packet 
stream. However, in the SFOC/NOC era, the DSN NOC is to 
receive the required information from the SFOC’s demulti- 
plexing process. This appears to be a desirable feature of the 
SFOC/NOC concept, particularly for support of asynchronous 
packetized data streams. 

Full and effective use of the packet data concept requires 
very low error rates in the end-to-end system. The error prob- 
ability requirements are two to three orders of magnitude 
more stringent than generally needed with previous deep space 
imaging data systems. It is not evident that any system changes 
to the DSN Complexes or the GCF will be needed solely to 
support the more stringent requirements. 

Surely, however, more rigorous performance validation will 
be required. The costs have been estimated for providing auto- 
matic retransmission or forward error correction on the GCF 
present or future 56- to 600-kbps wideband data lines. This 
capability will be appropriate for real-time transmission of 
high-rate compressed images or packet telemetry. However, 
the costs (S360K) were not included as a marginal cost for 
adopting packet telemetry. This is because the capability 
would only be provided in response to a specific mission 
requirement and could be driven by missions not using 
packets, such as missions using data compression (Ref. 3). 
This low-error capability could also be required by a mission 
using packets but not in conformance with the CCSDS 
Recommendations. 


The estimated DSN implementation and operational costs 
to adopt the Packet Telemetry Recommendation are also given 
in Table 1. The cost impact, without the GCF implementation 
noted above, is small ($60K). Thus, the Packet Telemetry 
Recommendation was endorsed to OSTDS as desirable for the 
DSN and NASA. 

D. Considerations for Missions Not Using Telemetry 
Standards 

There are no spacecraft currently designed or in flight that 
are fully compatible with the Recommendations. For their 
lifetimes, those spacecraft must be supported by DSN operat- 
ing modes, or equipment, different from that for standardized 
spacecraft. Also, it should be anticipated that some future 
spacecraft designs, for reasons of overall economy, will use 
“inherited” flight equipment that will result in data systems 
incompatible with the standards. 

For those reasons, for a long time into the future, capabil- 
ity for support of nonstandard spacecraft data systems will 
have to be retained by the DSN. Even so, adopting the stan- 
dards is expected to lead to fewer, not more, telemetry designs 
to support during the period. The basic approach to planning 
for removal of no longer needed capability in the DSN will be 
to obtain concurrence of affected projects, both present and 
potential future, and of affected NASA Centers, before re- 
moval of existing capability. 

E. Advantages and Disadvantages 

Table 2 lists various advantages and disadvantages of the 
Telemetry Standards that were uncovered in the course of 
the TDA study. They are listed in the table without quanti- 
fication. Nevertheless, it is seen that advantages far outweigh 
disadvantages. 


III. TDA Perception of Effect of Standards for 
Space Data Systems 

A. Introduction 

This section presents the TDA assessment of how Data 
System Standards based on the CCSDS Recommendations 
will affect the way JPL and NASA conduct missions in the 
future. Issues of cross-support, mission design, mission oper- 
ations, and space program growth are assessed in this light. 
The overall conclusion is that the Standards will promote JPL 
and NASA programs by reducing costs and increasing coopera- 
tive ventures. 
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B. Cross-Support 

The greater interoperability that the Standards enable will 
make it easier for NASA to give and get cross-support. A 
greater number of cooperative missions between space agen- 
cies, such as the International Solar Terrestrial Physics Pro- 
gram (Ref. 13), can be expected. The common performance 
standard must be maintained high enough so that NASA can 
receive cross-support as well as give it. This is because NASA’s 
deep space missions usually require performance near attain- 
able limits. 

C. Mission Design 

The Standards can benefit both the U.S. and International 
space programs in the mission design phase. It is expected that 
the Packet Telemetry Standard will facilitate mission- 
independent instrument integration including re-use and re- 
flights. Adaptive telemetry and data compression become 
easier. 

D. Mission Operations 

The benefits can be substantial, as the SEASAT A opera- 
tions experience has shown. It will become much easier to 
handle data from acquisition to distribution and archiving. 
Common program data, e.g., planetary data system, will be 
easier to provide. Likewise, coordination of data from multi- 
ple spacecraft, such as the ISTP Program and Space Station, 
will be easier. It will become much easier to have distributed 
operations and science teams. Common hardware and software 
are expected to reduce costs in the mission operations system. 
Overall, the Standards will benefit mission operations. 

E. Space Program Growth Potential 

The international Space Data System Standards can be a 
strong factor in the growth of the world’s space program. The 
Standards should cut lead time and cost of experiment integra- 
tion. Multinational experiments on a space agency’s mission 
will be easier, as well as multispacecraft multiagency missions. 
Cross-support will be cheaper to obtain or provide, and this in 


turn can encourage participation by other nations in NASA 
programs. The availability and utility of space-derived data 
will increase. 


IV. Future Work to Support Standards 

For the U.S. to be a full participant in and beneficiary of 
the international Space Data System Standards program re- 
quires a continuing effort. This effort is to evolve the Stan- 
dards to take advantage of new spacecraft and ground technol- 
ogy and mission concepts. The possibilities for enhancing 
space programs require a broad planning perspective for 
NASA, JPL, and TDA within JPL. The TDA responsibility is 
to determine performance and interface specifications and 
verify system performance against theory and specification. 
One change to be encouraged early is to promote the inter- 
national standardization of the common JPL/GSFC constraint- 
length 7 rate 1/3 convolutional code as an alternative to the 
constraint-length 7 rate 1/2 code of the current Guidelines. 
This can benefit missions which are not bandwidth constrained 
and can reduce costs or increase science value from the 0.5 dB 
increase in performance (Ref. 6). 


V. Conclusions 

The adoption of appropriate standards for space data sys- 
tems will be of long term benefit to NASA. It will reduce the 
cost of mission design, implementation, operations, and data 
analysis. These benefits will be obtained by more thorough 
development of spacecraft and ground mission and facility 
designs and operating procedures. Joint ventures between 
NASA and the space agencies of other nations will be made 
easier, better, and less costly. Space data will be more available. 

Two of the CCSDS Recommendations, Telemetry Coding 
and Packet Telemetry, have been concurred in by TDA at JPL. 
They have been well developed and have sound bases. They do 
not cause noticeable expense to NASA in the near term and 
can only save in the long term. Concurrence in them was there- 
fore recommended to OSTDS by the TDA Office at JPL. 
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T able 1 . Estimated DSN marginal costs 8 for implementing Telemetry 
Channel Coding and Packet Telemetry Recommendations 


Standards 
Cost Elements 

Channel Coding 

Packet 

Telemetry 0 

Viterbi 

Decoding 

R-S 

Decoding 13 

Hardware implementation 

47.5 

0 

0 d 

Software implementation 

23.9 

0 

57.5 d 

Operations 

0 

0 

0 


a All costs in thousands of FY 84 dollars, ±30% estimate. 
b Assumes use of current R-S Decoder hardware/software now 
available in SFOC. If decoding function were shifted to DSN 
SPC’s in future, there would be additional hardware and training 
costs. However, those costs are independent of adoption of 
standards. 

c Assumes fixed length transfer frames as per Recommendation. 

1:1 Auto recall of GCF WBDL (wideband data lines): cost estimated 
to $360K - $193K for hardware, $167K for software. Needed 
only for missions requiring real-time high data rate with low 
deletions. This would be a mission requirement, if it arises, inde- 
pendent of whether the CCSDS Packet Telemetry Recommenda- 
tions were adopted. 



Table 2. Characteristics and effects of telemetry guidelines 


Guideline 


Characteristic 


Effect 


Coding 


Positive 


Common channel coding for all DSN-supported missions 


Cross-support of U.S. missions by other Agencies and conversely 
Enables full use of Packet Telemetry 

Coding proposed is essentially that evolving as a de facto DSN 
standard for deep space mission support 


Allows mission-independent decoder software/hardware 


Reed-Solomon/Viterbi provides low error probability for all 
missions using it without bit power penalty 

Enables mission-independent frame sync 


Simplifies all phases (planning through design through 
operation) of DSN support of all missions (Deep 
Space, HEO, NASA, Non-NASA) 

Wider base, greater depth of understanding by users of 
coded channel capability and characteristics 

More cross-support potential 

More cross-support potential 

Builds on considerable DSN investment 

Technology sufficiently mature to support providing 
a stable standard 

Appropriate time frame for initial establishment of 
standard 

Reduced hardware/software at Complexes 

Reduced effort to test and verify performance for 
mission set 

Reduced effort for sustaining in Net 

Low-error-probability engineering telemetry may 
make isolation of spacecraft vs ground anomalies 
easier 

Reduced implementation and test costs and increased 
confidence 


Negative 


Retention of old nonconforming designs in overlap period Increased sustaining engineering costs 


Packet Enhanced ability to obtain cross-support 

Enhanced ability to provide cross-support 
Makes delivery of data to Projects easier 

Eliminates potential need for multiple subcarriers 


Positive 


Reduced NASA costs if facilities of other agencies 
are needed 

Broader base of DSN customers; more joint missions 

Reduced end-to-end system test and operations 
costs, increased user satisfaction, and greater cross- 
support potential 

Fewer Subcarrier Demodular Assemblies and Symbol 
Synchronizer Assemblies needed at Signal Processing 
Centers 
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Table 2. (cont’d) 


Guideline 

Characteristic 

Effect 


Negative 


Packet 

Required packet system performance not specified 
Performance of full adaptive options not specified or understood 

Possible complication of project end-to-end testing involving DSN 

Need to retain non-packetized telemetry for old designs 

Potential user surprise and dissatisfaction 

Potential need to implement still better coding or 
provide more signal-to-noise some other way 

May increase Signal Processing Center time needed for 
testing 

Possible enlarged complexity in Telemetry Processing 
Assembly and increased sustaining 
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Fig, 1. Baseline 1988 - 90 era DSN Telemetry System 
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The Goldstone R/D High Speed Data Acquisition System 

L. J. Deutsch, R. F. Jurgens, and S. S. Brokl 
Communications Systems Research Section 


A digital data acquisition system that meets the requirements of several users (initially 
the planetary radar program) is planned for general use at Deep Space Station 14 
(DSS 14). The system, now partially complete, is controlled by a Digital Equipment 
Corporation (DEC) VAX 11/780 computer that is programmed in high level languages. 
A DEC Data Controller is included for moderate-speed data acquisition, low speed data 
display, and for a digital interface to special user-provided devices. The high-speed data 
acquisition is performed in devices that are being designed and built at JPL. Analog IF 
signals are converted to a digitized 50 MHz real signal. This signal is filtered and mixed 
digitally to baseband after which its phase code (a PN sequence in the case of planetary 
radar) is removed. It may then be accumulated (or averaged) and fed into the VAX 
through an FPS 5210 array processor. Further data processing before entering the VAX is 
thus possible (computation and accumulation of the power spectra, for example). The 
system will be located in the research and development pedestal at DSS 14 so easy 
access by researchers in radio astronomy as well as telemetry processing and antenna 
arraying is anticipated. 


1. Overview 

A High Speed Data Acquisition System is being constructed 
for research and development (R/D) use in the Deep Space 
Network (DSN). This system will, when completed, be a sub- 
system of the DSN Ground-Based Planetary Radar System. Its 
primary purpose will be to collect the 15 MHz bandwidth RF 
signal from the DSN Radar Receivers and perform all neces- 
sary data rate compression (filtering, phase correction, and 
correlating) to get radar data on a general purpose computer 
for further processing. This function is currently performed by 
an array of hardware devices that are, for the most part, inade- 
quate for the type of high bandwidth radar experiments that 
will be performed over the next decade. A good idea of the 
capabilities of the current system can be seen in Refs. 1 and 

2. This article provides a system level understanding of the 
new system. 


The High Speed Data Acquisition System is a collection of 
advanced digital signal processing devices that are connected to 
a host minicomputer. The system will be configured to operate 
in a wide range of modes, allowing a user to interconnect these 
devices and control them completely under software control. 
The set of processing devices that will be available to users 
enables many types of radar experiments as well as pulsar 
observations, general radio astronomy, and some telemetry 
demodulation and decoding capabilities. 

A Digital Equipment Corporation (DEC) VAX 1 1/780 com- 
puter is the heart of the system. The standard DEC operating 
system (VMS) will be used to run the machine while all of the 
general purpose peripheral devices will be connected to the 
VAX through conventional VMS driver routines. Conse- 
quently, nearly all data acquisition programs (i.e., software for 
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the control of experiments) may be written in high level 
languages such as FORTRAN or Pascal. This design philosophy 
provides a means for low cost data collection and ease of docu- 
mentation. It will also enable more people to operate the sys- 
tem with only a small amount of specialized training. 

The configuration of the High Speed Data Acquisition Sys- 
tem at completion is shown in Fig. 1. The top part of the fig- 
ure represents the VAX computer and its peripherals while the 
bottom part shows the special purpose devices that will be 
built. 

The VAX 11/780 computer is also configured with all the 
necessary peripherals for program development and testing. In 
addition there are commercial peripherals that are essential to 
the radar program, including a high resolution color video 
monitor for real time image display and a high resolution 
graphics printer for data monitoring during experiments. Data 
will be collected and saved on a 500 Mbyte Winchester disk or 
on two high speed magnetic tape drives. The storage capability 
on this system will be nearly a factor of ten better than that 
of the current radar system. 

In order to process the high rate data imposed by antici- 
pated signal bandwidths, special purpose devices are also re- 
quired. These high speed digital signal processing devices will 
be connected to the VAX through a Floating Point Systems 
(FPS) array processor. Subsequently, the data may be further 
reduced and processed before being transferred to the VAX. 
This will relieve the typical throughput problems that one 
encounters when using a host computer to collect and process 
real time data. 

In addition, there will be low rate data input and output 
available on the VAX by means of a set of standard DEC 
analog-to-digital and digital-to-analog converters. Data col- 
lected in this manner by the VAX can be transferred to the 
array processor for processing and then back to the VAX for 
handling or storage. The digital-to-analog converters can be 
used to drive instruments or real time displays in order to 
enhance certain experiments. 

The radar correlator assembly consists of a block of sixteen 
correlator modules that will run at a clock rate of 10 MHz. 
Each module will be capable of performing 256 delay lags in a 
one by four bit correlation function in a continuous stream. 
Each individual module will be able to operate in a cross- or 
auto-correlation mode. The selection of input signals, refer- 
ence signals, and interconnection of the modules will be under 
complete control by the host VAX. This means that a user can 
reconfigure the correlator assembly to suit the needs of each 
individual experiment, or change their function during an 


experiment in a controlled manner, all from software written 
in a high level language. 

The radar system typically uses up to three complex input 
channels. However, a fourth complex channel is included in 
the new system. This extra channel provides for expansion as 
well as a hardware backup. 

Since some radar experiments will require pulsed or spread 
spectrum radar as well as CW radar, pseudonoise (PN) se- 
quence generators must also be part of the High Speed Data 
Acquisition System. One generator is required for each chan- 
nel. The output of these generators will be supplied to the cor- 
relator assembly, the IF processing devices, and the transmitter 
driver. These PN generators are programmable (from the VAX) 
so that a large set of PN codes of various lengths can be used 
for different experiments. 

A unique feature of the new system is the IF processor 
hardware. Input signals are sampled directly at the output of 
the IF amplifier which is centered at 7.5 MHz. With this center 
frequency, a channel bandwidth of 15 MHz is available. Thus, 
a sample rate of at least 30 MHz will be needed to handle the 
information band. A sampling rate of about 50 MHz will be 
used to digitize the IF signal as the input to the system. All 
signal conditioning and processing following this will thus be 
accomplished digitally. The conditioning processes will consist 
of antidispersion filters, baseband mixers, and baud inte- 
grators. The resulting baseband signal will then be fed to a set 
of signal multiplexers that will be used to distribute the vari- 
ous signals to the correlator assembly. 

A wide range of preprocessing can be achieved, ranges from 
none to the full capability of filtering, correlating, and trans- 
forming. For this reason, it is better to view the High Speed 
Data Acquisition System as a four channel digital receiver sys- 
tem. Among the nonradar applications that this system will be 
capable of supporting are general radio astronomy, SETI 
(Search for Extra-Terrestrial Intelligence), and telemetry 
reception. These applications are addressed in more detail in 
the following section. 


II. Applications 

The primary mode of operation of the High Speed Data 
Acquisition system is the collection and processing of radar 
data. Since the system was designed with this purpose in mind, 
it will do this job particularly well. The system will support 
CW, pulsed, and spread spectrum radar experiments, and thus 
it will support experiments targeted at the outer planets as 
well as the terrestrial planets. A number of radar experiments 
anticipated over the next five years are listed below along with 
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their system requirements. Figure 2 shows a time table of the 
various experiments. 

Mars. This will be the first experiment that will use the 
High Speed Data Acquisition System. The observations are 
planned for June and July of 1984. Only one station (DSS 14) 
will be used. This experiment does not require the IF proces- 
sor portion of the system. It will require at least one quarter of 
the correlator assembly, a PN sequence generator, and all of 
the computing facility. In addition, a baseband mixer and 
sampler assembly must be fabricated (mostly out of existing 
parts) to interface the system to the IF output of the radar 
receiver. 

Mercury. Mercury can be observed several times each year. 
This work will require either a dual polarization receiver or the 
use of two stations operating at different polarizations. In 
either case, at least two channels of the IF processor must be 
in place. Mercury will be observed with both CW and spread 
spectrum radar. 

Venus. Venus will be observed in the first six months of 
1985. These observations require three stations. The entire IF 
processor must be working by this time. Also, the full capa- 
bility of the correlator assembly will be needed. 

Asteroids. These observations will be made on a semiregular 
basis. CW radar with frequency hopping and a dual-polarized 
receiver will be used. 

Comets. These observations are similar to those of asteroids. 
Of particular note are the proposed observations of Halley in 
early 1986. 

Moons of Jupiter. These observations will resume in April 
of 1985. They will also use CW radar with frequency hopping 
and dual-polarization. 

Rings of Saturn. The first observation will be in May of 
1985. CW radar with frequency hopping and dual-polarized 
signals will be used. 

The High Speed Data Acquisition System will also be a use- 
ful tool in general Radio astronomy. It could be used to mea- 
sure Stokes parameters (polarization parameters) and average 
the results in real time. It will also be possible to run the sys- 
tem as a real time interferometer by using the multiple chan- 
nels in the IF processor. 

The flexibility that is being designed into the High Speed 
Data Acquisition System will make it a useful tool in DSN 
receiver development. It will be the first system at a DSN site 
to have a fully digital IF processor. Since all aspects of this IF 


processor are to be programmable from the VAX computer, it 
can be used to experiment with new types of mixing and filter- 
ing techniques for telemetry and radiometric reception. It will 
be possible to control the loop filters dynamically during 
reception to provide improved receiver performance. The four 
channels in the IF processor can be used to perform real-time 
arraying of up to four antennas as well as carrier arraying for 
telemetry reception. The correlator assembly could also be 
used to demodulate a subcarrier and track data in a telemetry 
signal. Telemetry decoding could then be done on the VAX if 
the data rates are low enough. If the system were to be used 
for extensive high rate telemetry experiments, then a Viterbi 
decoder (MCD) could be attached to the VAX as another high 
speed peripheral. 

The extreme flexibility in the High Speed Data Acquisition 
System and the ability to control this flexibility from simple 
high level programming languages is certain to make the sys- 
tem a major research tool for many projects well into the 
future. 

III. Computer Facility 

The High Speed Data Acquisition System is based on a 
Digital Equipment Corporation VAX 11/780 computer. This 
computer was chosen to meet the need for high-performance 
real time data collection, data processing, and program devel- 
opment. This machine, and all its associated peripherals and 
hardware add-ons, will be located at the Goldstone facility in 
the pedestal of DSS- 14. 

The VAX computer uses a 32-bit architecture that is based 
on the PDP-11 family of 16-bit minicomputers. Although the 
VAX uses addressing modes and stack structures that are 
similar to those on PDP machines, the VAX provides 32-bit 
addressing. This is enough to access a four gigabyte program 
address space. The 32-bit data paths in the VAX also offer 
increased processing speed and accuracy over the PDPs. The 
VAX processor hardware and instruction set were specifically 
designed to support high-level programming and still produce 
efficient machine code. The VAX also has a floating point pro- 
cessor that augments the primary processor by performing 
floating point operations with a substantial performance 
improvement. There is also a high speed memory cache that 
provides a large reduction in memory access time. The VAX 
main memory uses error correcting codes to correct single bit 
errors and detect double bit errors in each byte. 

The VAX operating system, VAX/VMS (or simply VMS) 
is a multi-user virtual memory system that supports an easy-to- 
use interactive command interface. The virtual memory 
feature allows the execution of programs that can be larger 
than the physical memory allocated to them. Virtual memory 
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paging (the operation that swaps memory and disk space to 
allow this feature) is handled automatically by VMS. Conse- 
quently, the programmer need not be concerned with the 
working details of memory management. Memory manage- 
ment facilities can be controlled by the user in the event that 
such control is necessary to optimize the real-time perfor- 
mance of a particular program. 

There are two hardware buses on the VAX backplane to 
which peripherals may be attached. One of these is the MASS 
Bus. It is capable of high throughput communications. The 
second, the UNIBUS, provides a standard interface that can be 
used with a large set of peripherals. Figure 3 shows the con- 
figuration of the High Speed Data Acquisition VAX. The 
blocks that have shaded corners represent equipment that is 
already installed and working in a JPL Laboratory. The 
peripherals that will be attached to the VAX are listed below 
(but do not include the special purpose hardware devices that 
will be described in the following sections): 

2 RK07 28 Megabyte Disk Drives. These drives use remov- 
able hard disk cartridges. These are useful for storing the oper- 
ating system and most user software. Since they are remov- 
able, they provide a fast method for transporting software or 
performing backups. 

1 RP07 500 Megabyte Disk Drive. This device, which is to 
be attached to the MASS bus, will be a principal storage area 
for collected data. 

2 TU78 Magnetic Tape Drives. These are high speed (1600 
or 6250 bpi at 125 ips) tape drives that can be used in real 
time data acquisition or high speed data transfers. 

1 DR 780 High Speed 32-bit Parallel Interface. This device is 
a 32-bit parallel interface that is capable of transferring data at 
rates up to two megabytes per second. The DR780 uses sepa- 
rate data paths for control information and data. The control 
path is 8-bit and bidirectional. The data path provides a 
synchronous 32-bit transfer which may be clocked by either 
an LPA11-K clock or by a clock provided by an external 
device. 

1 LPA1I-K Analog/ Digital Data Controller. This device is 
an intelligent controller that transfers data between the VAX 
computer and a set of interface modules. These modules 
include analog-to-digital and digital-to-analog converters and a 
high speed 16-bit parallel interface (see Section IV). The 
LPA11-K subsystem uses silo buffers and direct-memory- 
access data transfers to allow real-time data acquisition with a 
very limited interrupt load and very limited burden in the 
VAX itself. The LPA1 1-K is implemented using two micro- 
processors, and it contains its own programmable real-time 


clock. The VAX can initiate a chain of data transfers through 
the LPA1 1-K with only a single I/O request. This silo structure 
allows data to be transferred continuously regardless of the 
status of the UNIBUS or system software latencies. Data rates 
up to 15,000 samples per second can be handled by the 
LPA11-K. There is also a mode where the LPA11-K services 
only the analog-to-digital converter and can handle 150,000 
samples per second. 

6 DR11-C 16-bit Parallel Interfaces. These will be used to 
provide control to the various special purpose hardware that 
will be described in Section IV. 

2 VT100 CRT Terminals. These will be the primary user 
interface with the VAX for both program development and 
experiment control. 

1 AED512 Color Graphics Terminal. This will be used as a 
real time high resolution display for monitoring experiments 
and reviewing collected data. The AED512 provides a 512 X 
512 pixel display with 8-bit resolution in each of the three 
primary colors for a total for 24 bits for each pixel. 

1 LA 120 Operator Console. This device will be used to 
monitor the VAX itself and keep hardcopy logs of system per- 
formance. Hardware problems will be reported to the console 
so that the operator may act upon them when necessary. 

1 LXY12 Printer/Plotter. This device will be the primary 
hardcopy output generator. It will be used to examine quick- 
look data analysis and to aid in software development. 

1 FPS-5210 Array Processor. This device is equipped with a 
programmable I/O processor (GPIOP) that can be used to 
collect data from an external device and pass it through the 
array processor independently of the VAX. Maximum data 
rate through the GPIOP is three million 38-bit words per 
second. The FPS 5210 is rated at 30 million floating point 
operations per second. It will be used to format and process 
incoming data before it is sent to the VAX. 

1 DEC PDS Power Distribution System. The PDS will pro- 
vide filtered and regulated power to all the VAX peripherals. 

The VAX provides the High Speed Data Acquisition with 
the required flexibility and expandability that are needed to 
support experimental work at Goldstone well into the future. 

IV. General Purpose Acquisition Hardware 

The High Speed Data Acquisition System will have a set of 
peripheral devices that perform conventional data acquisition 
functions. These functions include digital-to-analog eonver- 
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sion, analog-to-digital conversion, and high throughput digital- 
to-digital interfaces. These can be used to interface the system 
with existing data acquisition hardware, laboratory instru- 
ments, or specialized display devices. They will also be instru- 
mental in the testing of the system itself. 

The general purpose acquisition hardware will reside on its 
own dedicated bus so that a continuous data flow may be 
achieved regardless of the status of the buses on the VAX 
itself. This bus will be attached to the VAX through the 
LPA11-K Analog/Digital Data Controller that was described 
in Section III. Figure 4 shows the configuration of this part of 
the system. The shaded boxes, once again, represent hardware 
that is already installed in the laboratory. 

There will be four modules attached to this LPA11-K bus 
that will comprise the general purpose acquisition hardware. 
These are described below: 

AA11-K 4-Channel D/A and Display Control. This device 
provides four independently buffered digital-to-analog con- 
verters. Each of these has 12-bit accuracy and a slew rate of 
5V/s. In addition the AA1 1-K includes a display control driver 
that supports a 4096 X 4096 dot array display. This can be 
used to drive a high-resolution graphics CRT or an oscilloscope. 

ADI 1-K Analog to Digital Converter. The AA1 1-K contains 
a high speed 12-bit analog-to-digital converter and an associ- 
ated multiplexer. The multiplexer can be used to support 8 
full channels or 16 single-ended channels of input. The timing 
for the conversions can be provided in one of three ways: 
under software control, from a real-time clock, or from an 
external input. A multichannel throughput rate of 50 KHz will 
be possible. 

DR 1 1-K Interface. This is general purpose 16-bit parallel 
interface for interconnection with other digital devices. 

KW11-K Dual Programmable Real Time Clock. This module 
provides two real-time clocks that will be used to control data 
flow in and out of the LPA11 subsystem and provide accurate 
time tagging. Each clock is independently programmable and 
may operate at several different clock frequencies, as may be 
required for various applications. 

These devices, together with the LPA11-K controller, will 
constitute a subsystem of the High Speed Data Acquisition 
System that will provide a general compatibility with other 
data acquisition and data processing hardware without disturb- 
ing the other portions of the system. The inclusion of this sub- 
system will make it easier to append hardware devices to the 
system when such devices are needed without having to install 
these permanently. 


V. High Speed Data Acquisition Hardware 

In addition to the computing facility and the general pur- 
pose data acquisition hardware, the High Speed Data Acquisi- 
tion System will include special purpose hardware that will 
enable it to collect and process data directly from the IF por- 
tion of the radar receiver. This same hardware could be used to 
interface with other receivers as well. Figure 5 shows the 
special purpose hardware and its relationship to the computing 
facility. Most of this special purpose hardware is not available 
commercially and hence it must be designed and built in- 
house. The hardware with shaded boxes in the figure has been 
fabricated and tested. 

The high speed data acquisition hardware comprises three 
subsystems. These are the IF processor subsystem, the corre- 
lator assembly, and the PDPG coders (PN sequence genera- 
tors). Figure 5 shows these subsystems and their relationship 
to the High Speed Data Acquisition System computing 
facility. 

The IF processor will serve as the principal path for input 
data in the completed system. It will be designed to accept, as 
input, an analog IF signal with a carrier of 7.5 MHz and a two 
sided bandwidth of 15 MHz. Thus the analog-to-digital con- 
verters (A/Ds) must sample at no less than 30 MHz and will do 
so with four-bit accuracy. However, in order to assure that 
there is a sufficient rejection of frequencies outside of this 
bandwidth, the converters will actually be capable of a 
50 MHz sampling rate. 

The outputs from the converters will be passed through 
Anti-Dispersion Filters (ADFs) before mixing so that phase 
corrections can be made on a controlled basis. These filters 
will be programmable from the VAX. 

The Digital Complex Mixers (DCMs) will be high speed 
(50 MHz) digital multipliers, which will heterodyne the IF sig- 
nal down to baseband. The reference frequency for the conver- 
sion can come from the VAX (predictions stored in memory 
and sent to a PLO) or from feedback within the system. 

The signal will then be passed through Complex Baud 
Filters (CBFs). These devices will serve two purposes. They 
will act as loop filters that can be controlled by the VAX. In 
addition, they will smooth the signal so that it may be sampled 
at a lower rate for subsequent processing. 

The correlator assembly consists of a baseband signal multi- 
plexer, correlator-accumulator modules, and a device for data 
formatting and interfacing to the array processor. 

The baseband multiplexer will allow the baseband signals 
from the four complex channels of the IF processor to be 
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routed to any combination of correlator-accumulator modules. 
This function will be under the control of the VAX. In this 
way, the correlator assembly can be configured many ways, 
depending on the needs of a particular experiment. The corre- 
lators can have a variable number of lags (selectable under 
computer control) and operate as cross- or auto-correlators. 

Each correlator-accumulator module will include a com- 
plete 4-by-l 256 lag correlator and a pipelined adder that will 
produce a 12-bit output. These will operate at a clock rate of 
10 MHz. The accumulators will be capable of summing the 
correlator outputs for up to one full second. The output of 
each module will be a 28-bit word. 

The Interface and Data Formatting will collect the output 
from the correlator-accumulator modules and feed that data 
through the General-Purpose Input/Output Port (GPIOP) of 
the array processor. This collecting function will also be under 
the direct control of the VAX. 

The PN sequences that will be used for spread-spectrum 
radar work will be supplied by the polynomial driven PN Gen- 
erators (PDPGs). Each of these will run at a 10 MHz chip rate 
that will be synchronized with an accurate reference clock. 
The PDPGs will also provide time domain compensation for 
the Doppler effects caused by the relative motion of the 
observer and the target. The PDPGs will be controlled by the 
VAX and will be capable of generating a large range of PN 
sequences from a single pulse to a sequence 2 24 - 1 bits long. 

This special purpose hardware will compose the first high- 
speed digital IF receiver to be installed at Goldstone. It will be 
capable of serving as a test bed for digital receiver work as well 
as radar and radio astronomy work. 

VI. Interfacing to the Remainder of the Radar 
System 

In order to perform the experiments that were described in 
Section II, the High Speed Data Acquisition System must 
interface with other DSN systems, including transmitters, 
receivers, communication links, and frequency reference gen- 
erators. These systems are also undergoing fabrication or up- 
grading. It is essential that they are operating at the time that 
they are needed for the various experiments. We describe here 


the major requirements that the High Speed Data Acquisition 
System will put on these other systems: 

(a) Radar experiments must be performed with transmit- 
ters and receivers that operate on the same frequencies. 
Both S- and X-band signals will be required. 

(b) The bandwidth of the transmitter, receiver, and 
antenna systems must be adequate to pass a 2 MHz PN 
phase code without a serious degradation of shape for 
any length 2 N - 1 where 1 <N <30. 

(c) The X-band bandwidth must also support the Rings of 
Saturn experiment. At least 10 MHz must be passed for 
this. 

The receiver must provide both orthogonal-circular polari- 
zations simultaneously at X-band. Two identical channels are 
required, each having a 10 MHz bandwidth centered on a fre- 
quency of 8495 MHz. A system noise temperature of less than 
20 K is preferable for this work. A similar S-band capability 
should also exist although it is not mandatory for the pro- 
posed radar experiments. 

Some experiments require that stations other than DSS-14 
be used during an observation. To this end, interstation links 
must exist between DSS-14 and both DSS-12 and DSS-13. 
These links must be capable of supporting the PN signal 
described above. 

Accurate timing signals must be provided to the acquisition 
system. Most of the work can be supported by rubidium 
clocks. However, any observations that will use more than a 
single station will require hydrogen clocks. 


VII. Conclusion 

The High Speed Data Acquisition System will serve the 
needs of the DSN’s Planetary Radar Program at the Goldstone 
facility over the next decade. The system will be useful in 
other areas of research including radio astronomy, SETI, and 
real time antenna arraying. The ability to control the system in 
a high level software language will enable it to be used by a 
larger group of people than the current system, and this flexi- 
bility will also facilitate its use as a tool for advanced receiver 
and demodulator development. 
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Fig. 1. High Speed Data Acquisition System overall block diagram 
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Fig. 2. Time table of anticipated radar experiments 
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Fig. 3. The VAX 11/780 configuration for the HSDAS. Shaded corners represent equipment that is installed and 

working at JPL. 



Fig. 4. General purpose data acquisition hardware for the HSDAS. Shaded 
corners represent hardware that is already installed. 
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S. S. Brokl 

Communications Systems Research Section 


In support of nonflight project data acquisition, a new High-Speed Data Acquisition 
System was designed. The demodulator and accumulator modules are at the heart of this 
new system and are described in this document. 


I. Introduction 

The High-Speed Data Acquisition System (HSDAS) is cur- 
rently being implemented principally to support the DSN Plan- 
etary Radar Program (Ref. 1). In the HSDAS, the Polynomial- 
Driven PN Generators (Ref. 2) will be used to bounce a coded 
waveform off the radar target. The demodulator and accumu- 
lator subsystems are to be used to detect the reflected signals 
and produce ranging information for further processing. 

Each demodulator consists of 16 TRW TCC1023J digital 
correlator chips configured as 4 by 1 bit auto or cross corre- 
lators, with 256 lags. Included on the module is all of the 
logic and adders required to represent the 256 lag data as a 
12 bit offset binary number. 

The accumulator module follows the demodulator and 
allows accumulations of 12 bit offset binary numbers from 
each lag to a maximum 28 bit number. Each accumulator 
module supports two demodulators. For the final system there 
will be 16 demodulator and 8 accumulator modules. 

II. The Demodulator 

Figure 1 is a block diagram for the demodulator. There 
are 16 TRW TCC1023J correlators which are configured as 
4X 1 auto or cross correlators. The autocorrelation function is: 


256 T 1 

C M = £ c « K J 

2=0 L J 

where G s is a one bit representation of the original data, G s+n 
is a four bit representation of the original data time-shifted, 
and £ is the lag number. 

The cross-correlation function is: 

C(n) ^2 r^2+nl 

2 = 0 L J 

where G e is a binary signal, H &+n is a four bit representation 
of the returned signal time-shifted, and £ is the lag number. 

The normal modes of operation for the demodulator in the 
radar acquisition system will be as a cross correlator, when 
doing delay Doppler gating using the Polynomial Driven PN 
Generator, or as an autocorrelator when doing continuous 
wave (CW) detection or radio astronomy. 

It will be possible to configure several of the modules 
within the system to run as a 512 lag or 1024 lag correlator, 
with both real and imaginary channels. It may even be possible 
to run 8 by 1 bit correlations. 



The correlator chips accomplish the multiply with 
exclusive-nor logic (Ref. 3). The truth table is shown below. 


have also been operated successfully in a system test with the 
Polynomial Driven PN Generators. 


A B A B 

0 0 1 

0 1 0 

1 0 0 

1 1 1 

This can be thought of as a multiplier if one replaces 0 by - 1 
and 1 by +1. This implies that the A/D converter data input to 
the demodulator should be expressed in offset binary. A con- 
version table for a four bit 2’s complement to offset binary 
number is shown in Table 1. As can be seen, the transforma- 
tion only requires toggling the sign bit. 

In the demodulator Fig. 1 there are 16 correlator chips 
broken into 4 groups of 4 chips each. The adder tree that fol- 
lows sums the results of the correlators into a 12 bit number 
representing a running total of correlations for four bits over 
256 lags. 


III. The Accumulator 

Figure 3 is a block diagram of one half of the accumulator 
module. Each module contains two complete 28-bit wide 
accumulators. Each accumulator has two buffer memories 
(double buffered), so that one memory may be unloaded by 
the computer while the other memory is accumulating data. 
The accumulation function is: 


A W 


M 

m = 1 


where 2 = 1,256 is the lag number and M is the number of 
accumulations. M is programmable from 1 to 2 16 . 


The purpose of the accumulator is to allow a correlated sig- 
nal to be integrated for a period m times longer than the 
storage length of the correlator. 


The correlators also have a mask input so that any number 
of bits may be masked off and will not be used as part of the 
total output calculation. This is extremely useful when run- 
ning short codes since it prevents the multiple copies of 
unused lags from being summed. 


The control of the accumulator is derived in the computer 
interface and formating module which is being designed. The 
accumulator is designed with Fairchild FAST tm logic and 
30 ns access time memories. The accumulator will run with a 
maximum clock rate of 10 MHz. 


The demodulator uses one master clock which can run as 
high as 20 MHz. The code, mask, and signal have their own 
enable inputs to the module. This allows the clock to bring the 
data as required. 

The correlator chips contain two shift-register chains with 
one isolated from the exclusive-nor logic by a register which 
can be parallel loaded with the 64 bit serial bit stream in one 
of the shift-register chains. The signal LDR into the demodula- 
tor does this function for 256 lags. This allows a new 256 bit 
code string to be serially loaded while correlating the first 
stream, which is being held in the registers. 

Figure 2 is a photograph of the demodulator module. The 
modules for the radar system have been built and tested. They 


Figure 4 is a photograph of the accumulator module. The 
design and construction are complete. 

IV. Conclusion 

The demodulator and accumulator modules perform the 
detection and prefiltering of data which will be further pro- 
cessed by a Floating Point Systems FPS 5210 Array Processor. 
The array processor is also used to form the complete sums 
when multiple units are used to extend the number of lags or 
the number of bits in the product table. The High-Speed Data 
Acquisition System was designed for High Resolution Radar 
Astronomy. The system has applications in astronomy and 
possible application for the Search for Extraterrestial Intelli- 
gence (SETI) and perhaps RFI monitoring. 
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Table 1. Conversion table for a four bit 2’s complement to offset 
binary number 


Range 

2’s Complement 

Offset Binary 
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1 
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1 

0 

0 

1 

0 

0 

0 

1 

-MAX 

1 

0 

0 

0 

0 

0 

0 

0 
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Fig. 1. Demodulator block diagram 













Fig. 2. Demodulator module 
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This article describes a method for calculating the increase in received telemetry signal 
power required to compensate for the use of a radio frequency interference (RFI) filter 
in front of the DSN receiving system low-noise amplifier. The telemetry system for which 
the degradation is calculated is an uncoded PCM/PM system in which the NRZ data 
directly modulates the carrier at a modulation level which leaves a discrete carrier. A 
phase-locked loop in the RF receiver tracks the discrete carrier and coherently demodu- 
lates the PCM data. The RFI filter may be a series of Butterworth, Tchebychev, or Bessel 
low-pass, high-pass, band-pass, or band-reject filters, each with arbitrary bandwidth, 
number of poles, and, for band-pass or band-reject filters, resonant frequency. The only 
restriction is that the RFI filter must have only simple poles. Numerical results are 
presented for the RFI filters the DSN plans to place in front of the S-band FET, S-band 
maser, and X-band maser low noise amplifiers. The main conclusion is that the filters will 
produce negligible degradation at data rates below 4 Mbps. 


I. Introduction 

In conjunction with the Networks Consolidation Program 
(NCP), the DSN is developing band-pass filters to be placed in 
front of the S-band FET (field-effect transistor) low-noise 
amplifier (LNA) and the S-band and X-band maser LNAs. The 
purpose of these filters is to avoid LNA saturation by out-of- 
band radio frequency interference (RFI) signals. The DSN 
itself produces some of this interference, such as that from the 
DSN transmitters. The remainder comes from other systems 
operating in adjacent frequency bands. Reference 1 reports the 
results of a study to define the requirements for these filters. 

The use of band-reject filters has also been considered. Such 
filters could be used to eliminate both in-band and out-of-band 
RFI signals that have fixed frequency and narrow bandwidth. 

An unavoidable consequence of the use of these filters will 
be the degradation in telemetry performance caused by the 


telemetry signal distortion produced by these filters. The term 
“telemetry degradation” as used in this report is the increase 
in received signal power necessary to compensate for the 
addition of the RFI filter. Thus, the degradation will depend 
on the basis (i.e., bit error probability) on which telemetry 
performance is measured. The degradation will be greatest for 
very wide-band telemetry channels operating near the edge of 
the RFI filter pass-band. The objective of this analysis is (1) to 
develop a method of calculating this telemetry degradation for 
a variety of RFI filters, and (2) to apply this method to the 
RFI filters being developed by the DSN. 

The telemetry system considered in this report is an un- 
coded PCM/PM system. A very simplified block diagram of 
this system is shown in Fig. 1. NRZ (non-return-to-zero) data 
directly phase modulates the RF carrier at a modulation level 
which leaves a discrete carrier component. At the DSN, a 
phase-locked loop in the RF receiver, which follows the RFI 
filter and LNA, tracks the received carrier and demodulates 
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the uncoded NRZ PCM signal. The RF receiver is followed by 
an integrate-and-dump circuit, which integrates the RF receiver 
output over one bit time. A bit synchronizer establishes the 
timing for the integrate-and-dump circuit. At the end of the 
integration time for each bit, a binary decision is made as to 
whether a PCM “0” or “1” was transmitted, based on the 
polarity of the integrate-and-dump circuit output. 

The primary alternatives to the telemetry system analyzed 
in this report are (l)an uncoded PCM/PM system where the 
carrier is completely suppressed and a Costas loop recovers the 
carrier reference, and (2) an uncoded PCM/PSK/PM system 
where the NRZ PCM data phase shift keys a subcarrier and the 
resulting signal phase modulates the carrier. While some work 
was done on both of these systems, it was decided to concen- 
trate on the discrete carrier PCM/PM system because this is the 
system most likely to be used, at least initially, by the DSN for 
the very wideband telemetry channels for which significant 
degradation is expected from signal distortion in the RFI 
filter. 

While this analysis considers only uncoded PCM/PM chan- 
nels, it would be desirable to analyze systems such as Viterbi 
or Reed-Solomon/Viterbi channels, where the raw data are 
coded and the resulting coded symbol stream phase modulates 
the carrier. However, with the exception of Manchester cod- 
ing, to which the results obtained herein can be easily ex- 
tended, an analysis for coded channels, similar to that pre- 
sented here for uncoded channels, does not appear feasible. 
Note, however, that the results presented here are equally valid 
for coded channels, as long as the channel symbol error prob- 
ability is used as the measure of telemetry channel perfor- 
mance and the probability of symbol sequences in the coded 
PCM is the same as that in random data. 

This analysis employs two important approximations. First, 
the analysis assumes that the RF receiver phase-locked loop 
tracks the received carrier with some constant phase error $. 
The RF receiver phase-lock loop phase errors caused by ther- 
mal noise and the portion of the telemetry sidebands that fall 
within the bandwidth of the RF receiver phase-locked loop are 
neglected. Thus, any interaction between the signal distortion 
caused by the RFI filter and the time-varying RF receiver 
phase-locked loop phase error will not be uncovered by this 
analysis. Second, this analysis does not determine the actual 
displacement in bit synchronizer bit timing caused by the 
signal distortion. The analysis supplies a method of calculating 
the degradation as a function of the displacement in bit 
synchronizer timing. A first estimate of the displacement in bit 
synchronizer timing is made using the slope at zero frequency 
of the phase response of the RFI filter’s baseband equivalent, 
and the telemetry degradation is calculated for this bit timing. 
This timing displacement will be that seen by very narrowband 


telemetry channels. Then, using this estimate of bit synchro- 
nizer timing displacement as an initial value, the bit synchro- 
nizer timing which minimizes the telemetry degradation and 
the degradation for this bit synchronizer timing is calculated. ■ 
The actual bit synchronizer timing displacement may not be at 
the optimum value. Of course, it would always be possible to 
manually adjust bit synchronizer timing to the optimum value, 
but this may not be practical for normal telemetry operations. 

II. Analysis 

The objective of this analysis is to provide a method of 
computing the telemetry channel performance degradation 
caused by the RFI filter. Each of the subsections into which 
this section is divided provides the equations necessary for one 
step in the computational process. Because of the complexity 
of the problem, no attempt will be made to analytically 
combine the results presented in these subsections. 

A. RFI Filter Poles and Residues 

A basic assumption of this analysis is that the RFI filter has 
only simple poles. The transfer function of such a filter can be 
expressed in the form 

N Qu 

m - Qo * £ t — (1) 

k = \ k 

where Q n is the filter response at infinite frequency,/^ is the 
kth pole of the filter, and Q k is the residue at that pole. The 
first step in calculating the telemetry degradation caused by 
the RFI filter is to calculate the response at infinity and the 
poles and residues for the RFI filter. 

This article considers an RFI filter consisting of an arbitrary 
number of cascaded elements. Since the RFI filter has only 
simple poles, each of its elements must have only simple poles 
and therefore a transfer function of the form shown in Eq. (1). 
Each of these elements can be a Butterworth, Tchebychev, or 
Bessel low-pass filter or the high-pass, band-pass, or band-reject 
filters which can be synthesized from these low-pass filters 
using simple frequency transformations. Thus, this analysis 
provides the capability for assessing the telemetry degradation 
produced by using both a band-pass RFI filter for rejection of 
out-of-band RFI and one or more band-reject filters for elimi- 
nating in-band RFI. 

The first step in calculating the RFI filter response at 
infinity, its poles, and the residues at these poles, is to calcu- 
late these parameters- for the low-pass equivalent of each RFI 
filter element. The response at infinity and the poles and 
residues for RFI filter elements that are high-pass, band-pass, 
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or band-reject filters will be calculated from the response at 
infinity and the poles and residues of these low-pass equiva- 
lents. The low-pass equivalent’s frequency response will be 
adjusted to have unit response at zero frequency and either 
amplitude response A at some frequency f A = (2 rr) _1 oj a or 
noise bandwidth B N and will have zero response at infinite 
frequency. The transfer functions of the Butterworth, Tcheby- 
chev, and Bessel low-pass filters considered in this article have 
no zeros. Therefore, if P Q , . . . , P N are the poles of the 
low-pass equivalent, the residue at the fcth pole will be 


e * = - ? *nRr 

2=1 

fi=£/c 


( 2 ) 


1. Butterworth low-pass filters. The poles of a Butterworth 
low-pass filter with N poles will be 


P k = " sin 


ico D cos 

ts 


j 2 *- ■) fv] 
_("->) 5v]. 


1 <k<N 


(3) 


where f B = (2tr) 1 co B is the half-power bandwidth of the low- 
pass filter, f B can be calculated from f , and A or from B 


using 


or 


f B = (a~ 2 - 1 r l/2N f A 


, B .&) 

/„ = B., sin 

\27V/ 


'N 


(4) 


(5) 


2. Tchebychev low-pass filters. The poles of a Tchebychev 
low-pass filter with N poles will be located at 


P R = - sinh 


- ioJ T cosh 


inhjyv -1 sinh 1 (e 1 >] sin £(2 k- 

[AT 1 sinh -1 (e“ 1 ) ] cos |(2fc- 1)^J, 


Kk<N (6) 


where f T = (2 it) 1 oj t is the cut-off frequency of the Tcheby- 
chev filter and 


e = (p 2 - 1) 


1/2 


(7) 


where 20 log 10 (p) is the Tchebychev filter ripple factor. 
When f A is specified, f T can be computed using 


f. 


T coshj^ 1 cosh 1 [e 1 (A 2 - 1) 1/2 ]} 

1 


fj, N = odd 


' f A » 


coshj/V^ 1 cosh 1 [e 1 (A 2 (1 + e 2 )- l) 1 / 2 ]} 

N = even (8) 

When B N is specified, the poles and residues,/*^ and Q' k , 1 < k 
< N, of a Tchebychev filter with unit cutoff frequency (f T 
= 1) are calculated. Then the poles and residues of the desired 
filter are 


B 


p k = #n> 


B 


(9) 


N 


and 




(10) 


N 


where 


B 


N 




k~l 


k ^2 

P' +P’ 

k = 1 2=1 \ k » 


(ID 


is the noise bandwidth of the Tchebychev filter with unit 
cutoff frequency. 

3. Bessel low-pass filter. An yV-pole Bessel low-pass filter 
with unit response at zero frequency has transfer function 


H u ,(s) = 


NO 

P N (cs) 


( 12 ) 


where 


Pf/W = E a Nk 


(13) 


k = 0 
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is the Mh order Bessel polynomial and c is a parameter which 
is adjusted to give either the required bandwidth f A at attenu- 
ation A or the required noise bandwidth B N . The ATh order 
Bessel polynomial is defined by the recursion formula 


/'( co) = -2A 2 < 


N 


i + £(-i) fc/2 ^ 


k-2 
A:=even 


Nk k 
co 

NO 


P N (s) = (2N- l)p N _ l (s) + s 2 p n _ 2 (s) (14) 


where p 0 (s) = 1 andpj(s) = s + 1. The corresponding recursion 
relationships for the Bessel polynomial coefficients are 


N a 

k= 2 a N0 

fc=even 


a N0 ^ U N~ 1,0 

(15) 

a Nl = ( - 2N - ^ a N-l,l 

(16) 


Nk 


{2N- 1 )a 


N- 


1 ,k + a N-2,k-2’ 


2<k<N - 1 


(17) 


and 


AW 


= A 


N-2.N-2 


(18) 


Used recursively, Eq. (12) and Eqs. (15) through (18) provide 
a simple method of calculating (for c = 1) the coefficients in 
the transfer function of a Bessel filter of any number of poles. 


When f A and A are specified, one first adjusts c to obtain 
the desired bandwidth and then employs a polynomial root 
funding subroutine to calculate the filter poles (the roots of 
p N (cs)/a N0 ). First, one finds the frequency/^ = co' A at 
which \a N0 lp N (s)\ is equal to A. This is equivalent to finding 
the zero of 


/(co) = 1 - A : 


i + £(-D* /2 

a m 

k= 2 


Nk 

'no 


*=eve n 


V (- i)(*-D/2 a Jil ^ 

C=1 a N0 


k= 1 
/c = odd 


(19) 


k = 1 a N0 

/c=odd 


N 


£ (-D (fe 

k= 1 
k=odd 


-11/2 U Nk 
°N0 


k co 


k 1 1 


( 20 ) 


Having determined cojj such that /(eo^) = 0, 


c 



( 21 ) 


When the noise bandwidth B N is specified, a polynomial 
root-funding subroutine can be used to calculate the poles P' K , 
1 < k <N, of a NQ /p N (s). Then, the residues Q' k of a N0 /p N (s ) 
at these poles can be calculated using P' K instead of P K in Eq. 
(2). Finally, Eqs. (9), (10) and (11) above can be used to 
calculate the poles and residues of a Bessel filter with the 
desired noise bandwidth. 


4. High -pass filters. A high-pass filter can be synthesized 
from a low-pass filter using the transformation 

( v ) <22) 

If the low-pass filter has unit response at zero frequency, 
amplitude response A at frequency f A , and zero response at 
infinite frequency, the high-pass filter will have zero response 
at zero frequency, amplitude response A at frequency f A , and 
unit response at infinite frequency. If P K and Q K are the poles 
and residues of the low-pass filter, the response at infinity of 
the high-pass filter is 


Go 


£ 




Qi 


(23) 


This can be solved using a series ^f linear extrapolations. To 
use this method an expression for / (co) = dI(<jo)/doj is needed. 
Differentiating Eq. (19), 
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The location of the fcth pole of the high-pass filter will be and 



(24) 



( 30 ) 


and the residue at that pole will be 


6. Band-reject filters. A band-reject filter can be synthe- 
sized from a low-pass filter using the transformation 



« a 

~^Q k 

k 


(25) 


~ ^ IP 




(31) 


5. Band-pass filters. A band-pass filter can be synthesized 
from a low-pass filter using the transformation 



where f R = (2tt)~ 1 oj r is the band-pass filter resonant fre- 
quency. The resonant frequency will be approximately equal 
to the filter center frequency. If the low-pass filter has unit 
response at zero frequency and zero response at infinite fre- 
quency, the band-pass filter will have zero response at zero 
frequency, unit response at f R , and zero response at infinite 
frequency. Furthermore, if the low-pass filter has attenuation 
A at frequency f A , the band-pass filter will have bandwidth 
f A between points at which the band-pass filter attenuation 
is ^4. 


If the low-pass filter has unit response at zero frequency and 
zero response at infinite frequency, the band-reject filter will 
have unit response at zero frequency, zero response at f R , and 
unit response at infinite frequency. Furthermore, if the low- 
pass filter has attenuation A at frequency f A , the band-reject 
filter will have bandwidth f A between points at which the 
band-reject filter attenuation iszl. 

If P K and Q k , 1 < k <Af, are the poles and residues of the 
low-pass filter, the response at infinity of the band-reject filter 
will be 

N Q, 

k= 1 k 

and the band-reject filter poles will be 


If P K and Q k , 1 < k <Af, are the poles and residues of the 
low-pass filter, the response at infinity of the band-pass filter 
will be zero and the band-pass filter poles will be 


K = 2 P * +iu * 
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1 <k<N 


(34) 


P k ~ ~2 P k~ ,C0 R 
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1/2 


1 <k<N 


(28) 


Note that the band-pass filter has 2N poles. The residues at 
these poles will be 


Note that the band-reject filter has 2 N poles. The residues at 
these poles will be 


and 


^ / K \ 

a >-- n a ‘ V - n ) 


(35) 


p" v-i 

Qk = Q k \ '-y- 


(29) 


ai~ -—oji-y) 


(36) 
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7. Cascaded elements. Having developed equations for cal- 
culating the response at infinity, the poles, and the residues at 
these poles for the individual elements of the RFI filter, the 
next step is to calculate the response at infinity, the poles, and 
the residues at these poles for the combination of a series of 
such elements. Consider two filters, one with poles and 
another with N 2 poles, which have no common pole. Let Q 0I 
be the response at infinity and P 1K and Q 1K be the poles and 
residues of the first filter and Q 02 be the response at infinity 
and P 2K and Q 2K be the poles and residues for the second 
filter. The cascade of the two filters has N 1 + N 2 poles, P IK’ 
1 < k < vVj , and P 2K , 1 < k < N 2 . The response at infinity of 
the cascaded filter will be 


and rj, 0 < rt < 1, is the fraction of this noise introduced after 
the RFI filter. 

Referring again to Fig. 2, the signal and noise components 
of the RFI filter input will be 

s 0 (t) = (2 P) 1/2 sin [oj c t + <P C + flx(t)] (40) 

and 

n Q (t) = n Ql (t) 2 1/2 sin (c o c t + 4> c - <P) 


Q 0 ^01^02 


(37) 


+ n o Q (0 2 1 12 cos (u c t + <p c -<t>) (41) 


The cascaded filter residue at P K will be 
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and that at P K will be 


^2 k ^2 k 
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^oi + S 
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12 
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Starting with two RFI filter elements, the above process can 
be repeated until the response of the RFI filter at infinity, the 
RFI filter poles, and the residues at these poles have been 
calculated. 


P is the received signal power, f c = (2jr) -1 aJ c is the carrier 
frequency, and <t> c is the carrier phase. 0 is the telemetry 
modulation level andx(r) is the binary-valued (±1) telemetry 
modulation. This analysis assumes x(t) is zero. n QI (t ) and 
n 0 Q(t) are statistically independent gaussian processes with 
constant one-sided power spectral density (1 - t))N 0 over the 
frequency band occupied by the telemetry data signal. Note 
that for analytical convenience the quadrature noise compo- 
nents cf n 0 (t ) are referenced to a carrier with instantaneous 
phase c o c t + <p c - 0, where 4> is the static phase error in the RF 
receiver phase-locked loop. 

For the input signal and noise specified in Eqs. (40) and 
(41), the signal and noise output of an RFI filter with impulse 
response h{t), amplitude response A (to), and phase response 
<p( to) will be 


s { (t) = (2 P) U2 cos (P)A (to c ) sin [to c r + <P C + <t>( oo c )] 


B. RF Receiving System 

Figure 2 shows the analytical model for the demodulation 
of the telemetry data signal from the RF carrier. Note that 
noise is introduced at two different locations in Fig. 2. n 0 (t) 
represents the noise, such as atmospheric and antenna noise, 
introduced prior to the RFI filter. n 2 (t ) represents the noise, 
such as the FET or maser LNA noise, introduced after the RFI 
filter. Two noise sources are required because the RFI filter 
will reduce the effective telemetry detector noise bandwidth 
for noise introduced prior to the RFI filter, but not for noise 
introduced after the RFI filter. The analysis assumes that the 
one-sided power spectral densities of n 0 (t) and n 2 (t) are 
(1 - r?)/V 0 and r/N 0 where N 0 is the one-sided power spectral 
density of the receiving system noise produced by all sources 


+ (2P) i/2 sin (P) cos [oj c t + 0 C + <t>( oj c )\ 



h c (r) x(t - t) dr 


+ (2P) l/2 sin (/3) sin [oj c t + <t> c + 4>(oj c )} 



h s (r)x(t - r) (It 


(42) 


109 



and 


Hj(0 = 2 1/2 sin [c o c t + <t> c + 0(to c ) - 0] 


Assuming the RF receiver phase-locked loop tracks the 
carrier component of the RF1 filter output with phase error 0, 
the reference signal for carrier demodulation shown in Fig. 2 
will be 


X f h c (T) n Ql (t - t) dr 


- 2 1/2 cos [w c r + 0 C + 0(cu c ) - 0] 


s 3 0) = 2 1/2 cos [co c t + 0 C + 0(w c ) - 0 ] (47) 


Thus, neglecting terms with zero frequency and terms in 
2 co c t, the signal and noise outputs of the carrier demodulator 
will be 


xf h s (r) « 0/ (t - r) dr 


s(t) = P 1 ' 2 sin (0) 
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[/z c (r) cos (0) 
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h c (j) = h(r) cos [go c r + 0(ce> c )j 


and 


co 
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(43) 



^(f) n QI (t - t) dr 


(49) 


Using the above results, Fig. 3 shows a simplified analytical 
model for the DSN RF receiving system and telemetry 
detector. The transfer functions H'(s),H c (s), and H s (s) shown 
in Fig. 3 are those corresponding to the impulse responses 
(44) h'(T), h c (j), and h s (r) where 


h s (T) = h(r) sin [w c r + 0(w c )l (45) 


h' (t) = h c (T) cos (0) + h s (r) sin (0) (50) 


The noise n 2 (t ) introduced after the RFI filter will have the 
form 


For an RFI filter having a transfer function of the form 
shown in Eq. (1), 


n 2 (0 = n 2 j(t) 2 1/2 sin [w c f + 0 C + 0(co c ) - 0| 


+ « 2e (r)2 1/2 cos [w c r + 0 c + 0(w c )-0| (46) and 

where, as in Eq. (41), the carrier instantaneous phase to which 
the quadrature noise components are referenced is chosen for 
analytical convenience. 
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where 


P Ck P Sk P k * 0j c> I m (P k ) > 0 

= P k +i(jJ c , Im (P k ) < 0 (53) 

where Im(z) denotes the imaginary part of a complex variable 

z, 

Q Ck = \ Q k ex P [-*0(w c )], Im(P k ) > 0 

= \ Q k exp [+ J0(w c ) ] , Im(P k ) < 0 (54) 

2^- = - Yi Q k ex P [- i(p(oj c ^ > Im{P ^ > 0 

= ^jQ k exp [z0(oj c )] , /m(P fc ) < 0 (55) 

Q co = Q o cos ( 56 ) 

and 

Q so = Q o sin W« c )l ( 57 > 

Thus, 

«'M = c; + E ~7 (58) 

fc = 1 r k 

where 

K= p ck =p sk’ m 

and 

Q' k = Q Cfe cos (0) + sin (0), 0 < * <^V (60) 

Equations (53) through (57) and Eqs. (59) and (60) provide a 
means of calculating the response at infinity and the poles and 
residues of the transfer functions H’(s), H c (s), and H s (s) 
shown in Fig. 3, given the response at infinity and the poles 
and residues of the RFI filter. 

The effect of the noise filtering shown in Fig. 3 on telem- 
etry system performance will be to increase telemetry perfor- 


mance by 1/6 jy where 6 V is the factor by which the variance 
of the noise at the integrate-and-dump circuit output is 
reduced by the filtering. The next section of this report 
provides a means of calculating the ratio of the bandwidth of 
an integrate-and-dump circuit preceded by a filter with only 
simple poles to that of the integrate-and-dump circuit alone. 
Thus, if is this ratio for the filter with transfer function 
H c (s ) and 5 ns is this ratio for the filter with transfer function 
H s (s), 

8 n = ri + (\- V )(5 nc + 8 ns ) (61) 

Thus, as expected, if all of the receiving system noise is 
introduced after the RFI filter, r?= 1 and 8 N = 1. When the 
noise is introduced after the RFI filter, the RFI filter does not 
reduce the effect of the noise on telemetry performance. 

C. Reduction in the Telemetry Detector Noise 
Bandwidth Produced by a Filter Preceding the 
Integrate-and-Dump Circuit 

As discussed in the preceding section, a method is needed 
for calculating the reduction in telemetry detector noise band- 
width caused by a filter with only simple poles preceding the 
integrate-and-dump circuit. Since the filter has only simple 
poles, its impulse response will have the form 

N 

V r ) = G 0 S(r) + X) 2* ex P( F * 7 ) ( 62 ) 

k = l 

where Q 0 is the filter response at infinite frequency, P K and 
Q k , 1 < k < N, are the filter poles and residues, and 6 (r) is the 
Dirac delta function. The impulse response of the integrate- 
and-dump circuit will be 

h 2 (j) = 1, o < t < r, 

= 0, t>T (63) 

where T is the integration time. The impulse response of the 
cascaded filters will be the convolution of h^(t) and h 2 (t ) or 

N Qk 

h (f) = Q q + ^ — [exp (P fc r) - 1 j , 0 < r < T 

k= l k 

N Q k 

= ^ -p- [1 -exp(-R fc r)] exp (P k r), 
k - 1 

(64) 
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The one-sided noise bandwidth of a filter with impulse 
response h(j) will be 


B 


N 


2 



h 2 (j) dr 


(65) 


The noise bandwidth of the integrate-and-dump circuit above 
is 772. Thus, using Eq. (64) in Eq. (65), the reduction in noise 
bandwidth caused by the filtering in front of the integrate-and- 
dump circuit will be 



The next step is to calculate each of the 2 Af+2 possible 
output signal levels. Refer to the simplified block diagram in 
Fig. 4 and the sample x{t) waveform in Fig. 5. Note that the 
«th data bit inx(t) is received during [(X + n - 1 )T, (X + n)T\. 
Assume that because of the distortion of x(t) the bit synchro- 
nizer causes the integrate-and-dump circuit to integrate over [(X 
+ n - 1 )T, (X + n)T\, during detection of the nth data bit, 
where X = X + e. Then 


(X+a-1 


(X+n + e)?’ 

dt 

(X+n-l+e)r 


(67) 


or 


J - ( X+n)T /* 

y(t)dt + I 

0+H-l+e)r J(\+n)T 


( \+n + e)T 

y(t)dt , 0 < e < 1 (68) 


and 


( 66 ) 

Having computed the response at infinity and the poles and 
residues for the transfer functions H c (s ) and H s (s), using the 
equations derived in the previous section, Eq. (66) with appro- 
priate substitutions, can be used to calculate & NC and S NS . 


~(\+n-l)T -n + c 

I y(t)dt + I y(t)c 

J(\+n-l+e)T J(\+n- 1)' 


(A. + « + e)r 
dt, 
)T 


■ l<e<0 (69) 


If h'(r), the impuse response corresponding to the transfer 
function has the form 


D. Probability Density of the Integrate-and-Dump 
Circuit Output Signal Level 

As a result of the filter with transfer function H'(s) shown 
in Fig. 3, the signal integrated by the integrate-and-dump cir- 
cuit will depend not only on the data bit being detected, but 
also on preceding bits and, in this analysis, on one succeeding 
bit. The dependence on preceding bits occurs because the 
impulse response of H'(s ) will, in general, be of infinite dura- 
tion; that is, the current output of H'(s) depends on all 
previous inputs. The dependence on the succeeding bit is a 
consequence of the displacement in bit synchronizer timing 
caused by H'(s). This analysis assumes the displacement in bit 
synchronizer timing does not exceed one bit. 


N 


h'(r) = Q q 5 (r) + Q k exp(P k T) 


(70) 


k = l 


N 


Qu 


y(t) = x n ^2 -5- [ ex P i p k [f-(^ + «- 1 )?"] ) - 1] 

k = 1 * 


n-1 


E X m TjT 


k 
) 

&=i k 


X [exp {P [t - (X + m - 1) T\] 


Each different signal waveform that can occur during a bit 
integration can result in a different integrate-and-dump circuit 
output signal level. This analysis assumes that the impulse 
response of H'(s) is significant for only M bit times. If the 
impulse response extends over M preceding bits, and one 
succeeding bit is significant because of bit synchronizer timing 
displacement, there are 2 M+2 possible integrate-and-dump cir- 
cuit output levels. This analysis assumes that all data sequences 
occur with equal probability. Hence, all of the 2 M+ 2 possible 
integrate-and-dump circuit output signal levels occur with 
equal probability 1/2 M+2 . 


-exp {P k [t- (X + ffOr]}], 

(\ + n-l)T<t < (X + n)T (71) 

Using Eq. (71) in Eqs. (68) and (69), 

00 

z = c x , x = ±1 (72) 

n m n-m n-m v 7 

m~~ 1 
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where 


and 


(1-e) 


N 


Qo- E 


0, 


fe= 1 


c 


m 


■ N Q , 

^ {exp [(m + 1 + e^T] - 2 exp [(m + e^T] 
k= 1 


+ 


E 




Q k 

-r— {exp [(1 + e)F r] 

p2 p * 

k 


+ exp [(m - 1 + e)/^]}, -l<e<l,w>2 (79) 


- 2 exp (eP fc T) + 1} , 0 < e < 1 (73) 


(1+e) 





E 



+ 


TV 


E 


[<l 'lei/',/] -I 

“k 1 


- 1 <e<0 


(74) 


c 


-l 


N 


e»- E 


k = 1 



+ 2 ~~ [exp (eP k T) - 1 ] , 0 < e < 1 (75) 

c_j = 0, -1 <e<0 (76) 


Thus, approximating the infinite summation in Eq. (72) with 
one whose upper limit is M, the probability density function 
of the integrate-and-dump circuit output signal level, normal- 
ized to P 1 ! 2 sin (p), is 


P(z„) 


1 

2 m+ 2 




=± 1 


X 5 




(80) 


E. Telemetry Performance Degradation 

Having determined the probability density function of the 
integrate-and-dump circuit output signal level and the reduc- 
tion in the noise bandwidth of the telemetry detector, it 
remains to calculate the telemetry performance degradation. 
In a system with no filter or RF receiver phase error (0 = 0), 
the channel bit error probability will be 


N 

= E ~T W [(2 + e)P k T 1 
*“1 f*T 



(81) 


-2 exp [(1 + e)P k T\ + exp (eP k T)}, 0 <e < 1 

(77) 




where E B is the received signal energy per bit and iV 0 is the 
receiving system noise spectral density. In such a system, z n 
and & N will be one. When z n and b N are not one, the effect on 
the bit error probability is equivalent to that of reducing 
( e b/N 0 ) 1/2 by z N^jJ 2 ■ Therefore, in the presence of signal 
distortion and RF receiver phase error, the bit error probabil- 
ity will be 


+ 



A 


{exp [(2 + e)P k T] 





1/2 


(x • c) 


- 2 exp [(1 + e)P k T] + 1}, -l<e<0 (78) 


(82) 
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x is the vector (x n _ M , . . . ,x n _ 1 , l,x„ +1 ),cis the vector (c M , 
,), x • c is the dot product of the two vectors, and 

E 

X 

denotes the summation over all the 2 M+1 possible values ofx. 
Note that to simplify the above analysis (without any less of 
generality) the «th transmitted bit was assumed to be a binary 
zero ( x n = 1). If Pb is the allowable telemetry channel bit 
error probability and ( E B /N Q ) B is the E B /N 0 that yields this 
bit error probability in the absence of filtering and RF receiver 
phase error, the telemetry channel performance degradation 
8 t will be the solution of 


the performance degradation 8 T will be a function of e, the 
displacement in bit synchronizer timing expressed as a fraction 
of the bit time. While analytical prediction of the bit syn- 
chronizer timing displacement appears feasible, such an analy- 
sis would be quite difficult. In this report two alternate 
approaches are used. 

1. Narrowband approximation. The first and simplest 
approach is to calculate an estimate of e from the phase 
response of the filter with transfer function H'(s ) shown in 
Figs. 3 and 4. For narrowband telemetry modulation, the 
displacement in bit synchronizer timing will be equal to the 
time delay through H'(s) at very low frequencies. Thus, if 


1 



1 -Erf 




(83) 


H'(ico ) = A '(co) exp [z0'(co)] (86) 


the narrowband approximation to the displacement in bit 
synchronizer timing will be 


This equation can be solved for 8 T using a series of linear 
extrapolations starting with some initial estimate for 8 T (i.e., 
6 r = 1). Since P B will be almost an exponential function of 
5 r , convergence of this iterative process will be accelerated by 
working with the natural logarithm S.n(x ) of both sides of Eq. 
(83). Thus, the slope needed for the linear extrapolation will 
be 


_ 1 d<p'(oj) 
e i T dco 


U) = 0 


(87) 


If the response at infinity Q 0 and the poles and residues, P k 
and Q k (1 <k <jV) of H'(s ) are known, and 


dtn(P B ) _ i dP B 
d&r Pb d8 r 


(84) 


where 


dp L 

rf5~ 


N„ 


n 8 N 8 T 


1/2 


a k = [R (P k )] 2 HKP k )] 2 

(88) 

K = - 2/ (V 

(89) 

= -R(Q k )R(P k )-I(Q k )I(P k ) 

(90) 

d k = 7 (e*) 

(91) 


X 


1 

2 m+ 1 


E 


(x • c) exp 



(85) 


e k = R(Q k )/(P k )-I(Q k )R(P k ) (92) 

f k = -R ( Q k ) (93) 


where R(z) and/(z) denote the real and imaginary parts of a 
complex variable z. 


F. Bit Synchronizer Timing Displacement 

Examining Eqs. (82) and (83) and Eqs. (73) through (79), 
one notes that the telemetry channel bit error probability and 


= 1 -T(O) y' (0) - T(0) X' (0) 
€/ T [*( 0)] 2 + [ F ( 0)] 2 
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where 


where 


X(0) 



N e , 

no) = E t 


k = 1 


*'(o) = E 


N ( a k d k - b k c k ) 


fc=i 


and 



>"( 0 ) = 


E 


a Jk - 


(98) 


2. Bit synchronizer timing displacement for minimum 
telemetry degradation. The second and somewhat more com- 
plex approach is to find the bit synchronizer timing displace- 
ment which minimizes the telemetry performance degradation. 
If the bit synchronizer does not automatically adjust the 
timing to minimize the telemetry degradation, the timing 
could be manually adjusted to achieve the minimum perfor- 
mance degradation. Of course, requiring manual adjustment of 
bit synchronizer timing might create some operational 
difficulties. 




(x- c) x- — 


dc\ 2 

del 


X exp 



6 


N 




(103) 


A second iterative procedure, also using a series of linear 
extrapolations, can be employed to find the value of e (e 0 ) 
that minimizes 5 r . Finding the value of e that minimizes S r is 
equivalent to finding the zero of d& T /de. To find a zero 
crossing of a function, the iterative process requires an initial 
estimate of the solution and a means of calculating the func- 
tion and its derivative. The initial estimate used herein (ef) is 
that obtained using the narrowband approximation discussed 
above. dd T /de and d 2 d T /de 2 can be calculated using 


dd T 

de 



(99) 


and 



' 6„ /E n \ 

X exp 

T ( B \ , ^2 

“s-IaT ( x ' c ) 


L N ' °'b J 


(104) 




N Q , 

+ Ef - ( ex P K 1 + e ^k 

k = 1 k 

-2 exp (ePJ)}, 0<e<l (105) 


d 2 5 


T 


de 2 


a + b T dot _ b r a dp_ 
ft de P de ^ 2 de 

(100) 


d 2 c N 

7 = L QJ < ex P [(1 +e)PJ} 

de 

-2 exp (eP k T)}, 0<e<l (106) 
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de 


Qo - E 


N Qu 


k = 1 


£?Cj 

de 


jl Qk 

k=l k . 


N Q 


-_k 
} 

k=\ k 


+ E P~ ex P K 1 + e ^ P k T ^ > - 1 < e < 0 


(107) 


rf 2 c n " 

T = £ Q k r ^p[(l+e)P k T], -1<£<0 

k= 1 


de 2 


(108) 


N Q v 

+ E F^ {ex P [(2 + e )^ 

fc=l k 

-2 exp [(1 + e)P k T] }, -l<e<0 


c? 2 c N 

1 = E e,r{exp[(2 + e)/’ fc 7’] 


(114) 


de 2 


k = 1 


dc -i 

r = e «- £?- 

fc=i fc 


de 


N Qk 

+ E ^exp(eP A r), 0 < e < 1 


(109) 


d 2 c_, " 

— f 1 = E e fc r ex P ^ eP k n 0 < e < 1 

rfe A = 1 


( 110 ) 


6?C_ tf 2 C 


= 0, -1 <e<0 (111) 


-2 exp [(1 + e)P.T] }, -l<e<0 
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dc * Q. 

IT = E P ^ ex P [(^ + 1 + t) p J] 


k= 1 


-2 exp [(m + e )P jT ] 

+ exp [( m - 1 + e)P k T ] }, -1 < e < 1, m > 2 


(116) 


<7 2 c iv 

m — 


de 2 


E Qk T ^ ex P C( w + 1 + e)VJ 


fc=i 


-2 exp [(w + e)P k T] + exp [( m - 1 + e)P k T\ , 
-l<e<l,m>2 (117) 


dc JL Q k 

= E ~T {ex P K 2 + e K T 1 ~ 2 ex P K 1 + e K T ^ 

k= 1 k 


+ exp (eP t T)}, 0 < e < 1 


( 112 ) 


< 7 2 c , ^ 


cfe 2 


1 = E <2*E{exp [(2 + e)R k T] 


k= 1 


-2 exp [(1 + e)/> fc T] + exp ( ePJ )}, 

0 < e < 1 (113) 


III. Numerical Results 

Data are presented in this section for the telemetry channel 
degradation for the S-band FET amplifier, S-band maser ampli- 
fier, and X-band maser amplifier RFI filters. For each of these 
RFI filters telemetry channel degradation is plotted as a func- 
tion of the telemetry channel carrier frequency for telemetry 
channel data rates of 4 Mbps, 10 Mbps, and 20 Mbps. The 
4 Mbps rate is the highest data rate currently planned. The 
data for the 10 Mbps and 20 Mbps data rates were included to 
show the impact of the RFI filters on future systems. Plotting 
the telemetry channel degradation as a function of carrier fre- 
quency shows the impact of the RFI filter on telemetry chan- 
nels whose center frequency does not fall at the RFI filter 
center frequency. 
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Two sets of curves are shown for each RFI filter. One set 
assumes that all of the receiving system noise is introduced 
after the RFI filter (17= 1). The other set assumes that all of 
the receiving system noise is introduced before the RFI filter 
(r] = 0). The telemetry degradation 8 T1 that occurs when all 
the noise is introduced after the RFI filter will always be 
greater than the telemetry degradation & T0 that occurs when 
all the noise is introduced after the RFI filter, because the RFI 
filter reduces the telemetry system noise bandwidth for noise 
introduced prior to the RFI filter. Given the two sets of data 
(77 = 1 and 77 = 0) for each RFI filter, the telemetry degrada- 
tion for any other value of 77 is 

6 n 6 ro 

S T = T . rr — P*— (118) 

T (1 - v)S Tl + V& T0 

A. S-band FET Amplifier RFI Filter 

Figures 6 and 7 show the telemetry channel degradation as a 
function of carrier frequency for the S-band FET amplifier 
RFI filter. Figure 6 assumes all of the noise is introduced after 
the RFI filter (1?= 1). Figure 7 assumes all of the noise is 
introduced prior to the RFI filter (1? = 0). The results pre- 
sented in Figs. 6 and 7 assume that the RF receiver phase- 
locked loop static phase error is zero and that the allowable 
telemetry channel bit error probability is 10~ 5 . 

The data shown in Figs. 6 and 7 are for 7th and 9th order 
Tchebychev bandpass filters with 270 MHz bandwidth be- 
tween points at which the filter attenuation is -?0dB and 
0.01 dB ripple factor. The filter resonant frequency, which is 
approximately equal to the filter center frequency, is 2250 
MHz. 

The telemetry channel degradation shown in Figs. 6 and 7 
increases significantly with data rate, even when the carrier 
frequency is at the center of the RFI filter passband. The 
degradation increases significantly as the edge of the passband 
is approached. While the filter amplitude response may cause 
significant degradation, it is suspected that much of the degra- 
dation shown in Figs. 6 and 7 results from the non-linear phase 
response of the Tchebychev filter. While the degradation 
shown in Figs. 6 and 7 for the 4 Mbps data rate, especially that 
for the 9th order filter, seems acceptable, it appears that the 
use of this filter will cause significant telemetry degradation at 
much higher bit rates. This degradation will be significant even 
for telemetry channels whose center frequency is at the RFI 
filter center frequency. 

B. S-band Maser Amplifier RFI Filter 

Figures 8 and 9 show the telemetry channel degradation as 
a function of carrier frequency for the S-band maser amplifier 


RFI filter. Figure 8 assumes that all of the noise is introduced 
after the RFI filter (77 = 1). Figure 9 assumes all of the noise is 
introduced before the RFI filter (7? = 0). The results presented 
in Figs. 8 and 9 assume the RF receiver phase-locked loop 
static phase error is zero and that the allowable telemetry 
channel bit error probability is 10 -5 . 

The data shown in Figs. 8 and 9 are for a 5th order 
Tchebychev bandpass filter with 60 MHz bandwidth between 
points at which the filter attenuation is -0.5 dB and 0.01 dB 
ripple factor. The filter resonant frequency is 2292 MHz. 

As in the case of the S-band FET amplifier RFI filter, the 
telemetry channel degradation increases significantly with data 
rate, even when the carrier frequency is near the RFI filter 
center frequency. The telemetry degradation does not seem 
unreasonable for the current 4 Mbps maximum data rate, but 
the degradation at higher data rates may not be acceptable. 

C. X-band Maser Amplifier RFI Filter 

Figures 10 and 11 show the telemetry channel degradation 
as a function of carrier frequency for the X-band maser ampli- 
fier RFI filter. Figure 10 assumes that all of the noise is 
introduced after the RFI filter (77 = 1). Figure 1 1 assumes that 
all of the noise is introduced before the RFI filter (77 = 0). The 
results shown in Figs. 10 and 11 assume that the RF receiver 
phase-locked loop phase error is zero and that the allowable 
telemetry channel bit error probability is 10 -s . 

The data shown in Figs. 10 and 11 are for a 5th order 
Tchebychev bandpass filter with 160 MHz bandwidth between 
the points at which the filter attenuation is -3 dB and either 
0.01 or 0.05 dB ripple factor. The filter resonant frequency is 
8450 MHz. The 0.05 dB ripple factor curves in Fig. 11 were 
omitted for the 4 and 10 Mbps data rates, because they were 
essentially the same as the 0.01 dB ripple factor curves for the 
corresponding cases. 

As in the cases of the S-band FET amplifier RFI filter and 
the S-band maser amplifier RFI filter, the telemetry channel 
degradation increases significantly with data rate. However, 
the telemetry channel degradation caused by the X-band maser 
amplifier RFI filter seems to be appreciably less than that 
produced by the S-band RFI filters. For 4 Mbps the degrada- 
tion for 17= 1, the worst case, is less than 0.1 dB over essen- 
tially all of the 8400 to 8500 MHz frequency band. For 
10 Mbps the degradation is less than 0.2 dB over almost all of 
this frequency band. Even at 20 Mbps, the degradation is only 
about 0.25 dB for a carrier located at the RFI filter center 
frequency. Finally, the Tchebychev filter ripple factor does 
not appear to strongly affect the telemetry channel degrada- 
tion produced by the X-band maser amplifier RFI filter. 
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D. Effect of Bit Synchronizer Timing Displacement 

As noted in Section II.F of this article, this analysis does 
not calculate the actual bit synchronizer timing displacement. 
The “narrowband” approximation is used to obtain an initial 
estimate of bit synchronizer timing displacement and then an 
iterative process is employed to find the bit synchronizer 
timing displacement that minimizes telemetry channel degra- 
dation. It is this minimum value of the telemetry channel 
degradation that is plotted in Figs. 6 through 11. 

The assumption in plotting this minimum telemetry chan- 
nel degradation is that bit synchronizer timing can always be 
manually adjusted to optimize performance. However, as man- 
ual adjustment of bit synchronizer timing may be very cum- 
bersome, let us assess the impact of non-optimum bit synchro- 
nizer timing displacement by comparing both the initial 
estimate of bit synchronizer timing displacement, obtained 
using the narrowband approximation, with the optimum value 
of bit synchronizer timing displacement and the telemetry 
channel degradations for these two estimates of bit synchro- 
nizer timing displacement. 

Figure 12 compares the initial estimate and optimum value 
of bit synchronizer timing displacement, and Fig. 13 compares 
the corresponding telemetry degradations as a function of 
carrier frequency for the S-band maser RFI filter. Clearly, for 
carrier frequencies near the edge of the RFI filter passband, 
there is a substantial difference between the initial estimate 
and the optimum value of bit synchronizer timing displace- 
ment and also a substantial difference between the degrada- 
tions. The degradation for the initial estimate is significantly 
greater, for carrier frequencies near the edge of the RFI Filter 


passband, than the minimum degradation. This suggests that 
the telemetry degradation could be somewhat greater than 
that shown in Figs. 6 through 11, if actual bit synchronizer 
timing displacement differs from the optimum value as much 
as the initial estimate does and manual adjustment is not 
feasible. 


IV. Conclusion 

This article has presented a method for determining the 
effects of telemetry channel filtering on the performance of a 
discrete carrier PCM/PM telemetry channel and has applied 
this technique to calculate numerical results for the S-band 
FET amplifier, S-band maser amplifier and X-band maser am- 
plifier RFI filters being considered by the DSN. This method 
can be adapted to calculate telemetry channel degradation for 
a very wide variety of RFI filters. The method can be easily 
extended to filters other than the cascaded Butterworth, 
Tchebychev, and Bessel low-pass, high-pass and band-reject 
filters considered herein, as long as the RFI filter has only 
simple poles and its poles and the residues at these poles can 
be calculated. 

The numerical results presented for the S-band FET ampli- 
fier, S-band maser amplifier, and X-band maser amplifier RFI 
filters show that the telemetry channel degradation should be 
negligible for bit rates less than the current 4 Mbps maximum 
capability of the BBA (baseband assembly) telemetry pro- 
cessor. However, at data rates significantly greater than 
4 Mbps, the telemetry degradation produced by the RFI fil- 
ters, especially the S-band RFI filters, may be unacceptable. 
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Fig. 1. Uncoded PCM/PM telemetry link block diagram: (a) transmitting system; 
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Fig. 2. Analytical model for the RF receiving system 










Fig. 3. A simplified analytical model for the DSN RF receiving system and telemetry detector 



Fig. 4. Analytical model for determining the possible 
integrate-and-dump circuit output signal levels 
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Fig. 5. Sample telemetry data signal x(f) 
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Fig. 6. Telemetry degradation as a function of carrier frequency for the S-band FET amplifier 
RFI filter assuming ail noise originates after the RFI filter 
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Fig. 7. Telemetry degradation as a function of carrier frequency for the S-band FET amplifier 
RFI filter assuming all noise originates before the RFI filter 
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FILTER PARAMETERS: 

BANDWIDTH: 60 MHz AT -0.5 dB 
RESONANT FREQUENCY: 2292 MHz 
TYPE: 5TH ORDER TCHEBYCHEV 



2260 2270 2280 2290 2300 2310 2320 2330 2340 

CARRIER FREQUENCY, MHz 

Telemetry degradation as a function of carrier frequency for the S-band maser 
amplifier RFI filter assuming all noise originates before the RFI filter 
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FILTER PARAMETERS: 

BANDWIDTH: 160 MHz AT -3 dB 
RESONANT FREQUENCY: 8450 MHz 
TYPE: 5TH ORDER TCHEBYCHEV BANDPASS FILTER 
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BIT SYNCHRONIZER TIMING DISPLACEMENT, fraction of a bit time 



Fig. 12. Comparison of the initial estimate and optimum value of bit Fig. 13. Comparison of telemetry degradation as a function of 

synchronizer timing displacement as a function of carrier frequency carrier frequency for the initial estimate and optimum value of bit 
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RFI filter 


125 




TDA Progress Report 42-77 


January- March 1984 


Time Interval Errors of a Flicker-Noise Generator 

C. A. Greenhall 

Communications Systems Research Section 


Time interval error (TIE) is the error of a clock at time t after it has been synchronized 
and syntonized at time zero. Previous simulations of Flicker FM noise have yielded a 
mean-square TIE proportional to t 2 . This study shows that the order of growth is actually 
t 2 log t, explains the earlier t 2 result, and gives a modified version of the Barnes-Jarvis 
simulation algorithm. 


I. Introduction 

Let the time base of a dock be an oscillator whose output 
is sin(27 iv Q t + <p(t)), where v 0 is a nominal frequency and 4>(t) 
the phase. Then we call x m (t ) = <p(t)l(2nv 0 ) the measured raw 
time deviation of the clock. Its derivative y(t) = dx m !dt is the 
instantaneous fractional frequency deviation. The r- average 
fractional frequency deviation is defined by 


y(tj) = }f 

•'t-T 


y(s)ds 


x (t)-x (t-T) 

m 7 m v 7 


0) 


in which E is mathematical expectation and t is arbitrary. 
Thus, the two-sample variance gives the mean-square change in 
r-average frequency during a time interval of length r. 

On the other hand, the time interval error (TIE), discussed 
by Kartaschoff (Refs. 1,2), directly measures the error of the 
clock after it has been synchronized and syntonized. Let the 
clock be calibrated at time zero, that is, we assume that x m (0) 
and Y 0 are available, where Y Q is some estimate of initial fre- 
quency. We could take Y 0 =^(0) if the latter were observable, 
or Y 0 = 37 ( 0 , 7 ) for some t. The TIE at a later time t is defined 
by 


The standard measure of frequency stability is the two-sample 
(Allan) variance o 2 (t), defined by 

l m 

a y( T ) = lim 2 (m-\) ly(F,r) ~ y((r 1 )t,t)\ 2 (2) 
m-+°° ’ i= 2 

if the limit exists. For most models of clock noise, the limit 
does exist and can be computed by an ensemble average, so 
that 


x(t) = x m (t)-x m ( 0)- Y Q t (3) 

Since y = dx m /dt, we have 

x(t) = f (y(s) - Y 0 )ds (4) 

The ensemble RMS average of x(t) is called o x (t) (Ref. 1); 
thus, by definition, 


o 2 y (T) = z^E \y(t,f) - y(t-T,r )\ 2 


o 2 x (t) = Ex 2 (t) 


(5) 
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One must keep in mind that o£(f) depends on the choice of 
Y 0 . Also, it is stipulated that Y 0 must not vary with t. 

Kartaschoff (Ref. 1) investigated the behavior of TIE for 
the three random frequency-noise models whose names and 
properties are given below: 


Name 

Sy(f) 

a y( r ) 



h n 

White FM 

K 

0 
2 r 

Flicker FM 

h_Jf 

h_ x In 4 

Random Walk FM 

hjf 2 

, 27 r 2 

h - 2 —T 


where S y (f) is the one-sided spectral density of the process 
y, and h a is a constant. Having conjectured the approximation 


o\(t) = C t 2 oj(0 (6) 

for some constant C, Kartaschoff verified Eq. (6) with C = 1 
by a discrete-time computer simulation of the three FM 
processes. 

On the other hand, a recently-developed structure theory 
for clock noise processes (Ref. 3) yields the following exact 
formulas, whose derivation is given in the Appendix: 

White FM, Y 0 = 0. 

(7) 

Flicker FM, Y 0 = y(0,r 1 ). 

(i 4) ' n 7 l + ( 1+ ?;) 1 "( 1 4) 

( 8 ) 

= h_ x t 2 In 1 +0^yj asf/Tj -►«> (9) 

Random Walk FM, Y Q = y(0). 

c 2 x (t ) = h_ 2 t 3 ( 10 ) 


Thus, provided that the calibration frequency Y Q is properly 
defined, Eq. (6) holds exactly for White and Random Walk FM 
with: C= 1. For Flicker FM, however, the theoretical result 
Eq. (8) differs from Kartaschoff’s simulation result 

o 2 J,t) = (h_ l In A)t 2 (11) 

by a factor that grows logarithmically with time. 

Our purpose here is to explain the discrepancy between 
Eqs. (8) and (11). Kartaschoff used a flicker-noise generaiion 
algorithm of Barnes and Jarvis (Ref. 4) for his simulations. A 
study of this generator, which is a certain linear filter acting on 
white noise, shows that the discrepancy is caused, not by any 
design defect of the filter, nor by any defect of theory, but by 
the procedure for initializing the generator. In our results, the 
initialization used by Barnes-Jarvis and Kartaschoff (setting all 
variables to zero) yields Eq. (11) (although with a different 
constant factor), while a more complex procedure given in 
Section IV (making the output process stationary) yields an 
asymptotic result like Eq. (9), with e replaced by a different 
constant. 

In contrast, the initialization has practically no effect on 
the observed two-sample variance. No matter which of the 
two initializations is used, one observes 

o 2 (t) = h In 4 

for r greater than three times the sample time. 

II. The Barnes-Jarvis Generator 

The discrete-time n-stage Barnes-Jarvis flicker-noise genera- 
tor is a cascade of n first-order filters with transfer functions 

z- 0-37,) 

G/iz) =~ (1 _ 7 -y " , / = 1,2 02) 

where 

1 

y. = 

1 6-9 m 

For n - 5, the frequency response \H(e iu> ) I 2 of the overall 
transfer function 

H(z) = G i (z)...G n (z) (13) 

differs from the ideal two-sided spectrum h_ l nloj by at most 
0.25 dB over four decades of frequency, where h_ l = 0.2757. 
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The difference equations that implement the generator are 
yft+ 1) = (i -y j )yft)+y j _ l {t+\) 

- (1-3 j = 1,2, . . . , n 

(14) 

where y 0 (t) is the input, consisting of “standard” (mean zero 
and variance one) white noise, and y n (t) is the output. This 
implementation differs from that of Barnes-Jarvis and 
Kartaschoff only in the overall gain and in the order of appli- 
cation of the filter modules. 

The system (14) implies thaty^ and y 0 are related by an 
nth-order difference equation whose coefficients can be 
obtained by expanding the numerator and denominator of 
H(z) [Eq. (13)] in powers of z. Let the sequence i//(0), i//(l), 
i//(2), ... be the impulse response of the overall filter//; thus 

oo 

H(z) = ^2 i p(t)z~ { , \z\>l (15) 

f=0 

Assume that the standard white noise y 0 {t) is available for all 
integer t, positive and negative. The unique stationary solution 
y{t) =y n (t) of the nth order equation is 

t 

y(t) = ^2 ^^ 0 ^’ (16) 

S = — 00 

in which y{t) depends on 7 0 (s) for -°°< s < t. In the terminol- 
ogy of Box and Jenkins (Ref. 5), y is an ARMA (n,n) process. 
For a practical simulation, it can be assumed thatj> 0 (s) is avail- 
able only for s > 1 . Thus, it is natural to split y{t) into two 
independent parts, 

y(t) = yjfi+yjf) (17) 

where 

t 

y+if) = J2 Hts)y 0 (s), t > l 

= 0, t<0 (18) 

is called the present part of y(t) because it depends only on 
the random shocks y 0 after time zero, while 


o 

yjt) = i p(t~s)y 0 (s), t > 1 

S—~ °° 

= y(t), <<0 (19) 

is called the past part of y(t) because it depends only on the 
random shocks before time one. In another terminology, 
y_(t) is just the mean-square best linear prediction of y(t) 
from the past of y 0 , while y + (t ) is the prediction error. 
The stationary process y(t) will also be called the complete 
Barnes-Jarvis (BJ) process. 

The time interval error x{t) of the discrete-time frequency 
process y{t), defined by 

t 

x(t') = ^2 ~ • t>X 

5=1 

jc(0) = 0 ( 2 °) 

also breaks up into present and past parts x + (t ) and x_(t), 
obtained from Eq. (20) with_y replaced by^ + and y_. 

The BJ generator, Eq. (14), is initialized by assigning values 
to jj(0), . . . , J„(0). The usual way of doing this, called the 
zero initialization here, is to set all these values, plus_y o (0), to 
zero (Ref. 4). By an induction on n, one can prove the follow- 
ing assertion: If the zero initialization is used, then the output 
y n ( tj of the BJ generator is exactly y + ( t) for f > / . Moreover, 
since ,y + (0) = 0, x + {t) is exactly the TIE simulated by 
Kartaschoff. 

Two questions arise. First, are there any practical differ- 
ences between y{t) andy + (t)? Second, inasmuch as the infinite 
past of y 0 (t) is unavailable, can the complete process y(t) be 
simulated accurately? Since, for t > 1 , both y(t) and y + (t) 
(alias y n (t )) satisfy Eq. (14) with identical inputs y 0 (f), the 
difference y(t) - y + (t) is merely a transient of the filter H\ 
consequently, 

y(t) - y + (t) -»■ 0 as t -^°° (21) 

Nevertheless, the stationary output y(t) behaves like true 
nonstationary 1/f noise only over a time span of the order of 
magnitude 6*9"~ 1 , the longest time constant found in H. 
Over the useful time span of the generator, one cannot expect 
the transient to be small. As we shall soon see, its effect causes 
the variances of x{t) and x + (t) to differ greatly. 

The second question also has a yes answer. To simulate the 
complete BJ process y{t), one has only to set up the initial 
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vector y o (0), >^(0), . . . , J„(0) with random values having the 
correct joint distribution. The difference equations [Eq. (14)] 
do the rest. The stationary initialization algorithm is given in 
Section IV. 

III. Statistical Results 

The mean-square time interval error o 2 (t) and the two- 
sample variance Oy(f) were calculated theoretically for the 
complete five-stage Barnes-Jarvis process y(t) and for its 
present part j> + (r). Then, the same quantities were estimated 
by computer simulation, the zero initialization being used to 
generate jy + (f) and the stationary initialization (Section IV) 
being used to generate y(t). Gaussian pseudorandom numbers 
were used for the white input noise y 0 (t). The results are 
given in Figs. 1 and 2. Figure 1 plots o 2 (t)/t 2 (linear scale) 
against f/xj (logarithmic scale), where r 1 is the sample time of 
the simulation. The mean-square TIE of the complete process 
dominates that of the present part. To put it another way, 
as time goes on, more and more of the TIE (on the average) 
comes from the remote past. The results can be summarized 
by empirical formulas. For the complete process, 

o 2 x (t) = h_ t t 2 In (5.5 f/Tj) (22) 

and for the present part, 

o 2 x (t) = 2h_ 1 t 2 (23) 

valid for 16 < f/Xj < 16000. Equations (22) and (9) have 
different constants inside the logarithm. This discrepancy is 
caused by high-frequency spectral deviation of the Barnes- 
Jarvis process from pure sampled flicker FM. Equation (23) 
agrees with Kartaschoff’s result Eq. (11), except for the con- 
stant factor. This discrepancy is perhaps caused by statistical 
variations, since Kartaschoff generated only 100 sample 
functions. 

The mean-square TIE o 2 (t) of the present part can also be 
interpreted as the minimal mean-square prediction error of the 
complete x(t ) from the history of the process before time one. 
Our result, Eq. (23), is in approximate agreement with a result 
of Percival for this prediction error variance (Fig. 3.1 in 
Ref. 6) although our factor 2 h_ 1 appears to differ slightly 
from his. 

Figure 2 plots aj(f) (not o y (T )) on a linear scale against 
x/xj on a logarithmic scale. Here, the situation is quite differ- 
ent. It is apparent that the two-sample variances of the com- 
plete BJ process and of its present part are practically indistin- 
guishable. In fact, the past part of the BJ process, which 


accounts for most of the TIE for large t, contributes only 
9 percent of the complete o^(4096xj) and 0.00015 percent of 
the complete o 2 (t j). This explains why earlier BJ users, who 
relied on a constant o 2 (t) as an indicator of the success of the - 
simulations, were not aware that a large part y_ of the process 
y was missing. 

The error bars in Fig. 2 concern a side issue addressed also 
by Kartaschoff (Ref. 1), namely, the calculation of the vari- 
ance of the classical o 2 (t) estimator S 2 (T,m), which is just the 
right side of Eq. (2) without the limit. The solid bars are the ± 
one-standard-deviation error bars of S 2 (r,m) as computed for 
pure Gaussian Flicker FM by Yoshimura (Ref. 7), and later by 
the author (Ref. 3); the dashed bars come from the sample 
variances observed during the 2048 simulation runs. The 
agreement is satisfactory. 


IV. Simulation of the Complete BJ Process 

It is desired to simulate the stationary solution y(t) = y n (t ) 
of the Barnes-Jarvis system, Eq. (14). To this end, consider the 
stationary n-vector process Z(t) = ( Z x (t ), . . . , Z n (t)), -°°< t 
<°°, where Zfif) = yft) - y / _ l (t). (The Z ; - are used instead of 
the yj because the Zj are less highly correlated.) The /th 
component Zj is obtained by acting on y Q with the filter 
whose transfer function is 

Kfz) = Gj (z) . . . G._ j (z) [Gfz) - 1 ] (24) 

[Recall Eqs. (12), (13).] Each Zfit) is orthogonal to y 0 (t), and 
the covariance matrix R z is given by 

R z (ij) = E\Z.(f)Zft)\ = f KfzWfl I z)~~ z 

|zl = l 

(25) 

which can easily be evaluated by residues. 

Suppose now that we create new random variables 7 o (0), 
Zj (0), . . . , Z„(0) with the above covariances, set yj( 0) = 
J7_](0) + z /(0), 7=1, ... , n, and generate y„(t), t> 1, by 
Eq. (14) from the white noise inputs _y 0 (f), 1. One can 

then prove that these new y n (t) have the same covariance 
structure as the original stationary process. 

To set up the Z ; (0), one can use the Choleski factorization 
of R z , namely, R z = LL T , where L is an nXn lower-triangular 
matrix. The coefficients Ly are given in Table 1 for n = 6; for 
smaller n, one simply truncates the matrix. These coefficients 
give the Z ; ( 0) as linear combinations of a standard white noise 
vector ttj If the simulated process is to be Gaussian, 
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then it is important to make the Uj Gaussian because there are 
so few of them and their effect is large. Although one could 
then save time by using uniformly distributed noise inputs 
y 0 (X) (with mean zero, variance one), our simulations used 
Gaussian noise inputs throughout. 

Here, then, is the simulation algorithm for the complete 
Barnes-Jarvis process. Let GRAN be the routine that generates 
independent Gaussian random numbers with mean zero, vari- 
ance one. 

1 . Generate j> 0 (0), u 1 , ... ,u n with GRAN 

2. For i = 1 to n 

i 

LetZ. = y' L.u. from Table 1) 

/= i 

Let y ,(0) = y._ t (0) + Z. 

Next i 

3. For t = 0, 1, 2, . . . 

Generate ^ Q (t+1) with GRAN 

For / = 1 ion 


Compute j>.(f+l) by Eq. (14) 

Next / 

Next t 

Remark: To generate the present part of the BJ process, set 
Z o (°)>Zi(0), ■ . . ,y n ( 0) = 0, then go to Step 3. 


V. Concluding Remarks 

As we have seen, the Barnes-Jarvis generator, if given the 
stationary initialization, produces an output whose mean- 
square TIE a^{t) has a t 2 log t behavior, just as theory pre- 
dicts. The previously-observed t 2 behavior, which appears 
when the BJ generator is given the zero initialization, can be 
interpreted as the minimal mean-square error of a linear clock- 
time predictor based on knowledge of the clock’s behavior 
over the entire past, to 0. Naturally, this error is smaller 
than the mean-square error of a predictor x m (0) + Y 0 t 
for the time deviation x m (t). We conclude that if the diagnosis 
of the noise of a clock includes a Flicker FM component, then 
its al(t) must include a t 2 log t component. Ify(0,Tj) is used 
for calibrating the initial frequency, then this component is 
given by Eq. (8). 
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Table 1. Coefficients L, / for the stationary initialization of the Bames-Jarvis process 


i 



i 




1 

2 

3 

4 

5 

6 

i 

0.603023 






2 

0.214635 

0.512223 





3 

0.301626E-1 

0.241088 

0.494406 




4 

0.345089E-2 

0.358003E-1 

0.244953 

0.491688 



5 

0.384698E-3 

0.412554E-2 

0.366905E-1 

0.245520 

0.491287 


6 

0.427600E-4 

0.460283E-3 

0.423277E-2 

0.368209E-1 

0.245599 

0.491231 






Appendix 


Derivation of Formulas for Mean-Square Tie 


Given fixed numbers t, r > 0, define the second-order dif- 
ference operator L by its action on a function /(s): 

Lfis) = f(s+t)~ (l +^j f(s) +~f(s-r) (A-l) 

where s is used as the time variable because t is fixed. The 
transfer function of L is 

L(ioS) = e i0Jt - (l +^)+~e~ iwT (A-2) 

Then, with Y 0 = y( 0,t), Eq. (3) can be written as 

*0) = LxJ 0) (A-3) 

Now assume that the second differences of the time devia- 
tion x m are stationary. (The examples treated below satisfy 
this condition.) According to the structure theorem of Ref. 3, 
there is a function D(t ) (appearing as 2 Re C(t ) in Ref. 3) such 
that 

o 2 x (t) = Ex 2 (t ) = E Lx JO) Lx JO) = AD(0) 

(A-4) 

where A is the fourth-order operator LL*, andi* is the opera- 
tor with transfer function L*(iu>) = L(icj)*. Hence the transfer 
function of A is 

A(/co) = \L(ico)\ 2 = 2 |l -J + |l + t ^je iut + 



+ l| e 'w(f+T) + £-ioj(f+r)j (A-5) 

Thus, from Eq. (A-4), 

a 2 x (t) = 2 |l + ~ +~|/J(0) _ ( 1 +^)(Z)(/) +£»(-?)) 

- - t (l +^)(d(t) + D(-r)) + ^D(t+T) + A-i-r)) 

(A-6) 


Each noise process has its own .©-function to be used in 
Eq. (A-6). 

I. Flicker FM 

From Eq. (27) of Ref. 3, 

h . 

D(t) = -y t 2 In |f I (A-7) 

which, when substituted into Eq. (A-6), gives Eq. (8) with 


II. Random Walk FM 

Since A is a fourth-order operator, one may add a third- 
degree polynomial to D. Thus, Eq. (28) of Ref. 3 is equivalent 
to 

h it 2 

D(t) = \t\ 3 (A-8) 

and Eq. (A-6) gives 

o 2 x (t) = h_ 2 y- t 2 (t+r) (A-9) 

As t -» 0,y(0,T) -*>>(0) because y is a mean-continuous process 
(Brownian motion), and Eq. (A-9) reduces to Eq. (10). This 
does not happen for Flicker FM because it has too much 
power at high frequencies. Consequently, y(0,r) is well de- 
fined, but y(0) is not. 

III. White FM 

Since Y 0 = 0, and x m has stationary first differences, we rede- 
fine L as the first-order operator : 

Lfis) = f(s+t) - f(s) (A- 10) 

Then everything works as before, except that A is now a 
second-order operator. In place of Eq. (A-6), we have 
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o 2 x (t) = 2Z>(0) - D(t) - D(-t) 


(A- 1 1) 


(A- 12) 


K 

D(t ) = ~j W 

This time, a first-degree polynomial may be added to D(t). which combines with Eq. (A- 11) to give Eq. (7), itself a well- 
Thus, Eq. (26) of Ref. 3 is equivalent to known fact. 
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This article considers codes for use in personal computer file transfer as control char- 
acters, when only upper-case ASCII can be used to avoid dependence on unique machine 
features and promote portability. If ten control functions are needed, a number used in 
at least one protocol, we seek a subset of ten upper-case ASCII characters with good dis- 
tance properties. The. control functions form themselves naturally into three groups, one 
of two functions (ACK and NAK) and two of four. We wish to make ACK and NAK as 
antipodal as possible (distance 6), make the distances within each of the other groups as 
large as possible (4), and otherwise have as few 2 ’s in the distance table as possible, recog- 
nizing that only even distances can occur. We find the minimum and an assignment that 
attains the minimum. The code is essentially unique. We also solve the analogous problem 
for two groups of three control functions and one group of four. 


I. Introduction 

Here we solve a simple coding problem in computer proto- 
col design. People with personal computers want to exchange 
binary files, but everyone owns a different machine. The 
problem is that different machines interpret binary seven- 
tuples differently. So when setting up communication, it may 
be impossible to even pass the initial protocol characters to 
the protocol-handling software. But there is a way out of this. 
All personal computers agree on what to do with upper-case 
ASCII. So in this article we restrict protocol control functions 
to upper-case ASCII. This permits communication to be estab- 
lished between the two computers, after which the files can 
then be transferred as straight binary files. 


Currently a graduate student at Harvard University 


Thus, we are in the following situation. We want to choose 
upper-case ASCII subsets with good distance properties. 
Before attempting this, we will note that ASCII here is ac- 
tually 7 bits with even parity adjoined. However, two of the 
8 bits are always equal in the upper-case alphabet. This means 
that we really have a 6-bit code to choose (see Table 1). 

Let us formulate this as a precise coding problem. A typical 
protocol has 10 “block control functions.” Here we choose to 
divide them into three groups, as in Table 2, a “2, 4, 4, Table.” 

The rationale is that functions within groups need extra 
protection from each other. This is because within groups, 
the control functions are more likely to be confusable because 
the protocol states typically result in outputs that can be or 
tend to be in only one group. Furthermore, ACK and NAK 
need to be as unconfusable as possible, to prevent false file 
transfers. 
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Other groupings are possible, for example a “3, 3, 4 Table” 
with NAK, ACK, QRY in Group 1, WT, ABH, ABW in 
Group 2, and EOD, EOT, SOH, SOD in Group 3. We shall 
consider this case as well. 

Let us now formulate this as a coding problem. We are to 
choose a subset of 10 of the 26 upper-case ASCII odd-parity 
six-tuples (odd parity, because the invariable 10 has been 
removed), with the following properties: 

a) The distance of ACK and NAK is 6 . 

b) The minimum distance within Groups 1 and 2 should 
be as large as possible. 

c) There should be as few distances of 2 as possible. 

This is for the 2, 4, 4 case. For the 3, 3, 4 case, we want the 
distances within groups to all be at least 4 (this makes them all 
4, it turns out), with as few distances of 2 as possible. We will 
call a set of three code words of mutual distance 4 an equi- 
lateral triangle and of four code words of mutual distance 4 a 
regular tetrahedron. We note here that for both problems, it 
turns out that the restriction to the 26 upper-case ASCII out 
of the 32 odd-parity six-tuples did not hurt the distance table 
any, as it turned out in these two cases. Thus we will talk 
about the alphabet restriction no longer. Probabilistic argu- 
ments can be given that make this not too surprising. 


d(oi.,a.) > 4, 

d(M,) > 4 > 

d{l v lj) > 4, 

and the number of 2’s in the distance table of oi l , a 2 ,a 3 ; j3j', 
0 2 > h ; > 7 2 » 73 , 74 is minimal. 

Five propositions will be useful in our search. 

Proposition 1: Let a l , a 2 , a 3 e E and d(a j , a 2 ), d(a 2 , a 3 ), 
d(a j , a 3 )>4. Thend(aj , a 2 ) = d(a 2 ,a 3 ) = d(a j , a 3 ) = 4. 

Proof: Assume the contrary. Say d(a j , a 2 ) > 4, i.e., = 6. 
Then a' 2 = a'j where x means 1 - x for x = 0, 1. Thus, for 
any a e {0, 1 }, \a - a‘ 2 \ = la - a'j | (i.e., if one of them is 0, 
then the other is one and vice versa) and, hence, \a - a 2 | + 
\a - a'j | = 1. Then 

6 

8 <:d(a 1 , a 3 ) + c?(a 2 , a 3 ) = ^ (la' 3 - « 2 I + |a' 3 - aj |) 

1=1 

6 

-E>-« 

t=i 


i, j 1, 2, 3, if 1 } \ 

U = 1,2,3, i ¥= j \ (2) 

i, j — 1,2, 3,4, i f- j J 


II. Intermediate Results 

Here we will for convenience revert to even-parity six-tuples 
instead of odd parity. Call this set E. We shall use lower-case 
Greek for the elements of E. Note that the Hamming distance 
d between any two elements of E is even. This section presents 
five propositions needed in deriving the optimal codes. 


This contradiction proves the proposition. 

Proposition 1 shows that we can replace “> 4” by “= 4” 
everywhere in Problems (1) and (2). 

Proposition 2: Let A = { 7 j , 7 2 , 7 3 , 7 4 } be a regular tetra- 
hedron. Then each column of the matrix 


We want to solve the following two problems: 

Problem 1: Find «j, a 2 ; |3 1 , 0 2 , 0 3 , 0 4 ; 7 j , y 2 , y 3 
such that 

d(a v a 2 ) = 6 

d(0 r &j) >4, i, j = 1,2, 3, 4, i =^/ 

d{y v y j ) > 4, i, j = 1, 2, 3, 4, i =£/ 



M) 

m 

~(2) 
' 1 

,y( 6 ) 

'1 

7 4 e E 
\ 

~(D 

'2 


y(6) 

'2 

i" 

~(D 

'4 


y(6) 

' 4 


and the number of 2 ’s in the distance table of Oj, a 2 ;j3 1 ,0 2 > 
0 3 , 04 ; 7i , y 2 » 7 3 , 7 4 is minimal. 


has exactly two 0 ’s and two l’s in it, and is unique up to per- 
muting and conjugating (complementing modulo 2 ) columns. 


Problem 2: Find , a 2 , a 3 ; 0j , jS 2 , 0 3 ; y x , y 2 , y 3 , y 4 e E 
such that 


Proof: We can assume without loss of generality that y x = 
( 0 , 0 , 0 , 0 , 0 , 0 ); otherwise just conjugate the columns where 
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= 1, i.e., replace 0 by 1 and 1 by 0 in those columns. 
Then y 2 , y 3 , and y 4 have exactly two 0’s and four l’s. By 
permuting the columns we can make y 2 = (1, 1, 1, 1, 0, 0). 
Now y 3 and 7 4 must have their last two positions be 11. 
Otherwise their distance from y 2 would be at most 2. By 
permuting the first four columns we can make y 3 = 
(1, 1,0,0, 1, 1). Since d( y 3 , 7 4 ) = 4, 7 4 must be equal to 
(0,0, 1, 1, 1, 1). 

Thus any regular tetrahedron can be obtained from 

(0, 0, 0, 0, 0, 0) = 7j 

(1,1, 1,1, 0,0) = y 2 
(1, 1,0,0, 1, 1) = 7 3 
( 0 , 0 , 1 , 1 , 1 , 1 ) = 7 4 
by permuting and conjugating columns. 

Since these operations preserve the property we are inter- 
ested in (exactly 2 zeros in every column), it is sufficient to 
check it just for the regular tetrahedron above. It, indeed, has 
exactly two 0’s in each column. Hence, any regular tetrahe- 
dron also has this property. (Note that we can only conjugate 
an even number of columns if the tetrahedron is to remain a 
subset of E.) 

Proposition 2 also follows from a more general result on 
constant-distance codes, but we shall not do it this way. 


(0, 0, 0, 0, 0, 0) = 7 j 

( 1 , 1 , 1 , 1 , 0 , 0 ) = y 2 

(1, 1,0,0, 1, 1) = 7 3 

as we showed in the proof of Proposition 2. It is easy to see 
that 7j , y 2 , and 7 3 have the desired property. 


Proposition 4: Let C = {7 x , 7 2 , 7 3 , 7 4 } be a regular tetra- 
hedron of elements in E, and let 7 also be in E. Then 

4 

Y d {y,y.) = 12 
i=l 


Proof ■ Proposition 2 implies that for a e {0, 1}, k - 
1, 2 , . . . , 6, we have 

Y = 2 

i=l 

for each column k. Thus 
4 

Y l7 (fc) - y ( i k) \ = 2, 1 <£<6 

;=1 

Adding these six equalities together, we get, as stated, 


Proposition 5: Let B = {/3j , P 2 , P 3 
angle. Then the matrix 

} be an equilateral tri- 

11 

04 

ft It 18 
\ 1=1 

0 ( , 2) 


4 

-El 

[t lr<‘> 


4 6) 

i-1 

\ k=l 

0 ( s° 

4 6) 

4 

= Y ^>7, ) 
/= 1 


has either one 0 and two 1 ’s or one 1 and two 0’s in each col- 
umn, and is unique up to permuting and conjugating columns. 


Proof: We apply the same argument as in Proposition 2. 
The property we are interested in is again invariant under con- 
jugations and permutations of columns. This operation allows 
us to transform any equilateral triangle into 


7j 


7/ 


Propositions. Let A = (ctj , a 2 , a 3 },B = { /3j , (3 2 , 0 3 }be 
two equilateral triangles of elements of E. Then 


Y d(a v 0 ,) < 30 

/,/=!, 2,3 
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If the sum is 30, then each 


Thus 


E | a ( fc >-0( fc >| = 5 

!’,/= t ,2 ,3 


/,/'= 1,2,3 


for 1 < k < 6 (used for exhaustive search). 
Proof '. First we show that 


for each k from 1 to 6. Adding these six inequalities together, 
we get, as stated, 


E | a W-j3j fc )|<5 

/,/= 1 ,2,3 


6 


30 

fc=i 


<',/=!, 2,3 


for each k - 1, 2, . . . , 6. The sum 


£ 

'= 1 ' 



does not change when we simultaneously conjugate all the kVa 
components 


= E 

i,j~ 1 ,2 ,3 \ fc=l / 

= E^M/) 


The second part of the proposition also follows. 


and permute 
or 

for each k= 1,2, ... ,6. 

Thus by Proposition 3 we can assume that a^, = 0 

and a ^ = 1. Then (again by Proposition 3) one of the j3^ 
(; = 1, 2, 3) is 0 and one of them is 1. Number the 0 first 
(i.e., 0^ = 0) and the 1 last (/3^ = 0). Now we only have two 
choices: = 0 and j3^ = 1. In the first case 

E i^-^i = 4 

/,/= t ,2,3 

and in the second case 

/,/= 1,2,3 


III. The Optimal Codes 

Figure 1 gives the optimal odd-parity 8-bit ASCII upper- 
case alphabetic code for the 3, 3, 4 case (2) and Fig. 2 an 
optimal code for the 2, 4, 4 case (1). Why are they optimal, 
even when we drop the alphabetic restriction and would be 
willing to allow any even-parity six tuples? 

Suppose we have a distance table satisfying (2), with 
equalities instead of inequalities by Proposition 1 . Break it up 
into blocks as shown in Fig. 3. By Proposition 4 each row sum 
within Blocks II and III must be 12. Since every entry is 2, 4, 
or 6, this means that each row in Blocks II and III contains at 
least two 2’s. Thus the total number of 2’s in Block II is at 
least 6 and the total number of 2’s in Block III is at least 6. 

By Proposition 5 the sum of the nine entries in Block I is 
at most 30. This implies that Block I contains at least three 
2’s (7 X 4 + 2 X 2 = 32 > 30). Blocks T, Il\ and III' are just 
transposes of Blocks I, II, and III, respectively, and, hence, 
always have the same number of 2’s as Blocks I, II, and III. 
We see that Fig. 1 has the minimal possible number of 2’s in 
each block. This shows that Fig. 1 is a solution for Problem 2. 
There are thirty 2’s in this table. 

The same argument (see Fig. 4) shows that any distance 
table satisfying (1) must have at least four 2’s in Block I, 
four 2’s in Block II, and eight 2’s in Block III (7 X 2 + 9 X 4 = 
50 > 48). This proves that Fig. 2 has the minimal possible 
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number of 2’s in each block and is a solution for Problem 1. 
There are thirty-two 2’s in this table. 

If 20 control functions are desired instead of 10, we can 
use pairs of upper-case letters. We have found a constant dis- 
tance 6 code that does this for 20 functions. This will be 
reported elsewhere in a more general context. 


IV. Uniqueness Results 

The code of Fig. 2 with its three groups has a unique auto- 
morphism crossing group boundaries. In this, O, H exchange 
with K, L (or Q, V with U, R), with O corresponding to L and 
H corresponding to K. Similarly, in Fig. 1, exchange X and Z 
between groups. This forces V to interchange with 0 and L 
with K. Also, in each problem, the two groups of equal size 
can be swapped. 


Are our solutions unique up to the obvious operations? The 
answer is yes for the 2, 4, 4 problem. This is easier to prove 
than to write up, and we will merely assert it. 

The answer is no for the 3, 3, 4 problem. In fact, a differ- 
ent distance table can even be achieved, although it can be 
shown that the two triangles are unique up to obvious trans- 
formations. Merely interchange U and A in Fig. 1, and observe 
that Block II (upper right) becomes the inequivalent (because 
of the column of all 4’s) block 

2 4 2 4 

4 2 2 4 

2 2 4 4 

This is also inequivalent to Block III, so it is a really different 
distance table. We shall say no more about these uniqueness 
problems. 
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Table 1. Upper case ASCII with even parity 


A 

0100 

0001 

B 

0100 

0010 

C 

1100 

0011 

D 

0100 

0100 

E 

1100 

0101 

F 

1100 

0110 

G 

0100 

0111 

H 

0100 

1000 

I 

1100 

1001 

J 

1100 

1010 

K 

0100 

1011 

L 

1100 

1100 

M 

0100 

1101 

N 

0100 

1110 

O 

1100 

mi 

P 

0101 

0000 

Q 

1101 

0001 

R 

1101 

0010 

S 

0101 

0011 

T 

1101 

0100 

U 

0101 

0101 

V 

0101 

0110 

w 

1101 

0111 

X 

1101 

1000 

Y 

0101 

1001 

Z 

0101 

1010 


Table 2. Three groups of block control functions in the 2,4,4 case 


Function 

Meaning 


Group 0 

ACK 

Acknowledgement 

NAK 

Negative acknowledgment 


Group 1 

EOD 

End of data 

EOT 

End of text 

SOH 

Start of header 

SOD 

Start of data 


Group 2 

WT 

Wait 

QRY 

Query 

ABH 

Abort and hang up 

ABW 

Abort and wait 



FUNCTION 


CHARACTER 


8-BIT ASCII 


F 


GROUP 0. 

ACK 

F 

01000110 


nak 

...... X 

01011000 


WT 

A 

11000001 

GROUP 1 , 

.QRY 

O 

01001111 

\ " 

ABH 

Z 

11011010 


ABW 

U 

11010101 

GROUP 2 

EOD. 

V 

11010110 

\ 

EOT 

Q 

01010001 


. _,_$OH 

L 

01001100 


SOD 

K 

11001011 




DISTANCE MATRIX 

GROUP 0 

GROUP 1 

GROUP 2 


ACK (F) 

0 4 4 

2 4 4 

2 4 2 4 

GROUP 0 

NAK (X) 

4 0 4 

4 2 4 

4 2 2 4 


WT (A) 

4 4 0 

4 4 2 

4 2 4 2 


QRY (O) 

2 4 4 

0 4 4 

4 4 2 2 

GROUP 1 

ABH (Z) 

4 2 4 

4 0 4 

2 4 4 2 


ABW (U) 

4 4 2 

4 4 0 

2 2 4 4 


EOD (V) 

2 4 4 

4 2 2 

0 4 4 4 


EOT (Q) 

4 2 2 

4 4 2 

4 0 4 4 


SOH (L) 

2 2 4 

2 4 4 

4 4 0 4 


SOD (K) 

4 4 2 

2 2 4 

4 4 4 0 


3 TWO's 
PER ROW 


3 TWO's 


Fig. 1. Optimal code for the 3,3,4 case 


FUNCTION CHARACTER 8-BIT ASCII 


GROUP 0 ACK E 01000101 r 6 

\ nak Z 11011010 u 




DISTANCE MATRIX 

GROUP 0 

GROUP 1 

GROUP 2 


ACK (E) 

0 6 

2 4 2 4 

4 2 2 4 


NAK (Z) 

6 0 

4 2 4 2 

2 4 4 2 


EOD (O) 

2 4 

0 4 4 4 

2 2 4 4 


EOT (H) 

4 2 

4 0 4 4 

2 2 4 4 


SOH (Q) 

2 4 

4 4 0 4 

4 4 2 2 


SOD (V) 

4 2 

4 4 4 0 

4 4 2 2 


WT (K) 

4 2 

2 2 4 4 

0 4 4 4 


QRY (L) 

2 4 

2 2 4 4 

4 0 4 4 


ABH (U) 

2 4 

4 4 2 2 

4 4 0 4 


ABW (R) 

4 2 

4 4 2 2 

4 4 4 0 


4 TWO's 3 TWO's 


4 TWO's 
PER ROW 


3 TWO's 
PER ROW 


Fig. 2. Optimal code for the 2,4,4 case 
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The relative probability of detecting randomly distributed CW transmitters as a func- 
tion of the fraction of the sky which is searched (in a fixed time) is given. It is shown that 
the probability of detecting such a class of transmitters with a given receiving system is a 
maximum if the entire sky is searched. The particular case of a search in which the num- 
ber of directions searched is equal to the telescope gain and the integration time per beam 
element is equal to the reciprocal of the channel bandwidth is discussed. 


I. Introduction 

It is important in designing a search program for SETI 
(Search for Extraterrestrial Intelligence) to understand the fac- 
tors which effect the probability for detecting signals from dis- 
tant transmitters. An optimum search strategy cannot be 
identified at the present time since the statistical properties of 
the transmitters and the motivations of the senders are not 
known. Nevertheless, a simple model, in which the transmit- 
ters are uniformly distributed in space and have a power law 
distribution of intrinsic transmitted power, appears to have 
sufficient generality to provide some guidance in designing a 
search program. We examine, in this note, the relative proba- 
bility of detecting signals from such a class of transmitters in a 
fixed period of time, as a function of the fraction of sky which 
is searched. The scenario envisaged is that the observer may 
choose between scanning slowly (and achieving high sensitivity 
in a few directions) or scanning the entire sky (at reduced 
sensitivity). 

We show that the probability of detecting such a class of 
CW transmitters is maximized if the entire sky is searched 
when the received signal-to-noise is low. This means that one is 


more likely to detect a weak CW signal by scanning the entire 
sky than by concentrating on a few areas. Since strong CW sig- 
nals can be detected either by focusing the search in a few 
areas or by scanning the entire sky, the primary interest in this 
paper is on the weak signal case. 

The detection of pulsed signals is not considered in this 
note. The detection of these signals may favor one scanning 
technique over another depending on factors such as the duty 
cycle of the transmitters. 

II. The Probability Relationship 

Derivation. Following Drake (Ref. 1), we assume that sig- 
nals from civilizations radiating an equivalent isotropic radi- 
ated power (EIRP), P a , are originating from uniformly spaced 
locations having a density of ( n Q ) transmitters per unit vol- 
ume. We ignore the fact that the galaxy is a highly flattened 
disk and consider the spherical volume in the vicinity of the 
sun where stars are distributed more or less at random except 
for a tendency to cluster. This volume has a maximum radius 
of approximately 1 kiloparsec. The maximum range of a given 
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search system is related to the minimum detectable flux, S m , 
by the expression 


R = 



( 1 ) 


and the total number of detectable signals within the spherical 
volume is given by 



The minimum detectable flux, S m , depends on the diam- 
eter, D, of the receiving antenna, the system noise tempera- 
ture, T s , and the receiver channel bandwidth, B s . Oliver and 
Billingham (Ref. 2) show that the minimum detectable (coher- 
ent) flux, defined as being equal to the rms fluctuations due 
to the noise, is equal to 


_ kT s B s 1+(1+») 1/2 
m 7T& n 

4 


(3) 


and the range at which a signal can be detected with a given 
receiving system is given by 


/d\ 

P o n 

w 

kT s B s (1+(1 +H) 1 ' 2 ) 


In these expressions 

k - Boltzmann’s constant 
T = system noise temperature 
D = receiving antenna diameter 
B = receiver channel bandwidth 

S 

t - integration time 

n = (Bt) = number of independent samples averaged 

Drake (Ref. 1) shows that the probability of success of an 
observation program is proportional to (a) the total frequency 
searched, (b) the total solid angle searched, and (c) the spheri- 
cal volume defined by the radius R. These proportionalities 
lead to the following expression for the probability of success, 
P s of a search over M different directions using a system having 
C channels (total bandwidth = CB S ), and an antenna beam 
solid angle equal to £2^: 


In this expression, K is a constant of proportionality, 

(= 47r) is the total solid angle of the sky, and B t is the total 
bandwidth range in which the signal is confined. The ratio 
(MQ, J£l t ) is the fraction of the sky which is surveyed and the 
ratio ( CBJB t ) is the fraction of the total bandwidth which is 
surveyed. 

Noting that the telescope gain, solid angle, the diameter are 
related through the equations 



and the total observation time, 0, number of sky elements 
examined, M, and the number of independent samples at each 
sky element, n, are related through the equations 


n = B t = B 

s 


sM 


(7) 


the probability of success can be written as follows 


Sl\/B 

P_ = KnMC\-r\\^- 


aj\Bl\ 64 


D : 


L kT s B s 1 + V 1 +«J 


3/2 


( 8 ) 


,3 /B 


Kn MC 

o 


(4 7l) 


4 \BJ 2 


G 1/2 


B 0 

s 


kTB I — x 

s s M+\/M 2 +MB 0 


3/2 


(9) 


In the derivation of Eq. (9), we assumed that each transmitter 
had the same intrinsic power. It is easy to generalize the equa- 
tion for the case of a continuous distribution of powers by 
rewritting Eq. (9) as a differential probability of success. 

\ 3 

dP = KMC 
S (4?r) 4 


\kTB , — = 

* s M + y/M 2 +MB s (t> 

If the distribution of powers is a power law of the form 


3/2 

P 3 / 2 P(P)dP 

( 10 ) 
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p(p ) = K^r* 


(ii) 


(13) 


M<B 0 
s 

Equation (10) may be integrated between the limits P u and P L This limit is in addition to the constraint set by the number of 

to yield the following expression for the probability of resolvable directions in the sky. If an attempt is made to in- 
success: crease M beyond this critical product, beam smearing will 

occur, and it will not be possible to differentiate spatial direc- 
tions on the scale of the antenna beam size. 

Figure 1 shows that it is advantageous to carry out a search 
by searching the entire sky (M = Gain) with the highest gain 
(largest) antenna available. This result holds provided that the 
telescope gain is less than the product B s <j>. Figure 2 shows 
how the probability of success varies with the telescope gain 
under the following assumptions: 

(12) G<B s cj> M = G (entire sky is surveyed) 

It is seen from this equation that the spectral index, a, deter- G >5^0 M = Bjj) (MSI <4 it) = 10 8 

mines whether the strong sources or the weak sources will 

dominate the probability of success. If the spectral index is Figure 2 shows that the probability of success increases as 
less than 5/2, then the strong sources, although less numerous (Gain) 3 / 4 up to the critical point, B s 4> = G; at larger values of 

will dominate. If the spectral index is greater than 5/2, then gain, the probability increases as (Gain) 1 / 2 . This figure shows 

the weaker and more numerous sources will dominate. If the that it is always better to use a large antenna rather than a 
spectral index is exactly 5/2, then all sources regardless of small antenna. Of course there are mechanical considerations 

their power will contribute uniformly to the probability of which may favor the use of a smaller antenna; these are not 

success. discussed here. 



In Eqs. (9) and (12), we have implicitly assumed that the 
telescope gain and solid angle are constant over the receiving 
bandwidth of the system (C B s ). This is true in most high per- 
formance antenna receiving systems, but may break down if 
the receiving bandwidth is very large. In this case, an integra- 
tion over the wavelength needs to be performed. 

Interpretation. Figure 1 shows the dependence of the prob- 
ability on the number of different beam areas examined, M, 
for a number of different telescopes whose gains vary from 
10 2 to 10 8 . The value of B s 0 used is 10 s . We note in this 
figure that for any given telescope, the probability increases as 
M increases up to the point where the telescope gain becomes 
equal to the number of directions examined. M cannot exceed 
the gain since the gain sets the limit on the number of indepen- 
dent areas in the sky which can be examined. 


It is of interest to make a comparison between two identi- 
cal receiving systems which both operate for the same time but 
with different amounts of the sky being examined. The ratio 
of the probability of success of a search which examines M 
different directions, to that of a search which examines the 
entire sky (M = G = gain) is given by 



Figure 3 shows this relationship for the particular case of a 
60 dB gain antenna and a channel bandwidth-total observation 
time product, B s (j) = 10 8 . The figure clearly illustrates that the 
relative probability for detecting a weak signal increases as 
more directions are searched. 


Another constraint on the number of different directions 
which can be examined is superimposed by the system channel 
bandwidth, B s . Since independent measurements cannot be 
made more frequently than the reciprocal of the bandwidth, 
the product of the channel bandwidth and total observation 
time is an upper limit on the number of independent direc- 
tions in the sky which can be examined in a specified time. 


III. The Special Case of the Matched System 

It is of interest to consider the special case in which the 
integration time per beam element is matched to the channel 
bandwidth, B s , such that B s t = 1, and the telescope gain is 
matched to the total observation time such that G = M = B s <j>. 
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These two conditions define a matched system and lead to the 
following expression for the probability: 


P = 


K 3 
(47r) 4 




\B. 


v3/2 


kT 


(15) 


It is noticed that this expression is independent of the tele- 
scope gain. This at first sight appears to be in contradiction 
with our previous result which stated that larger antennas pro- 
duce a higher probability of success than smaller antennas. 
However, since the matched condition requires that channel 
bandwidth vary with the telescope gain, the previous result 
which assumed the bandwidth to be constant does not apply. 


In order for matched systems to have the same probability 
of success, they must cover the same total bandwidth. Since 
matched low gain antenna systems utilize smaller channel 
bandwidths than higher gain antenna systems (channel band- 
width is proportional to Gain in matched systems), they must 
compensate by using a larger number of channels. The table 
below shows parameters for three matched systems with the 
same probability of success. The total observation time is 
assumed to be 10 7 seconds. The highest gain system is 
assumed to have a channel bandwidth of 1 Hz and a total 
bandwidth of 1 KHz. 


Table 1 illustrates the ultranarrow channel bandwidths and 
large number of total channels which are required for a small 
antenna to have the same probability of success as a large sys- 
tem. Since the interstellar medium imposes a minimum usable 
bandwidth (i.e., Ref. 3). the bandwidth cannot be reduced 
indefinitely. Hence, in practice, a small antenna system cannot 
be made as sensitive as a large antenna system, even if the 
number of channels required imposes no limits. 


IV. Discussion 

We have shown in this note that the probability of detect- 
ing a weak signal from randomly distributed transmitters in- 
creases faster by increasing the fraction of sky that is searched 
than by increasing the sensitivity in a given direction. The 
reason why it is more advantageous to search a larger solid 
angle at reduced sensitivity than to search a smaller solid angle 


at increased sensitivity can be understood through the follow- 
ing argument. The number of detectable transmitters increases 
as the minimum detectable flux raised to the -3/2 power 
[Eqs. (1) and (2)] . Since the minimum detectable flux de- 
creases as the inverse square root of time for weak signals, the 
number of detectable transmitters in a given direction in- 
creases as the observation time raised to the 3/4 power. On the 
other hand, the number of detectable transmitters is also pro- 
portional to the solid angle searched. Since the solid angle 
searched can be made to increase in proportion to the time by 
searching in different directions, the number of detectable 
transmitters can be made to increase in direct proportion to 
the time (rather than to the 3/4 power). Hence, the probabil- 
ity of success increases faster by scanning the entire sky than 
by concentrating in the search to a few directions. The ratio of 
the number of transmitters which could be detected with a 
given telescope system in a given time if the entire sky is 
searched to those that could be detected in the same time if a 
single direction is observed is proportional to the observation 
time raised to the 3/4 power. 

Although this result was derived for the case of coherent 
radiation, it applies more generally to detection situations in 
which the minimum detectable flux decreases more slowly 
than the observation time raised to the -2/3 power. For exam- 
ple, the minimum detectable flux varies as the inverse square 
root of time for incoherent radiation. Hence the result holds 
for incoherent as well as (weak) coherent radiation. The result 
breaks down whenever the minimum detectable flux decreases 
faster than the observation time raised to the -2/3 power. One 
example where this occurs is in the detection of strong coher- 
ent signals. In this case, the flux varies inversely with time. 

Aside from mechanical considerations, a large telescope can 
always be used to produce a higher probability of success than 
a smaller telescope. Also, it is possible that a small telescope 
which surveys the entire sky will have a higher probability of 
success than a larger telescope that concentrates a search in a 
few directions. 

The probability of success of a “matched” system is the 
same for all telescopes regardless of their gain provided that 
they cover the same bandwidth. This requires that small tele- 
scopes use narrow bandwidths and thus a larger number of 
channels. 
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Table 1. Channel bandwidths and total of channels required for a 
small antenna to have the same probability of success as a large 
system 


Gain 

Observation Time 
per Sky Element 

Number of 
Elements 

Bandwidth 

Number of 
Channels 

10 6 

10 

10 6 

1 

10 3 

10 5 

100 

10 s 

0.1 

10 4 

1 

10 7 

1 

10~ 6 

10 9 
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Fig. 1. Relative probability of success for a number of different 
antennas whose gains vary from 10 2 to 10 8 as a function of the 
number of directions in the sky which are searched. The search 
time-channel bandwidth product is taken to 10 s . 



Fig. 2. Relative probability of success as a function of the antenna 
gain. The search time-channel bandwidth product is taken to be 
10 s . 



Fig. 3. Relative probability of success as a function of the number of 
directions in the sky which are searched for a 60 dB antenna. Curve 
is normalized to unit probability at 10 6 directions which corre- 
sponds to the entire sky being searched. The search time-channel 
bandwidth product is taken to be 10 s . 
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This article presents an analysis of the scalloping problem for the case of Gaussian 
noise statistics. We derive the optimal weighting strategy for linearly combining two ob- 
servations in adjacent beam areas, and compare the sensitivity and scalloping for this 
weighting strategy with that realized using a single observation or using equal weighting of 
two observations. We also calculate the variation of the probability for detecting ETI 
signals with scan separation for the various weighting strategies, assuming that the trans- 
mitters are of equal strength and are uniformly distributed throughout space. 


I. Introduction 

One component of the Search for Extraterrestrial Intelli- 
gence (SETI) will be a survey of the entire celestial sphere over 
a broad frequency range to a significantly low limiting flux 
(Ref. 1). The inherent advantage of this strategy is that all di- 
rections are observed, and thus any signal which exceeds the 
threshold of the search will be located. The SETI program 
must design an efficient sky survey strategy which realizes 
the survey goals within the constraints of available antenna 
resources. 

Time constraints and antenna dynamics dictate that the 
survey will be carried out by smoothly sweeping the beam 
across the sky. At the end of each scan, the motion must be 
reversed without exciting the antenna’s natural modes of 
mechanical oscillation. At the same time the pointing in the 
orthogonal direction must be stepped so that the subsequent 
scan traverses a neighboring strip of the celestial sphere. Over 
a period of time, all of the celestial sphere available to the 
antenna will have been observed. Unless care is taken in the 


design of the survey, however, the sensitivity to an ETI signal 
will be a periodic function of position with respect to the scan 
pattern. This feature is commonly referred to as scalloping and 
its magnitude is defined to be the quantity of minimum sensi- 
tivity minus maximum sensitivity. 

This article is one of a series of technical reports proceeding 
from the SETI sky survey definition studies and presents an 
analysis of the scalloping problem for the case of Gaussian 
noise statistics and ETI sources whose signal strength does not 
change with time. We derive an optimal weighting strategy for 
linearly combining data in two adjacent beam areas and then 
compare the sensitivity and scalloping achieved using this 
strategy with that realized utilizing practical weighting strate- 
gies. Finally, we derive the variation with scan separation of 
the probability for detecting ETI signals, assuming that the 
transmitters are of equal strength and uniformly distributed 
throughout space. Subsequent papers will extend our analyses 
to the case of non-Gaussian noise statistics, the realm in which 
the contemplated sky survey will be operating due to the small 
number of independent samples comprising an accumulation. 
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II. Derivation of the Optimal 
Weighting Strategy 

In this section, we derive an optimal strategy for linearly 
combining two observations to achieve the maximum sensitiv- 
ity to a signal. To simplify the calculation, we assume a one 
dimensional model (see Fig. 1). The x-axis represents the loci 
of possible source positions between scan tracks, and we 
search by stepping a symmetrical beam,/(x), in increments of 
x Q between the scans. The y-axis represents the gain of the 
beam, and the scan direction is normal to the plane of the fig- 
ure. Suppose that two neighboring beams straddle the source 
position, x s , and that the first beam is located on the origin 
of the x-axis and is the nearer of the two to x y Thus the gain 
of the first beam applied to the signal, f t , is larger than that of 
the second beam,/ 2 . 

The criterion for detection of the signal is an excess noise 
power observed in the receiver attached to the one dimen- 
sional antenna. If the beam is located far from the source posi- 
tion, the noise power distribution function is Gaussian with an 
expectation, <n>, and a variance, a 2 . If the beam is near the 
source position, the noise power distribution function is differ- 
ent in that its expectation is augmented by the product of the 
source strength and beam gain (see Fig. 2). We wish to design 
an efficient algorithm to detect this signal consistent with a 
previously set probability that it is caused by noise alone, a 0 , 
and probability of missing the signal, P Q . 


Consider two detection algorithms: 

Al. Treat each beam area independently, choosing a 
threshold Th 1 such that: 

> Pr {p y > Th l Inoise alone } (1) 

P Q > Pr (Pj < Th 1 Inoise + signal } (2) 

Since we are assuming that f x > f 2 , the weakest source which 
may be detected is one for which p x > Th x . Let s' be the 
strength of this weakest source, thus: 

= /, (3) 

A2. Linearly combine the powers observed in the two 
beams by means of a weighting function, w, and then apply a 
threshold Th 2 . Assume the same a Q and fi 0 as in algorithm Al : 

P 3 = +(1 - w)p 2 ( 4 ) 

a Q > Pr {p 3 > 77? 2 Inoise alone} (5) 

P 0 > Pr (p 3 < Th 2 Inoise + signal s 3 } (6) 

Thus a source is detected if p 3 > Th 2 . Let s" be the strength 
of the weakest source that may be detected by this algorithm, 
thus: 


In the discussion which follows we shall employ the 
notation: 


fi 

h 

n i » n 2 


S 


P 1 


P 2 
Th 


Po 

Pr{A\B } 


the gain of the (nearer) first beam area at x s 

the gain of the second beam area at x s 

power received due to background noise in 
each of the two beam areas 

time invariant strength of the source 

received power from first beam = n y +f x s 

received power from second beam = n 2 +f 2 s 

power threshold which must be exceeded to 
satisfy detection criterion 

probability of false alarm 

probability of missing signal 

probability of A gvien that B is true 


s 3 = w//+(l - w)f 2 s" (7) 


The analysis which follows derives the ratio of the sensitivi- 
ties of these two detection algorithms. Then the weight, w, is 
optimized so that s' /s" is a maximum. Since the noise is 
Gaussian, p t and p 2 are also Gaussian with probability distri- 
bution functions: 


P(P,) = ~~r~ 

y/2ir a 

if no signal is present; and 


exp 


Pip) =• 


1 


exp 


(p,.- <n>) 2 


2 a 1 


(p. - <n> - s.) 2 
2a 2 


\/2m a 

if signal s. is present (where, i = 1,2). Similarly, 

(p - <n» 2 


Pip ,) = 


1 


\Thio 


exp 


2a, 


( 8 ) 


(9) 


( 10 ) 


if no signal is present; and if signal s 3 is present, 
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1 


exp 


( 18 ) 


P(P 3 ) =-— 

\[2v a 


(P 3 -<n>~s 3 ) 2 ' 


2a l 


where we have defined: 

a\ = [w 2 +(1 - w) 2 ] a 2 


(ID 


( 12 ) 


We are now in a position to derive an expression for the 
ratio of the limiting sensitivities of the two detection algo- 
rithms, s' /s", which can then be maximized to find the optimal 
weight. The notation for the derivation is simplified if we 
define <b 0 (z): 


*„(*) = 




forz >0 


Assuming that >/ 2 , for strategy A1 we see from Eqs. (1) 
and (8) (the noise only case): 


[Th -<ri>\ 

% = Pr {P 1 >Th 1 \h =0} = j (13) 

and from Eqs. (2) and (9) (the signal present case): 


/s +<«>- Th, 

$o = Pr (Pi <Th l \s l >0} =$ o ( — 


+<«>- Th 1 j 


(14) 

On the other hand, for strategy A2 we see from Eqs. (5) 
and (10) (the noise only case): 


Th 2 -<n>\ 


% = Pr &3 > Th i I = 0} = $ 0 
and from Eqs. (6) and (11) (the signal present case): 

/s„ +<n>- Th „ 


(15) 


0 O = Pr (P 3 <^ 2 l S 3 >0 > = $ o 


(16) 


Now recall that we required a 0 to be the same for both 
algorithms. Thus Eqs. (13) and (15) are equal: 


Th 


l 




(17) 


Similarly, we require 0 O to be the same for both algorithms. 
Thus Eqs. (14) and (16) are also equal: 


(s x + <n> - Th A /s 3 + <n> - Th 2 


Combining Eqs. (17) and (18), we have Sj/s 3 = a/a 3 . We may 
use this result and Eqs. (3), (7), and (12) to find the ratio of 
the sensitivities of the two detection algorithms: 


f 2 

w +— (1 - w ) 

7 = m = ■ 1 

y/ W 2 t ( 1 - w) 2 


(19) 


We derive the optimal weighting strategy for combining the 
two responses so that the ratio of the sensitivities for the two 
detection algorithms is maximized by solving dF(w)/dw = 0: 


A 

(f 1 +/ 2 ) 


( 20 ) 


Since we have no foreknowledge of the position of the ETI 
source, we must employ our assumption that the signal 
strength is time invariant and rewrite Eq. (20) in terms of the 
received power and instantaneous noise power: 

_ fri ~ v 

W opt ~ (p x - n l ) + (p 2 - n 2 ) ^ 21) 


The noise powers contained within the observed powers are 
not available to the observer and can only be estimated by 
calculating the expectation value. A subsequent paper 
(Lokshin, in preparation) will show that the optimal weighting 
strategy cannot be used to increase the signal to noise ratio 
due to the presence of cross products between the error in the 
estimation of the weights and the difference between the two 
received powers. The optimal weighting calculation does, how- 
ever, give us a theoretical limit which we may use to judge 
practical weighting strategies. 


III. Some General Results 

We shall compare the theoretical limit to the enhancement 
of sensitivity and minimization of scalloping realized by the 
optimal weighting strategy to the sensitivity and scalloping (1) 
achieved by utilizing a single observation and (2) achieved by 
combining two observations with equal weights. 

Suppose the sky survey could detect an ETI source of 
strength, $ if the center of the beam passes directly over the 
source position. What would be the sensitivity if two adjacent 
scans^ bracket the source position? Obviously, if information 
from only one beam area is used (strategy Al) the detectable 
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signal is a function of the distance the source is offset from the 
nearer beam center (beam number 1 is assumed to be nearer): 

= S 0 lfix s ) = sjf 1 (22) 

Substituting Eq. (22) into (19) we find the expression for the 
detectable signal if one or two beam areas are used with arbi- 
trary weighting: 

Hr n _ s 0 Sqx/w 2 +(1 - w) 2 

* " fix) • F(w) = w/j + (1 - w)f 2 (23) 

This expression can now be employed to compare the peak 
sensitivity and scalloping as a function of x Q for arbitrary 
weighting strategies, wherein it is assumed that the sky survey 
is carried out by scanning along parallel tracks which are sepa- 
rated by x 0 . 

Note that if w = 1 , Eq. (23) reverts to Eq. (22), the single 
beam area case. This choice leaves the survey exposed to the 
full effect of scalloping inherent in the interscan separation. 
On the other hand, if we combine information from two beam 
areas in adjacent scans, the variation in detectable signal 
strength with x s can be decreased. In the case in which the two 
beam areas are weighted equally without giving attention to 
the strength of the signal in each we find: 


s"(x) 


/j + / 2 


(24) 


This choice of weighting function minimizes the scalloping in 
sensitivity at the price of degrading the peak sensitivity. Of 
course, application of the optimal weighting strategy achieves 



(25) 


Given the symmetrical beam shape, f{x), these expressions 
may be evaluated for all x s to determine the relative sensitivi- 
ties of the three strategies and their scalloping as a function of 
x Q . In fact, we have done so for the case in which the beam 
shape is a Gaussian of arbitrary half power beam width 
(HPBW): 


f(x) = exp 


-4 In (2) 


x \ 2 
HPBW/ 


(26) 


Figure 3 shows the relative sensitivities as a function of x s 
of the three strategies for a scan separation, x 0 = 1 HPBW. The 
x-axis is shown only over the range 0 < x s < x 0 since the 
response is symmetrical about either limit. For the case in 
which data from only a single beam area is used, scalloping is 


3 dB. For the case in which data from two neighboring beam 
areas are combined with equal weight, scalloping is reduced to 
0.34 dB. This is achieved at the cost of a loss in peak sensi- 
tivity of 1.17 dB relative to the single beam result due to 
averaging data containing almost no signal with data having the 
maximum signal to noise ratio. For the case in which the 
optimal weight has been used, scalloping is reduced to 1 .5 dB 
while peak sensitivity is enhanced 0.01 dB relative to the single 
beam result. 

Figure 4 shows the relative sensitivities as a function of x s 
for a scan separation, x 0 = 0.75 HPBW. Scalloping for the sin- 
gle beam case is reduced to 1 .7 dB whereas scalloping for the 
equal weight case is increased to 0.48 dB and the loss in peak 
sensitivity relative to the single beam result is decreased to 
0.2 dB. Note that the larger scalloping is due to a greater 
enhancement of sensitivity between the beam areas than at 
their centers. This effect was just barely noticable in the pro- 
ceeding figure. Scalloping for the optimal weight case is re- 
duced to 0.3 dB and the peak sensitivity is enhanced by 0.1 dB 
relative to the single beam result. 

Figure 5 shows the variation of scalloping as a function of 
x 0 for the three weighting strategies. For large x Q the improve- 
ment in scalloping over that resulting from a single beam is 
3 dB for equal weighting and 1.5 dB for optimal weighting. 
Note, however, that minima occur in the scalloping for the 
equal weighting and optimal weighting cases. For a scan sepa- 
ration smaller than 0.95 HPBW, the sensitivity to a source 
located between equally weighted scans increases faster than 
the sensitivity to a source located at the center of either scan 
as the scan separation is decreased. Thus the nonuniformity in 
sensitivity increases with decreasing scan separation until a 
substantial overlap is achieved. A similar phenomenon occurs 
for the optimally weighted scans after x 0 shrinks below 0.65 
HPBW. 

Figure 6 shows the variation of the peak sensitivity as a 
function of x 0 relative to that achieved using only a single 
beam. The optimal weighting strategy always achieves a better 
sensitivity than a single beam, albeit the improvement is not 
very great for x Q > 0.65 HPBW. The equal weighting strategy 
peak sensitivity has already dropped halfway to its asymptotic 
minimum at x Q 0.85 HPBW. 

The variations in scalloping and peak sensitivity shown 
in Figs. 5 and 6 are key considerations in the design of a sky 
survey strategy and will impact the survey sensitivity, given the 
constraints on antenna time and available memory and pro- 
cessing power. The sensitivity of the survey can be increased 
by (1) dwelling longer in each beam area, (2) decreasing the 
scan separation, (3) combining data from neighboring beam 
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areas, or (4) combining all three options. Each option has its 
price, however. 

In the first option, the time to complete the survey in- 
creases as the square of the ratio of new to old sensitivities, 
but the scalloping is not affected. In the second option, the 
time increases directly as the ratio of old to new scan separa- 
tions increases and the scalloping is reduced. In the third 
option, survey time is not affected and scalloping is reduced, 
but memory and processing requirements increase and peak 
sensitivity may suffer depending upon the manner of combina- 
tion chosen. The fourth option will always be chosen, but the 
mix will vary depending upon the constraints under which the 
survey will operate. 


IV. The Effect of Scan Separation in the Case 
of Uniformly Distributed Transmitters 
and Fixed Survey Time 

Plausible arguments may be advanced for any number of 
assumptions about the spatial, power, duty cycle, and trans- 
mitted frequency distribution functions for signals of ETI 
origin. Given a set of assumptions, a survey may be tailored to 
maximize the probability that it will detect a signal of that 
class. Many reasonable scenarios have been advanced in the 
literature, but the great advantage of an all sky survey lies in 
the fact that it incorporates the fewest a priori assumptions. 

We now consider the impact of scan separation and weight- 
ing strategy upon the probability of detecting an ETI signal. 
The calculation requires that some assumptions be made about 
the distribution of sources and the manner in which the survey 
will be carried out. In the light of real life constraints for a sky 
survey and our state of ignorance concerning possible ETI 
sources, we shall follow in the footsteps of Drake (Ref. 2) and 
Gulkis (Ref. 3) and assume that: 

(a) An M X N HPBW 2 area of the sky is to be surveyed in a 
fixed time, T. 

(b) The transmitters are of equal strength and are distrib- 
uted uniformly throughout space. 

Suppose that our hypothetical sources each have an effec- 
tive isotropic radiation power, P, and the boresight gain of the 
antenna is G. If a particular source is a distance, R, from earth 
and is displaced relative to the boresight by a distance, x, the 
flux seen by the receiving system in one beam area which is 
due to the source is 


s(x) = 


G • P 
4? tR 2 


•fix) 


(27) 


If the source is to be detected, this flux must be greater than 
or equal to the minimum detectable flux, s'. Thus, the source 
will be detected if it is closer than the distance R m : 

fa . p 11/2 

(28) 

Equation (19) may be substituted into (28) to find the 
expression for the maximum distance if the weighted data 
from two beam areas are used. We must keep in mind that 
f(x ) is really f l in Eq. (19), and we must explicitly show the 
dependence upon the integration time, r, since s 0 ~ r _1 1 2 : 

R m (x) ~ • fix) ' F(w) • \Zrj (29) 

The number of detectable transmitters (and thus the proba- 
bility of detecting one) is proportional to the volume of space 
observed: 

8Pr { detection} ~ 8V = R ^ (30) 

Substituting Eq. (29) into (30), we have: 

8Pr {detection} ~ r 3/4 {/(jc) • F(w)} 3/2 6S2 

Ls J (31) 

Assumption (a) may now be applied as a constraint so that 
the effect of the scan separation, x 0 , upon the probability of 
detecting a signal may be evaluated. The simplest survey 
strategy entails making scans which are N HPBWs long, step- 
ping by x 0 in the orthogonal coordinate until a distance of 
M HPBWs is covered. Thus the number of scans is equal to 
Mlx 0 . The integration time is set equal to some fraction of the 
amount of time required to scan through one HPBW, and we 
shall assume here that the fraction is unity. Thus the total time 
allowed for the survey (assuming zero time between scans) 
must be: 


T = 


M 


N • T 


(32) 


Solving for r, we may now integrate Eq. (31) over all x s to 
derive an expression for the probability of detecting a signal 
as a function of scan separation: 


/¥{detl* 0 } 


~MN 1 1/4 r GPy/r' 

!_*<> J L * J 


3/2 


f l ' 2x 0 

[w/C*) + (l - w)f(_x 0 -x) ] 

J 

1 

5 

to 

+ 

1 

to 

1 


(33) 
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Figure 7 shows the variation of a normalized Pr {detec- 
tion lx 0 ) with x 0 for the three weighting strategies. The 
normalization is chosen so that the maximum relative proba- 
bility for the optimal weighting strategy is unity. As x Q be- 
comes larger, the scans can be slowed down and still allow the 
survey to cover the same area of sky in the given time limit. 
However, the scalloping in sensitivity increases with x 0 and 
degrades the probability of detection. For very small values of 
x 0 , the optimal weighting and equal weighting strategy 
improve the probability over that achieved by the single beam 
strategy by the expected ratio of 2 3 / 4 due to the effective 
doubling of integration time on source. As the scan separation 
increases the single beam area strategy approaches the optimal 
result and the eqnal weighting strategy falls off to the ex- 
pected ratio of 1/2 due to doubling of the noise. 

The relative probabilities of detection for the single beam 
strategy and the equal weighting strategy are equal for a scan 
separation of about 0.8 HPBW, and are degraded relative to 
the optimal result by approximately 15%. The probability of 


detection achieved by the equal weighting strategy peaks at a 
scan separation of 0.6 HPBW, and it is degraded from the 
optimal peak by 8%. The probability of detection achieved 
by using a single beam area peaks at a scan separation of 1.3 
HPBW, and it is degraded from the optimal peak by 5%. It is 
clear that the scalloping allowable in a survey will depend 
upon the assumptions of the designers of the search strategy. 


V. Suggestions for Further Analysis 

It is possible to extend the foregoing analysis to cover the 
general case of combining N beam areas. Of more immediate 
concern, however, is an extension to the case of a non- 
Gaussian noise statistic. The contemplated high speed all sky 
survey will operate in this domain due to the small number 
(4 < n < 1 00) of independent samples which will be combined 
before thresholding. A series of papers is in preparation which 
will cover this topic. 
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POSITION, HPBW 


Fig. 1. One dimensional model of survey geometry 



Fig. 3. Relative sensitivity as a function of source position for a 
scan separation of 1.0 HPBW 



<n) Th <n)+f^s 


RECEIVED POWER 



Fig. 2. Probability distribution functions of received power due to 
noise alone and due to noise and source 


Fig. 4. Relative sensitivity as a function of source position for a 
scan separation of 0.75 HPBW 
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SCAN SEPARATION, HPBW 


. Scalloping as a function of scan separation 



Fig. 6. Relative peak sensitivity as a function of scan separation 



Fig. 7. Normalized relative probability to detect ETI as a function of 
scan separation 
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